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ABSTRACT 

The  flowfield  about  a  Vertically-Launched  Surface-to-Air  Missile  model 
at  an  angle  of  attack  of  50°  and  a  Reynolds  number  of  l.lxlO5  was  investi- 
gated in  a  low  speed  wind  tunnel  at  the  Naval  Postgraduate  School.  The  goal 
of  this  thesis  is  to  determine  the  location  and  intensity  of  the  asymmetric 
vortices  in  the  wake  of  the  VLSAM  model  and  to  display  these  vortices  by 
velocity  mapping  and  pressure  contours.  The  two  model  configurations  tested 
were  for  a  cruciform  missile  with  wings  and  tails;  one  at  0°  roll  angle  ("plus" 
aspect)  and  the  other  at  a  45°  roll  angle  ("cross"  aspect).  Two  flowfield 
conditions  were  treated:  the  nominal  ambient  wind  tunnel  condition  and  a 
condition  with  a  grid-generated  turbulence  of  length  scale  1.08  inches  and 
1.88%  turbulence  intensity.  The  turbulence  length  scale  is  61.7%  of  the  model 
diameter  and  4.7%  of  the  model  length.  The  following  conclusions  were 
reached:  1)  An  increase  in  turbulence  intensity  tended  to  reduce  the  strength 
of  the  asymmetric  nose-generated  vortices;  2)  the  two  asymmetric  vortices 
remained  in  approximately  the  same  position  for  an  increase  in  turbulence;  3) 
"cross"  aspect  vortices  were  more  diffused,  slightly  larger  and  centered 
further  away  from  the  model  surface  than  those  of  the  "plus"  aspect  body 
configuration,  which  correlates  with  the  differences  in  induced  side  forces  for 
these  configurations  observed  by  Rabang;  4)  the  top  vortex  of  the  two 
asymmetric  vortices  was  closer  to  the  model  surface  and  appeared  to  be 
stronger  for  both  configurations;  and  5)  the  addition  of  wings  and  tails  did  not 
greatly  alter  the  vortex  pattern  around  the  nose  of  the  missile  model. 
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I.  INTRODUCTION 

A.  BACKGROUND 

The  development  of  the  Vertical  Launch  Surface-to- Air  Missile  (VLSAM) 
has  provided  greater  weapon  system  reliability,  availability  and  flexibility 
over  its  predecessors.  Earlier  systems,  such  as  trainable  launchers  and  box 
launchers,  were  rather  cumbersome  and  required  excessive  deck  space,  thus 
limiting  the  number  of  missiles  per  ship.  Additionally,  these  systems  were 
slow  to  reload  since  trainable  launchers  required  cycling  of  rounds  to  get  the 
desired  one  in  position  for  launch,  while  box  launchers  were  reloaded  by 
hand.  The  vertical  launcher  provides  for  quick  access  to  any  round  without 
cycling,  rapid  reloading  and,  due  to  its  design  and  the  fact  that  the  missile 
blast  is  kept  in  a  concentrated  area  near  the  launcher,  it  requires  less  deck 
space.  [Refs.  1,  2  and  3] 

The  VLSAM's  trajectory  allows  it  to  point  to  its  target  after  launch  and 
subsequently  guide  itself  to  the  correct  heading.  Its  aerodynamic  characteris- 
tics may  significantly  change  as  it  operates  from  subsonic  to  supersonic 
speeds  during  the  launch,  midcourse  and  terminal  phases  of  its  flight.  [Ref.  2] 
The  varying  flight  control  requirements  and  the  airflow  about  the  missile 
during  these  phases  provide  for  different  flight  regimes. 

When  launched,  the  VLSAM  enters  an  open  ocean  environment  and  is 
subject  to  potentially  significant  cross-winds,  the  result  of  which  is  a  missile 
flying  at  relatively  low  velocities  at  a  high  angle  of  attack.  [Refs.  1  and  2]  In 
particular,  an  example  of  this  low  velocity/high  angle  of  attack  condition  is 
the  Standard  Missile  2-Block  4,  or  AEGIS  Extended  Range  version,  which 


has  an  added  sustainer  section  and  leaves  its  launcher  at  a  much  lower  speed 
than  the  unboosted  SM-2  version.  A  missile  flying  at  relatively  low  velocities 
at  a  high  angle  of  attack  may  bring  with  it  the  formation  of  asymmetric 
vortices  about  the  missile  nose  and  afterbody  which  can  generate  unpre- 
dictable and  undesired  side  forces,  thereby  affecting  overall  missile  flight  sta- 
bility. The  airflow  around  a  ship's  superstructure,  the  ocean  surface  condi- 
tions and  the  atmospheric  boundary  layer  conditions  during  launch  may  also 
provide  turbulent  flow  over  the  missile.  [Ref.  2]  The  advent  of  the  VLSAM 
and  the  desire  to  have  highly  maneuverable  supersonic  missiles  have 
increased  the  need  for  further  studies  in  high  angle  of  attack  research  under 
various  flowfield  conditions. 

This  thesis  continues  experiments  and  research  conducted  to  date  at  the 
Naval  Postgraduate  School  (NPS)  to  understand  the  effects  of  turbulence  on 
the  VLSAM  aerodynamic  characteristics.  Previously,  Roane  [Ref.  1]  developed 
a  system  model  to  measure  flowfield  turbulence  through  the  use  of  four  dif- 
ferent grids  which  generate  varying  turbulence  levels  in  the  NPS  low  speed 
wind  tunnel.  Rabang  [Ref.  2]  studied  the  effects  of  this  turbulence  on  the 
asymmetric  vortex  forces  on  the  missile  model.  Lung  [Ref.  3]  determined  the 
location  and  intensity  of  the  asymmetric  vortices  in  the  wake  of  the  model  by 
experimental  flowfield  measurements  about  a  body-only  missile  configuration 
with  and  without  freestream  turbulence.  Similarly,  the  goal  of  this  study  is  to 
determine  vortex  locations  and  intensities,  through  a  series  of  wind  tunnel 
experiments,  for  two  missile  configurations,  both  with  and  without  turbu- 
lence. The  configurations  considered  are  for  a  cruciform  missile  with  wings 


and  tails;  one  at  0°  roll  angle  ("plus"  aspect)  and  the  other  at  a  45°  role  angle 
("cross"  aspect). 

B.     HIGH  ANGLE  OF  ATTACK  AERODYNAMICS 

In  high  angle  of  attack  aerodynamics,  flow  separation  from  the  body, 
wing  and  tail  surfaces  is  important  due  to  the  strong  normal  and  side  forces 
which  may  be  generated.  Major  factors  which  influence  flow  separation  are 
nose  shape,  angle  of  attack,  crossflow  Reynolds  number  and  nose  fineness 
ratio.  Other  factors  include  roll  angle  and  roll  rate,  free  stream  turbulence, 
surface  roughness,  acoustic  environment  and  VLSAM  model  vibrations 
[Ref.  4]. 

1.      Aerodynamic  Regimes 

As  the  angle  of  attack,  °s  of  a  slender  body  of  revolution  ranges 
from  0°  to  90°,  there  are  at  least  four  distinct  airflow  regimes  through  which 
the  missile  body  transitions.  [Ref.  4] 

(1)  Regime  I  (0°<  <*  <  5°):  At  very  low  angles  of  attack,  the  axial  flow  domi- 
nates and  there  is  no  discernable  boundary  layer  separation  (flow  is 
attached). 

(2)  Regime  II  (5°<  «=  <  20°):  At  intermediate  angles  of  attack,  boundary 
layer  separation  occurs  on  the  lee  side  of  the  body.  This  becomes  a  free 
shear  layer  which  rolls  up  into  a  symmetric  vortex  pair  that  is  steady 
with  time.  No  side  forces  or  yawing  moments  are  induced. 

(3)  Regime  III  (20°<  «*=  <  60°):  At  high  angles  of  attack,  crossflow  effects 
start  to  dominate  and  vortices  are  now  shed  asymmetrically.  These 
vortices  induce  side  forces  (out-of-plane  forces)  and  yawing  moments. 
The  more  asymmetric  the  vortex,  the  greater  the  side  force  magnitude. 
There  are  some  flow  instabilities  toward  the  higher  end  of  this  angle  of 
attack  range. 


(4)  Regime  IV  (60°<  «  <  90°):  At  very  high  angles  of  attack,  the  crossflow 
completely  dominates  and  flow  separation  is  unsteady.  The  Reynolds 
number,  Mach  number  and  geometry  determine  whether  the  boundary 
layer  is  shed  as  a  von  Karman  vortex  street  or  a  random  wake-like 
flow.  [Refs.  1,  2  and  3] 

Figure  1  shows  sketches  of  all  four  vortex-generation  regimes. 
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Figure  1.  Airflow  Regimes  [Ref.  2] 


2.      Asymmetric  Vortex  Theory 

The  principal  cause  of  the  formation  of  asymmetric  vortices  is  still 
not  completely  understood  and  may  be  attributed  to  many  factors.  One  idea  is 
that  boundary-layer-induced  asymmetry  in  the  location  of  flow  separation 
causes  the  vortex  flowfield  to  become  asymmetric.  These  boundary  layer 
asymmetries  may  result  from  transition  and  separation  differences  on  each 
side  of  the  missile  body.  Another  proposition  is  that  a  hydrodynamic  (inviscid) 
instability  in  the  pair  of  initially  symmetric  vortices  causes  the  asymmetry. 
[Refs.  5,  6  and  7]  These  vortices,  which  increase  in  strength  with  angle  of 
attack,  interact  with  the  surrounding  potential  flowfield  to  provide  the 
asymmetric  configuration.  A  vortex-switching  phenomenon  has  also  been 
observed  in  which  the  vortex  pattern  rapidly  switches  from  an  almost  sym- 
metric to  a  highly  asymmetric  configuration,  which  may  possibly  relate  to  a 
second  inviscid  solution  in  the  leeward  flowfield.  [Refs.  8  and  9]  At  any  rate, 
even  though  their  major  cause  has  not  been  determined,  the  behavior  of 
asymmetric  vortices  has  been  well  documented  for  a  large  number  of  models 
and  shapes. 

Nose-generated  asymmetric  vortices  appear  in  the  flowfield  around  an 
ogive-nosed,  slender,  cylindrical  body.  The  vortex  formation  occurs  along  the 
entire  body  length  and  induces  a  significant  side  force  on  the  body.  Figure  2 
shows  this  vortex  flow  along  the  length  of  an  unyawed,  slender  nose  cylinder. 


FREE   STREAM 


VELOCITY 


At  CLE  OF  ATTACK 


Figure  2.  Vortex  Flow  About  a  Slender  Nose  Cylinder  [Ref.  1] 


3.      Two-Dimensional  Crossflow 

Airflow  over  a  slender  body  can  be  divided  into  normal  and  axial 
components.  Axial  flow  is  along  the  missile  body  length  while  crossflow  is 
essentially  a  two-dimensional  flow  normal  to  a  cylinder.  The  crossflow 
analogy  provides  information  for  cylinder  lift  and  drag  which  act  in  the 
crossflow  direction.  Depending  on  the  type  of  flow  separation  on  either  side  of 
the  cylinder,  side  forces  (at  right  angles  to  the  crossflow)  may  exist.  [Ref.  2] 

Boundary  layer  transition  and  sej^ration  mechanisms  may  provide 
an  explanation  for  flow  separation  and  subsequent  asymmetric  vortex  genera- 
tion. The  primary  factor  which  influences  the  separation  location  of  the 


boundary  layer  is  the  crossflow  Reynolds  number.  Other  factors  include  sur- 
face roughness  and  turbulence. 

Flow  around  a  cylinder  in  incompressible  flow  can  be  classified  into 
four  regions,  represented  by  differing  flow  separation  and  drag  behavior,  as 
depicted  in  Figure  3.  [Refs.  10,  11  and  21  In  the  subcritical  range,  the  bound- 
ary layer  is  laminar  and  flow  separation  occurs  close  to  the  lateral  meridian 
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Figure  3.  Two-Dimensional  Crossflow  About  a  Cylinder  [Ref.  2] 


where  the  angle  (<{>)  from  the  crossflow  direction  varies  from  80°  to  90°. 
[Ref.  12]  In  the  critical  Reynolds  number  region,  a  drag  bucket  is  produced. 
The  laminar  boundary  layer  separates  from  the  body  at  <j)  =  90°  and  reat- 
taches as  a  turbulent  boundary  layer  which  is  more  energetic.  Separation  is 
delayed  to  $  ~  140°,  resulting  in  a  reduction  of  the  drag.  [Ref.  3]  A  laminar 
bubble  is  formed  between  the  laminar  separation  and  the  turbulent  reat- 
tachment. At  this  point,  the  flow  separation  may  easily  fluctuate  from  critical 
to  subcritical  for  small  changes  in  Reynolds  number.  From  Figure  4  [Ref.  10], 
when  critical  separation  exists  on  one  side  of  the  body  while  subcritical  sepa- 
ration is  on  the  other,  a  large  difference  in  <})  is  possible.  Therefore,  vortices 

will  be  at  maximum  asymmetry  and  the  side  force  at  the  highest  magnitude. 
Since  maximum  side  force  occurs  at  the  critical  Reynolds  number,  it  is  a 
noteworthy  parameter.  [Ref.  13] 

As  the  Reynolds  number  is  further  increased,  turbulent  separation 
moves  forward  to  $  <  140°  and  the  laminar  bubble  no  longer  exists.  At 
<j)  >  140°,  asymmetric  vortices  are  ineffective  at  producing  a  significant  side 
force,  thus  the  sudden  decrease  in  magnitude  (Cy/Cn)  as  shown  in  Figure  4. 
For  supercritical  and  transcritical  Reynolds  numbers,  the  laminar  transition 
point  moves  towards  §  ~  0°  and  turbulent  separation  occurs  at  0  =  100°.  The 
asymmetric  transcritical  separation  point  moves  towards  the  lateral  merid- 
ian, where  the  vortices  once  again  produce  a  significant  side  force.  [Refs.  5 
and  2]  The  Reynolds  number  provides  the  greatest  influence  on  the  normal 
force  and  drag  characteristics,  particularly  within  the  critical  range  where 
the  maximum  normalized  side  force  and  maximum  vortex  asymmetry  occur. 
[Ref.  12] 
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Figure  4.  Side  Force  to  Normal  Force  Ratio  [Ref.  5] 


Another  study  by  Lamont  [Ref.  14]  describes  a  different  effect  of 
Reynolds  number  on  the  maximum  side  force  as  illustrated  by  Figure  5, 
where  the  side  force  at  an  angle  of  attack  of  55°  is  plotted.  The  maximum  side 
force  falls  from  a  high  value  at  laminar  separation  to  almost  zero  in  the 
middle  of  the  transition  region  before  climbing  again  to  higher  levels  at  fully 
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Figure  5.  Effect  of  Reynold's  Number  of  Maximum 
Side  Force  at  -  =  55°  [Ref.  14] 


turbulent  separation.  Thus,  there  appears  to  be  two  different  mechanisms  for 
producing  asymmetric  flow  and,  hence,  a  side  force  on  an  ogive-cylinder.  One 
mechanism  operates  in  both  the  laminar  and  the  fully  turbulent  separation 
regimes,  in  which  the  side  force  results  from  asymmetric  vortex  patterns  in 
the  wake  of  a  body.  The  other  mechanism  occurs  only  in  the  transitional 
separation  regime.  Here,  the  Reynolds  number  at  which  the  near-zero  side 
forces  were  recorded,  is  the  same  range  of  Reynolds  number  in  which  the 
minimum  drag  coefficient  on  a  2— D  cylinder  occurs  and  in  which  no  coherent 
vortex  shedding  can  be  detected. 
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4.      Three-Dimensional  Vortices 

The  missile  nose  geometry  is  an  important  factor  in  vortex  genera- 
tion and  disposition  since  vortices  shed  at  the  nose  tend  to  dominate  other 
vortices  along  the  body  length.  [Ref.  15,  8  and  16] 

Nose-generated  vortices  are  sensitive  to  the  nose  roll  angle  due  to 
surface  imperfections  and  nose  geometric  deviations  [Ref.  17].  Rabang  varied 
the  roll  angle  in  45°  increments  and  investigated  the  resulting  side  force  coef- 
ficients, shown  in  Figure  6.  The  vortex  system  generated  by  the  nose  domi- 
nates afterbody  vortices  for  body  configurations  with  and  without  wings 
regardless  of  the  turbulence  conditions  [Refs.  2,  18,  and  19]. 
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Nose  geometry  may  be  pointed  or  blunt  for  cones  and  ogives.  For 
pointed  noses,  angle  of  attack  for  the  onset  of  asymmetric  vortices  (^av)  is  a 
function  of  the  semi-vertex  angle  (9a).  Asymmetric  vortices  start  at  the  nose 
and  are  rapidly  shed,  yielding  unsteady  side  forces.  [Ref.  3]  At  all  Mach 
numbers,  asymmetric  vortex  shedding  starts  when  «av  is  greater  than  8a- 
For  a  conical  nose: 

^a  ~  "'av/  2  (approximation)  (1) 

For  a  tangent  ogive  nose: 


0A  =  tan"1 
A 


lJvd-0.25 

N 


(2) 


where  In  is  the  nose  length  and  d  is  the  base  diameter,  or 

©A  =  Wd  (approximation)  (3) 

Nose  fineness  ratio  also  affects  the  asymmetric  vortex  induced  side 
forces  in  that  as  this  ratio  increases,  the  side  force  also  increases.  With  an 
increasing  ratio,  both  the  nose  apex  angle  (0a)  and  the  angle  of  attack  for  the 
onset  of  asymmetric  vortices  (^av)  will  decrease,  making  the  missile  more 
susceptible  to  induced  side  forces  at  lower  angles  of  attack.  [Refs.  15  and  20] 
Decreasing  nose  fineness  ratio  has  been  found  to  be  more  beneficial  in  reduc- 
ing side  forces  than  blunting  the  nose  [Ref.  15].  Side  force  decreases  with  an 
increasing  Mach  number  to  nearly  zero  at  supersonic  speeds  [Refs.  18  and 
21]. 
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C.     TURBULENCE 

Turbulence  denotes  the  presence  of  random,  short  duration  variations  in 
a  flowfleld  with  a  given  mean  velocity.  When  calculating  turbulence  effects  on 
a  body  in  the  flowfleld,  a  comparison  between  the  scale  of  the  body  and  that  of 
the  turbulence  must  be  made.  The  energy  in  the  turbulence  flowfleld  should 
also  be  considered.  [Ref.  1] 

Turbulence  intensity,  Tu,  is  the  measure  of  the  relative  magnitude  of 

velocity  fluctuations  in  the  flowfleld.  For  a  horizontal  flowfleld  or  crosswind, 
it  is  the  ratio  of  the  root-mean-square  (rms)  velocity  fluctuation,  u',  to  the 
mean  velocity  component  in  the  flowfleld,  U«>. 

Tu  =  u'/Uoo  (4) 

Turbulence  length  scales  describe  the  time-averaged  measure  of  the  size 
of  the  constantly  changing  fluid  disturbance  eddies.  An  increase  in  the  spatial 
length  of  the  turbulence  corresponds  to  an  increase  in  the  time  the  body  is 
exposed  to  the  fluctuation.  Large  and  small  scale  turbulence  length  scales  are 
both  found  in  a  flowfleld. 

From  a  single  source,  the  cascade  effect  produces  turbulence  eddies  of 
different  length  scales.  This  "cascade"  effect  is  caused  by  a  strain  in  one  direc- 
tion (x,  y,  z  plane)  which  affects  the  orthogonal  components  due  to  the  conser- 
vation of  angular  momentum.  For  example,  an  increase  in  the  x  and  y  velocity 
components  of  a  vortex  rotating  in  the  x-y  plane  will  have  an  effect  on  the  ve- 
locity and  length  scales  of  the  y-z  components.  [Ref.  1]  Cascading  continues 
until  the  smaller  eddies  disappear  due  to  the  viscosity.  As  turbulence 
decreases,  the  energy  transfer  decreases  and  the  individual  intensities  of 
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each  eddy  will  decrease  at  a  faster  rate.  [Ref.  22]  Thus,  the  larger  scale  turbu- 
lence predominates. 

The  length  scale  to  body  size  ratio  may  determine  the  manner  in  which 
turbulence  affects  the  VLSAM  flowfield.  It  may  be  compared  to  missile 
length,  Lu:Li,  or  missile  diameter  Lu:Ld.  [Refs.  1  and  2] 

For  length  scales  much  greater  than  the  body,  Lu»Li,  the  effect  is  like  a 
steady-state  flowfield,  where  deviations  in  speed  and  direction  would  be  of 
long  duration.  The  flowfield  effects  on  vortex  development  are  mainly  domi- 
nated by  the  same  factors  and  conditions  as  for  a  two-dimensional  cylinder. 

In  contrast,  unwanted  rolling,  pitching  and  yawing  motion  of  the  body  is 
primarily  caused  when  the  turbulence  length  scale  is  comparable  to  the  body 
length,  Lu~L].  [Ref.  23]  The  flowfield  is  distinctly  non-steady  for  this  case. 

When  the  length  scale  is  of  a  dimension  much  smaller  than  the  body, 
most  significantly,  when  it  is  smaller  than  the  missile  diameter,  Lu«Ld,  it 
has  a  magnitude  comparable  to  the  boundary  layer  thickness  on  the  missile 
surface.  Thus,  boundary  layer  development  and  flow  separation  over  the  body 
may  be  affected  by  the  presence  of  small  scale  turbulence.  An  increase  in  tur- 
bulence intensity  with  a  length  scale  on  the  order  of  the  boundary-layer  scale 
tends  to  reduce  the  magnitude  of  induced  side  forces.  [Refs.  2,  24,  and  25]  A 
goal  of  the  current  investigation  is  to  determine  the  effect  of  vortex- scale 
turbulence  on  the  asymmetric  vortices  and  resulting  induced  side  forces. 

D.     LIFTING  SURFACE  EFFECTS 

The  complete  vortex  structure  of  the  missile  is  a  net  result  of  the  indi- 
vidual contributions  from  the  body,  wings,  strakes  and  tails.  In  general,  mis- 
siles use  low  aspect  ratio  wings  (when  compared  with  aircraft).  Since  some 
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missiles  have  wing  spans  that  approach  body  diameter,  it  is  important  to 
consider  the  joint  effects  from  a  wing-body  combination.  Nose  vortices  dictate 
flow  behavior  over  a  missile  body  at  high  angles  of  attack  and,  consequently, 
these  vortices  may  also  be  felt  by  the  wings.  Nose  and  body  vortices  move 
away  from  a  missile  body  without  wings  but,  when  wings  are  added,  they 
move  closer  to  the  body.  This  result  is  comparable  to  increasing  the  effective 
angle  of  attack  causing  unsteady  asymmetric  vortices.  For  wings  with  low 
aspect  ratio,  a  major  portion  of  the  lift  produced  by  the  wing  will  be  a  result  of 
vortex  lift.  The  net  effect  of  the  wing-body  combination  appears  to  be  a  reduc- 
tion in  the  effective  angle  of  attack  for  the  onset  of  asymmetric  vortices  and 
side  forces.  [Refs.  2,  3,  and  26] 

Vortex  lift  effects  are  improved  by  incorporating  strakes  with  low  aspect 
ratio  wings.  The  strakes  produce  additional  strong  vortices.  Some  researchers 
have  found  that  placing  long  strakes  on  a  missile  would  induce  interference 
with  the  crossflow  component  around  the  body,  thus  decreasing  the  effect  of 
the  forces  and  moments  generated  by  asymmetric  body  vortices.  [Ref.  6] 
Rabang  has  shown  that  the  addition  of  typical  VLSAM  wings  and  strakes 
tend  to  preserve  the  induced  side  force  for  all  levels  of  turbulence  intensities 
and  length  scales.  [Ref.  2] 

The  addition  of  tails  has  a  minor  influence  on  forebody  flowfields  and 
maximum  side  forces,  in  particular,  at  low  angles  of  attack.  At  higher  angles 
of  attack,  nose  and  wing  vortices  may  have  a  slight  effect  on  the  tailflow, 
depending  on  the  wing  placement  and  afterbody  length  between  the  wings 
and  tails. 
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E.     VLSAM  LAUNCH  ENVIRONMENT 
1.      Marine  Environment 

Turbulence  conditions  which  exist  within  the  atmospheric  boundary 
layer  (ABL)  may  significantly  impact  the  VLSAM.  The  ABL  is  the  lowest  por- 
tion of  the  atmosphere  and  is  formed  by  its  interaction  with  the  surface  over 
which  it  flows.  Turbulence  in  this  layer  is  the  result  of  the  transfer  of  heat, 
momentum  and  mass. 

The  surface  layer,  the  lowest  segment  of  the  ABL,  can  vary  in 
height  from  5  to  200  meters  but  is  typically  on  the  order  of  50  meters.  It  is 
also  characterized  by  mechanically  produced,  small-scale  turbulence  resulting 
from  surface  roughness  or  friction  from  waves  on  the  ocean  surface.  This 
small-scale  turbulence  is  larger  than  the  missile  length.  This  region  is 
described  by  variations  in  wind  speed,  nearly  vertical  heat  and  mass  fluxes, 
and  other  meteorological  fluctuations  with  height.  [Refs.  27  and  23]  Further- 
more, the  majority  of  the  flow  in  the  surface  layer  itself  can  be  considered 
horizontally  homogeneous  and  two-dimensional  [Ref.  23]. 

A  measure  of  the  roughness  of  the  surface  is  called  the  roughness 
length,  Z0,  which  is  a  function  of  the  mean  wind  velocity  at  various  heights 

above  the  surface.  By  combining  the  roughness  length  with  the  elevation  and 
wind  speed,  both  the  turbulence  intensity  and  length  scale  can  be  empirically 
determined.  [Refs.  2,  3  and  28]  Typical  open  ocean  surface  roughness  lengths 
are  in  the  range  of  0.001<Z0<0.01,  with  Z0  in  meters.  For  a  10  meter  elevation 
and  a  mean  wind  speed  of  25  m/sec  in  a  neutral  atmosphere,  turbulence 
intensities  may  range  from  13  to  17  percent.  [Ref.  28]  The  longitudinal 
turbulence  length  scale  would  then  range  from  262<LU<295  feet  (80<Lu<  90 
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meters).  Therefore,  for  a  typical  missile  with  a  1.1  foot  diameter,  the 
turbulence  length  scale  to  missile  diameter  length  scale  ratio  is  Lu:Ld~280:l. 
[Ref.  28]  This  represents  a  length  scale  very  much  larger  than  a  conventional 
missile  length  and,  therefore,  would  have  little  effect  on  its  boundary  layer 
development.  However,  the  cascade  effect  from  large  scale  turbulence  and 
from  crosswind  interaction  with  a  ship's  superstructure  allows  length  scales, 
initially  much  larger  (85  meters)  than  the  dimension  of  a  missile,  to  decrease 
(cascade)  to  scales  where  they  could  affect  the  missile  boundary  layer 
development  and  the  development  of  vortices  from  the  missile  nose.  The 
actual  amount  of  such  small-scale  turbulence  present  in  the  marine 
atmosphere  is  largely  unknown,  however.  [Ref.  23] 
2.      Launch  and  Crosswind  Velocities 

A  typical  VLSAM  at  launch  (vertical  velocity  =  164  ft/sec)  is  still 
well  within  the  surface  layer  environment  and  is  subject  to  both  crosswind 
and  turbulence  effects.  [Ref.  29]  Crosswind  velocity  depends  on  both  the 
ambient  wind  speed  and  the  speed  of  the  launch  platform.  A  ship  speed  of  20 
knots  with  a  head  wind  of  20m/sec,  combined  with  the  VLSAM  launch 
velocity,  results  in  an  effective  angle  of  attack  of  31°  at  191  ft/sec,  which 
places  the  missile  in  the  asymmetric  vortex  region  (Regime  III)  almost  imme- 
diately after  clearing  the  exit  plane  of  the  launcher,  0.2  seconds  after  launch. 
[Refs.  1  and  2] 

When  the  missile  first  leaves  its  launcher  it  will  experience  an  even 
greater  effective  angle  of  attack  due  to  its  slower  velocity.  A  ship's  hull  and 
superstructure  can  dictate  changing  flow  fields  and  turbulence  at  this  initial 
launch  altitude.  The  ship  airwake  may  increase  crosswind  velocities  and 
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cause  significant  crosswind  gradients  in  the  flowfield,  thus  increasing  turbu- 
lence intensities  while  decreasing  turbulence  length  scales.  [Ref.  2] 

Later  in  its  launch  profile,  when  the  VLSAM  pitches  over  towards  a 
target,  it  may  reach  effective  angles  of  attack  of  up  to  50°.  [Refs.  2  and  30] 
Thus,  during  the  launch  phase  there  exists  a  definite  possibility  of  asymmet- 
ric vortex  induced  side  forces  on  a  VLSAM. 

3.      Additional  Launch  Considerations 

During  launch,  a  missile  is  influenced  by  many  factors  which  are 
dependent  on  missile  design,  ship's  orientation,  launcher  mechanics  and 
missile  flight  control.  Shipboard  roll,  pitch  and  yaw  are  directly  transmitted 
to  the  launch  platform  and  must  be  taken  into  account.  Inherent  factors,  such 
as  the  plume  (jet)  effect  of  the  missile's  engine  and  blast  effects  of  the  vented 
exhaust  gases,  can  also  affect  VLSAM  aerodynamics.  Exhaust  gases  can 
impinge  directly  on  the  missile  surfaces  or  they  can  impact  the  flowfield  in 
which  the  missile  is  launched,  especially  if  the  gases  are  vented  upward  into 
the  vicinity  of  the  accelerated  missile.  [Ref.  1]  The  manner  in  which  control 
systems  respond  to  missile  orientation  changes  is  another  factor.  Should  the 
missile  change  its  flight  attitude,  the  flowfield  around  it  will  also  be  altered. 
Obviously,  the  considerations  discussed  in  this  section  are  not  all-encompass- 
ing. There  are  many  other  factors  that  affect  missile  flight  behavior  during 
launch,  but  they  are  beyond  the  scope  of  this  thesis  and  will  not  be  included. 
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II.  EXPERIMENT  AND  PROCEDURES 

A.     PURPOSE 

In  this  study,  the  location  and  intensity  of  asymmetric  vortices  in  the 
wake  of  the  VLSAM  model  were  determined  for  varying  levels  of  turbulence. 
This  turbulence  was  generated  by  the  placement  of  a  series  of  grids  in  the 
wind  tunnel.  The  vortices  were  displayed  by  velocity  mapping  and  pressure 
contours.  To  accomplish  this,  wind  tunnel  flowfield  pressure  measurements 
were  taken  for  a  specified  survey  grid  using  a  scanivalve/probe  data  acquisi- 
tion system. 

Figure  7  [Ref.  3]  shows  the  planar  survey  grid,  the  x-y  plane,  which  was 
perpendicular  to  the  freestream  velocity  and  located  10.5  inches  downstream 
from  the  missile  models  nose.  The  model  body  was  centered  on  the  y  dimen- 
sion. For  the  actual  data  acquisition  runs,  23  points  were  measured  in  the  y 
direction.  There  were  11  points  above  and  below  the  model  centerline,  with 
point  12  directly  at  the  centerline.  Along  the  x  axis,  13  points  were  measured, 
the  x-y  dimension  for  the  experiment  was  3x5.5  inches,  with  a  step  distance 
(increment)  of  0.25  inch  (22x0.25=5.5  and  12x0.25=3  inches).  This  dimension 
covered  the  main  portion  of  asymmetric  vortices. 

Pressure  measurements  were  obtained  by  the  5-hole  probe  throughout 
the  survey  plane.  The  data  from  the  pressure  probe  was  reduced  through  the 
use  of  computer  programs  to  obtain  isobars  of  total  pressure  coefficient  and 
static  pressure  coefficient,  and  to  map  the  crossflow  velocity  vectors.  These 
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Figure  7.  The  Planar  Survey  Grid  [Ref.  3] 
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results  were  correlated  with  the  force  measurements  of  Rabang  and  with  the 
previous  experiments  by  Lung  to  provide  a  greater  understanding  of  the 
vortex  flowfield.  The  following  sections  further  discuss  the  equipment  and 
software  used,  and  the  experimental  procedures  followed. 

B.     APPARATUS 

Information  about  the  construction,  specifications  and  configurations  of 
the  major  pieces  of  equipment  used  in  this  study  is  described  in  this  section. 

1.      Wind  Tunnel 

The  low-speed,  single  return,  horizontal-flow  tunnel  located  in 
Halligan  Hall  at  NPS  was  utilized.  (Figure  8,  [Ref.  31])  .  It  is  powered  by  a 
100  horsepower  electric  motor  coupled  to  a  three-blade  variable-pitch  fan  via 
a  4-speed  Dodge  truck  transmission.  Aft  of  the  fan  blades  are  a  set  of  stator 
blades  which  help  straighten  flow.  Two  fine  wire  mesh  screens  located  up- 
stream of  the  settling  chamber  plus  turning  vanes  at  all  four  corners  reduce 
turbulence.  A  heavy  wire  screen  behind  the  test  section  prevents  foreign 
object  damage  to  the  fan  blades.  [Refs.  1,  2  and  31]  The  tunnel  is  64  feet  long 
and  ranges  from  21.5  to  25.5  feet  wide. 

The  wind  tunnel  test  section  measures  45  inches  by  32  inches.  The 
walls  diverge  slightly  to  prevent  reduction  in  freestream  pressure  due  to 
boundary  layer  growth.  The  settling  chamber  area  to  test  section  area 
contraction  ratio  is  approximately  10:1.  Corner  fillets,  which  house  the  light- 
ing, reduce  the  section  area  from  10  ft2  to  9.88  ft2.  Fillets  are  found  at  wall 
intersections  throughout  the  tunnel  to  help  reduce  boundary  layer  effects. 
[Ref.  31] 
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Figure  8.  Naval  Postgraduate  School  Wind  Tunnel  [Ref.  31] 


A  reflection  plane  installed  in  the  test  section  reduces  the  available 
height  to  28  inches.  A  flush-mounted  turntable  allows  for  changes  in  model 
pitch  angle  or  angle  of  attack  via  a  remotely  controlled  electric  motor  beneath 
the  tunnel.  Since  the  test  section  operates  at  atmospheric  pressure,  breather 
slots  are  installed  around  the  tunnel  perimeter  to  replenish  air  lost  through 
leaks  and  to  ensure  a  uniform  test  section  pressure.  The  tunnel  was  designed 
to  provide  test  section  velocities  of  up  to  290  ft/sec.  [Ref.  31] 
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Wind  tunnel  temperature  is  measured  by  a  dial  thermometer 
extending  into  the  settling  chamber.  Dynamic  pressure  is  measured  by  the 
static  pressure  difference  between  the  test  section  and  the  settling  chamber 
using  a  water  filled  manometer.  The  static  pressure  is  measured  by  four  pres- 
sure taps  located  upstream  from  the  test  section  to  avoid  interference  from 
the  model.  These  taps  are  connected  via  a  common  manifold  prior  to  feeding 
into  the  manometer.  Pressure  differences  are  measured  in  centimeters  of 
water.  Equation  (5)  is  used  to  convert  to  the  actual  wind  tunnel  velocity. 
[Ref.  3] 


U    = 

m 


(2X2.0475)  (P^O) 


(K)(p) 


(5) 


where: 


Um  =  measured  velocity  (ft/sec) 

2.0475  =  conversion  factor 

Pcm  H2O  =  manometer  reading 

K  =  calibration  factor  (for  specific  grid) 

p  =  air  density  Qb/ft2) 

2.      Turbulence-Generating  Grids 

Four  grids  are  used  to  create  turbulence  of  varying  intensities  and 
length  scales.  Each  is  mounted  in  a  wooden  frame  and  placed  73  inches  for- 
ward of  the  pivot  axis  of  the  model  support  system  (see  Figure  9).  Three  of  the 
grids  are  constructed  from  wood  and  the  fourth  is  made  of  wire.  Grid  specifi- 
cations are  listed  in  Table   1  and  are  shown  in  Figure  10.  They  are 
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Figure  9.  Plan  view  of  VLSAM  Model  With  Pressure  Probe  and  Grid  In  the 
Test  Section  of  the  Wind  Tunnel  (not  drawn  to  scale)  [Ref.  2] 
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TABLE  1.  GRID  SPECIFICATIONS  [REF.  2] 


Grid 

Mesh  Width  (in.) 

Bar  Diameter  (in.) 

Mesh/Diameter 

Material 

One 

5.00 

1.00 

5 

Wood 

Two 

3.75 

0.75 

5 

Wood 

Three 

2.50 

0.50 

5 

Wood 

Four 

1.00 

0.0625 

16 

Wire 

m 


/%///////     _  Bar  Diameter  (d) 
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Figure  10.  Square-Mesh  Turbulence-Generating  Grid  [Ref.  3] 
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square-mesh  square-bar  biplanar  grids  which  generate  nearly  isotropic  homo- 
geneous turbulence.  [Ref.  29]  Roane  measured  turbulence  intensities  and 
estimated  length  scales,  shown  in  Figures  11  and  12.  [Ref.  1]  The  grid  turbu- 
lence parameters  taken  by  Roane  are  summarized  in  Table  2.  Grid  turbulence 
effects,  with  respect  to  changing  length  scales  at  constant  intensity  or  con- 
stant length  scales  with  changing  intensities,  can  not  be  investigated  with  the 
present  grid  geometries.  [Ref.  3]  Figure  13  shows  photographs  of  the  grids. 
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Figure  11.  Grid  Turbulence  Intensities  (dashed  line  indicates 
model  pivot  axis)  [Ref.  2] 
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Figure  12.  Grid  Turbulence  Length  Scales  [Ref.  2] 


TABLE  2.  GRID  TURBULENCE  PARAMETERS 
(AT  MODEL  PIVOT  AXIS)  [REF.  1] 


Grid 

Intensity 
(percent) 

Length  Scale 
(in.) 

Turbulence/ 
Model  Dia. 

Dynamic  Pressure 
(lb/ft2) 

One 

3.31 

1.84 

1.05 

15.35 

Two 

2.78 

1.56 

0.S9 

14.88 

Three 

1.88 

1.08 

0.62 

16.38 

Four 

0.47 

0.27 

0.15 

15.61 

None 

0.23 

- 

'   - 

15.85 
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Figure  13.  Turbulence-Generating  Grids 
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3.  VLSAM  Model  and  Support  Equipment 

The  model  was  designed  to  represent  a  current  cruciform  tail-con- 
trol missile  with  very  low  aspect  ratio  wings  (long  dorsal  fins).  It  was  con- 
structed from  6061  and  2024  aluminum  alloy  by  NPS  personnel.  [Ref.  1] 

The  hollow  cylinder  body  section  contains  locating  pin  attachment 
points  for  the  balance,  sleeve,  wings  and  tails.  The  machined  sleeve  provides 
a  close  tolerance  fit  between  the  balance  gage  and  the  interior  of  the  model. 
[Ref.  3]  Both  body  roll  angle  and  nose  roll  angle  may  be  varied  in  45°  incre- 
ments. The  wings  with  strakes  and  the  tail  control  fins  are  rigidly  connected 
to  the  model  body  by  countersunk  screws.  Figure  14  depicts  the  dimensions 
and  specifications  of  the  VLSAM  model.  [Ref.  2]  The  model's  surface  is 
polished  and  free  of  protruberances. 

The  model  support  is  rigidly  fixed  in  the  test  section  by  the  reflec- 
tion plane  turntable  at  the  base  and  an  aluminum  reinforced  clear  plexiglass 
section  at  the  top.  The  pivot  point  of  this  rotating  support  coincides  with  the 
approximate  center  of  the  VLSAM  model.  The  plexiglass  has  three  slots  (7-,  8- 
and  10-inches  long)  cut  in  it,  each  5/4  inches  wide.  These  slots  correspond  to 
the  positions  of  model  length  to  diameter  ratios  of  3,  6  and  9;  i.e.,  5.25,  10.5 
and  15.75  inches  from  the  nose.  [Refs.  2  and  3] 

4.  Velmex  8300  3-D  Traverser 

The  Velmex  8300  is  composed  of  a  motor  controller  assembly  and  a 
traversing  assembly,  and  uses  three  microcomputer-controlled  stepping 
motors  (one  for  each  axis  of  movement).  The  motor  controller  assembly  is 
capable  of  interpreting  motor  movement  commands  from  a  host  computer, 


29 


Total  length 

Base  diameter 

Length/diameter  ratio 

Ogive  nose  length 

Ogive/diameter  ratio 

Wing  span/root  chord 

Tail  span/root  chord 

Center  of  pressure 


22.85  in. 

1.75  in. 

13.06 

4.0  in. 

2.29 

3.13  in./ 13.55  in. 

5.50  in./ 1.70  in. 

13.5  inches  aft  of  nose  tip  (approx.) 
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Figure  14.  Specifications  of  VLSAM  Model  [Ref.  2] 


programmable  control  or  terminal.  Software  commands  allow  the  operator  to 
select  motor  variables  such  as  velocity,  acceleration,  increment  distance  and 
units  (motor  steps  or  inches).  [Ref.  32] 
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The  stainless  steel  and  aluminum  traverser  assembly  (Figure  15) 
consists  of  three  separate  motor/jackscrew  assemblies. 


Figure  15.  Velmex  8300  Traversing  Assembly 


The  traverser  was  mounted  to  existing  hardware  on  top  of  the  tun- 
nel so  as  to  minimize  tunnel-induced  vibrations.  A  5-hole  pressure  probe, 
attached  to  the  8300  control  drive,  can  be  accurately  and  effectively  moved 
through  the  test  section. 
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5.  5-Hole  Pressure  Probe 

The  three-dimensional  5-hole  probe,  Figure  16  [Ref.  33],  is  made  of 
corrosion  resistant  non-magnetic  stainless  steel.  It  is  0.125  inch  in  diameter 
and  24  inches  in  total  length  with  22  inches  of  reinforcement  tubing.  The 
probe  has  five  measuring  holes  located  on  its  prism-shaped  tip.  A  centrally 
located  hole  (Pi)  measures  total  pressure,  while  two  lateral  pressure  holes 
(P2,  P3)  are  used  to  determine  yaw  angle  of  flow.  Pitch  angle  is  determined  by 
pressure  holes  (P4,  P5)  located  above  and  below  the  total  pressure  hole.  The 
probe  is  usable  for  speeds  up  to  Mach  0.7. 

The  speed  of  reading  depends  on  the  length  and  diameter  of  the 
pressure  passage  inside  the  probe,  the  size  of  the  pressure  tubes  to  the 
manometer,  and  the  displacement  volume  of  the  manometer.  [Ref.  33]  For 
smaller  diameter  tubes,  the  time  constant  increases  rapidly.  For  this  experi- 
ment, the  tube  diameter  was  1/4-inch  O.D.  and  the  tube  lengths  were  three 
feet,  so  the  time  delay  was  about  0.15-0.26  second.  [Ref.  3] 

6.  Scanivalve  and  HP  Data  Acquisition  System 

One  48-port  scanivalve  was  used  to  measure  each  of  the  5-hole 
probe  pressures.  The  Hewlett-Packard  (HP)  data  acquisition  system  consists 
of  a  combination  of  hardware  and  software  that  enables  the  IBM  PC-AT  com- 
puter to  act  as  a  fully  automated  instrumentation  system.  [Ref.  34]  Individual 
instruments  include  the  Relay  Multiplexer,  Digital  Multimeter  and  Relay 
Actuator. 
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Figure  16.  The  5-Hole  Pressure  Probe  and  Measuring  Tip  [Ref.  33] 


The  scanivalve  mechanism  puts  out  a  7-bit  binary  coded  decimal 
(BCD)  signal  that  corresponds  to  the  port  (1-48)  currently  connected  to  the 
scanivalve  transducer.  This  allows  remote  electronic  monitoring  of  the  port 
assembly  configuration.  [Ref.  3]  The  scanivalve  consists  of  a  transducer, 
motor  drive,  port  assembly  and  solenoid  controller,  which  actually  regulates 
the  scanivalve.  Two  commands  allowed  by  this  solenoid  are  STEP,  which 
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moves  the  scanivalve  to  one  port  location,  and  HOME,  which  sends  the  scani- 
valve  to  port  number  48.  The  Relay  Actuator  is  used  solely  in  controlling  the 
scanivalve  to  STEP  or  HOME.  [Ref.  3] 

The  HP  Data  Acquisition  System  is  shown  in  Figure  17.  The  scanivalve 
signal,  containing  probe  port  pressure  information,  passes  to  the  Relay  Multi- 
plexer which  provides  one  common  output  channel  for  the  Digital  Multimeter 
(DMM).  The  signal  is  conditioned  by  a  low  pass  filter  prior  to  being  measured 
by  the  multimeter.  The  DMM  automatically  converts  input  analog  voltage 
signals  into  a  digital  (binary)  form  which  can  be  read  by  the  computer. 

C.     EXPERIMENTAL  CONDITIONS 

To  facilitate  data  correlation,  the  conditions  for  this  experiment  were 
similar  to  the  conditions  of  the  previous  studies  by  Lung  and  Rabang. 

(1)  Test  section  reference  dynamic  pressures  were  set  at  7.2  cm  H2O  for  the 
no  grid  run  and  at  10.0  cm  H2O  for  grid  #3,  which  yielded  a  subcritical 
Reynolds  number  of  Re=l.lxl05.  These  reference  pressures  are  the 
same  as  those  used  in  the  turbulence  mapping  by  Roane  [Ref.  1]  and 
were  duplicated  in  this  study  in  order  to  ensure  comparable  test 
section  velocities  and  turbulence  grid  length  scales  and  intensities. 

(2)  The  VLSAM  model  nose  geometry  was  held  fixed  at  nose  position  eight, 
which  Rabang  showed  gave  the  maximum  side  force  magnitude. 
[Ref.  2] 

(3)  Afterbody  roll  angles  were  as  follows: 


• 


Body  A 


wings  and  tails  at  roll  angle  <|>r=0o  in  a  "+"  configuration 


•  Body  B:    no  wings  or  tails  at  roll  angle  <j)R=45° 

•  Body  C:    wings  and  tails  at  roll  angle  <J)r=45°  in  a  "x"  configuration. 

Figure  18  shows  these  three  configurations.  Only  Bodies  A  and  C  were 
tested  in  this  study. 
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Figure  17.  HP  Data  Acquisition  System  [Ret  3] 
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Figure  18.  VLSAM  Model  Body  Configurations  [Ref.  2] 


(4)  Model  blockage  factor  corrections  were  calculated  by  Rabang  for  each 
body  configuration.  These  factors  are  a  function  of  the  model  angle  of 

attack.  For  this  experiment,  the  angle  of  attack  was  fixed  at  50°,  where 

the  total  blockage  correction  e  equalled  0.0123.  This  factor  was  imple- 
mented in  data  conversion  programs. 

(5)  The  longitudinal  position  for  data  acquisition  was  at  a  length/body 
diameter  ratio  of  6,  which  was  10.5  inches  from  the  nose  of  the  missile 
model. 

(6)  Wind  tunnel  temperatures  were  not  allowed  to  vary  by  more  than  20°F 
from  the  beginning  to  the  end  of  a  run.  Wind  tunnel  settling  chamber 
temperatures  tended  to  rise  quickly  due  to  air  friction,  particularly 
when  the  grid  was  added.  When  temperatures  were  excessive,  the  tests 
were  stopped  and  the  air  in  the  tunnel  was  circulated  until  it  cooled 
down  sufficiently  before  tests  were  continued. 
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D.     SOFTWARE  AND  PROCEDURES 

In  order  to  correlate  results  with  previous  data,  the  computer  programs 
used  (or  developed)  by  Lung  to  acquire  and  reduce  data  were  also  used  for 
this  study.  Figure  19  [Ref.  3]  is  a  schematic  flowchart  of  the  various  programs 
and  their  resulting  data  files.  The  following  sections  provide  further  elabora- 
tion on  these  programs. 
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Figure  19.  Program/Data  File  Flowchart  [Ref.  3] 


37 


1.      PPROBE  Program  (Data  Acquisition) 

The  BASICA  application  program  which  runs  the  VLSAM  experi- 
ment is  comprised  of  STATEFILE,  PGMSHEL  and  traverser  programs. 
STATEFILE  gives  the  computer  configurations  of  the  data  acquisition 
instruments,  while  PGMSHEL  informs  the  computer  of  all  the  functions 
available  at  each  of  these  instruments.  STATEFILE  and  PGMSHEL  perform 
initialization  chores  and  allow  communication  between  the  HP  instruments 
and  the  IBM  computer.  [Ref.  34] 

The  traverser  program  is  the  actual  application  code  which  allows 
the  operator  to  precisely  control  pressure  probe  movement  by  either  manual 
or  computer-controlled  input.  The  traverser  program  was  written  by 
Kindelspire  [Ref.  35]  in  the  Advanced  Basic  Language.  The  PPROBE  program 
is  shown  in  Appendix  A. 

Manual  control  was  used  to  initialize  the  pressure  probe  position 
prior  to  collecting  data  for  the  actual  run,  which  utilized  computer-controlled 
movements.  While  in  the  manual  mode,  the  program  asked  a  series  of  ques- 
tions which  enabled  the  operator  to  test  motor  movements  and  set  traverser 
motor  velocity  and  acceleration  default  values.  Through  manual  inputs,  the 
probe  was  positioned  such  that  the  Pi  (total  pressure)  hole  was  centered  on 
the  lengthwise  axis  of  the  VLSAM  model  body  and  placed  as  close  as  possible 
to  it.  From  this  point,  the  probe  was  moved  vertically  downward  to  the  posi- 
tion of  the  first  point  in  the  data  collection  plane.  For  this  experiment,  the 
origin  was  located  2.75  inches  below  the  model  axis,  and  the  data  field  dimen- 
sions were  2.75  inches  above  and  below  and  3  inches  outward  from  the  model. 
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After  this  initial  probe  position  was  set,  the  computer-controlled 
motor  movement  option  was  selected.  The  field  dimensions  (x,  y  coordinates), 
the  traverser  motor  step  distance,  and  the  input  file  name  were  entered  into 
the  program.  The  data  plane  was  3.0"  by  5.5"  with  a  0.25-inch  step  distance 
for  this  study.  From  this  input,  PPROBE  then  reiterated  the  total  number  of 
points  to  be  measured  (299  this  case)  and  assigned  filenames  for  each  column 
of  data.  An  example  of  how  data  files  were  named  is  as  follows: 

example:     R0A1A3 


where: 


R  =  run 

0  =  grid  number  (type) 

A  =  VLSAM  model  configuration 

1  =  field  dimension  (3,  5.5) 
A  =  step  distance  (0.25) 

3  =  test  number 


Thus,  from  the  example  above,  the  program  assigned  filenames 
R0A1A301.DAT  through  R0A1A313.DAT,  which  represented  the  13  columns 
of  data  (23  points  per  column). 

The  5-hole  pressure  probe  scale  wheel  was  then  adjusted  until  the 
P2  and  P3  (lateral)  pressures  were  nearly  equal  (nulling),  as  measured  by  a 

portable  digital  manometer/calibrator.  The  measured  yaw  angle  was  read  off 
the  wheel  and  typed  into  the  computer.  Once  this  was  accomplished, 
PPROBE  moved  the  scanivalve  from  port  1  to  port  4  via  the  Relay  Actuator. 
There  was  a  one-second  delay  to  allow  pressure  equalization  before  the 
Digital  Multimeter  sampled  the  output  voltage  from  the  scanivalve  trans- 
ducer via  the  Relay  Multiplexer.  After  ten  samples  were  taken  at  port  4,  the 
Relay  Actuator  stepped  the  scanivalve  to  the  next  port  (5),  where  another  ten 
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samples  were  taken.  This  process  was  repeated  until  all  five  channel  pres- 
sures (ports  4  through  8)  were  measured.  [Ref.  3]  Note  that  scanivalve  port  4 
represents  probe  pressure  Pi,  port  5  is  P2,  port  6  is  P3,  port  7  is  P4,  and  port  8 
isP5. 

The  Relay  Actuator  homed  the  scanivalve  to  port  48  after  all  the 
pressures  were  measured,  and  then  PPROBE  displayed  the  measurements 
and  average  values  for  each  channel  on  the  computer  screen.  The  program 
either  moved  the  probe  upward  one  step  (0.25")  or  remeasured  the  same 
point,  depending  on  whether  the  data  was  within  tolerance  or  not,  as  deter- 
mined by  the  operator.  For  this  study,  the  following  tolerances  were  used:  P2 
and  P3  differed  by  0.09  psf  or  less,  and  Pi  was  a  positive  number  (or  a  very 
small  negative  number  on  the  order  of  -0.5  psD.  Once  a  column  of  data  was 
measured,  PPROBE  would  store  the  average  values  for  each  point  in  a  file  (23 
pts)  and  move  the  traverser  to  the  next  column  position.  The  data  acquisition 
process  was  continued  until  all  299  points  were  completed. 

2.      CALP  Program  (Scanivalve  Transducer  Calibration) 

The  CALP  program  (Appendix  B)  was  the  other  data  acquisition 
program  utilized  in  this  study.  It  was  run  both  before  and  after  the  actual 
test  (PPROBE)  to  account  for  any  change  in  experimental  conditions  which 
might  have  occurred  over  the  8—10  hour  period  it  took  to  run  PPROBE. 

The  transducer  voltage  was  first  adjusted  to  approximately  zero 
millivolts.  Calibration  manometer  (Figure  20)  readings  were  then  entered 
into  the  computer.  The  manometer  provided  a  known  pressure  source  for 
scanivalve  calibration.  From  the  transducer  output  voltage  and  pressure  data 
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Figure  20.  Calibration  Manometer  [Ref.  3] 

provided  by  CALP,  two  calibration  curve  equations  were  calculated  (before/ 
after  PPROBE  runs).  An  averaged  slope  equation  was  then  used  in  data 
reduction  computations  to  minimize  errors.  This  slope  equation  was  used  by 
the  CONVERT  program  to  change  voltage  data  into  dynamic  pressure  data. 
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3.      CONVERT  Program 

This  program  (Appendix  C)  was  used  to  reduce  data.  It  read  the 
PPROBE  data  files  (ex:  R0A1A303.DAT)  containing  average  pressures  (PI- 
PS) and  converted  them  into  x-y  coordinates,  velocity,  yaw  angle,  pitch  angle 
(alpha),  total  local  pressure,  total  pressure  coefficient,  local  static  pressure, 
and  static  pressure  coefficient.  The  x,  y  coordinates  and  yaw  angle  were  sim- 
ply read  from  the  PPROBE  data  files  and  input  into  CONVERT.  The  pitch 
angle  coefficient  was  determined  by  a  ratio  of  the  pressures  measured  by  the 
5-hole  probe.  This  coefficient  was  then  used  in  an  equation  developed  by  Lung 
[Ref.  3]  using  commercial  curve-fitting  software  in  order  to  find  the  corre- 
sponding pitch  angle.  Calibration  curves  provided  by  the  probe  manufacturer 
were  used  for  the  curve-fit.  The  program  also  obtained  velocity  pressure  coef- 
ficients for  particular  pitch  angles,  from  which  the  local  velocity  was  calcu- 
lated after  addition  of  a  wind  tunnel  calibration  factor,  K,  and  other  unit  con- 
version factors.  Additionally,  the  total  and  static  pressure  coefficients  were 
obtained  for  different  pitch  angle  regions. 

Room  ambient  pressure  was  used  as  the  reference  pressure  for  these 
coefficients,  which  are  functions  of  total  and  static  pressures  and  are  non- 
dimensionalized  by  the  tunnel  dynamic  pressure: 

Cps  =  (PsL-Ps)/Q  (6) 

Cpr  =  (PtL-PtVQ  (7) 

where 

Cps  =  Static  pressure  coefficient 

Cpr  =  Total  pressure  coefficient 

Q  =  Freestream  dynamic  pressure 

Ps  =  Freestream  static  pressure 

Pt  =  Freestream  total  pressure 

Psl  =  Local  static  pressure 

PtL  =  Local  total  pressure 

The  actual  dynamic  pressure  is  (nearly)  the  same  for  the  two  cases. 
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For  this  experiment,  the  reference  dynamic  pressure  values  were  7.2  cm 
H2O  for  grid  0  and  10.0  cm  H2O  for  grid  3.  Temperature  input  for  the 
CONVERT  program  was  an  average  of  the  initial  and  final  wind  tunnel 
temperatures  for  the  entire  run  time.  Similarly,  barometric  pressure  values 
were  recorded  before  and  after  each  run.  CONVERT  also  added  yaw  and 
pitch  angle  (alpha)  corrections  to  the  data.  These  factors  are  +5.0°  for  yaw 
and  -17.942°  for  pitch.  They  were  determined  from  a  preliminary  run 
conducted  with  no  grid  and  no  missile  model  in  the  wind  tunnel.  Further 
explanation  of  this  run  is  discussed  in  the  Preliminary  Tests  section.  The 
output  of  the  CONVERT  program  was  stored  in  a  file  named  RESULT.DAT, 
which,  in  turn,  was  used  as  input  to  the  TECPLOT  system. 

4.      TECPLOT 

The  commercial  TECPLOT  software  system  was  used  to  generate 
crossplane  velocity  vector  plots  and  pressure  contour  plots.  These  plots  could 
be  tailored  in  many  different  ways  by  choosing  scale  factors,  arrowhead 
wedge  angles,  contour  levels  and  spacing,  and  many  other  parameters.  A 
Hewlett-Packard  7470A  x-y  pen  plotter  was  utilized  in  conjunction  with 
TECPLOT  to  provide  both  the  vector  and  contour  plots. 

E.     PRELIMINARY  TESTS 

1.      Dynamic  Pressure  Calibration 

All  of  the  turbulence  grids  were  previously  calibrated  in  the  wind 
tunnel.  Readings  from  the  tunnel  calibration  manometer  and  from  a  pi  tot- 
static  tube  inserted  in  the  center  of  the  test  section  were  recorded  over  a 
speed  range  and  wind  tunnel  calibration  factors  were  obtained.  These  factors 
were  used  to  adjust  the  tunnel  flow  velocity  to  the  expected  experimental 
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condition  for  the  different  grids.  [Ref.  3]  The  calibration  factors  K  are  0.8891, 
1.5084,  1.6487,  1.6545  and  1.1167  for  no  grid,  and  grids  1  through  4  respec- 
tively. These  values  were  used  in  the  CONVERT  program  to  calculate  the 
pressure  and  velocity  in  the  test  section. 

2.      Yaw  and  Pitch  Angle  Corrections 

A  test  was  conducted  to  find  correction  factors  for  yaw  and  pitch 
angle.  From  previous  arrow  plot  data  by  Lung  [Ref.  3],  inconsistent  crossflow 
velocity  magnitudes  and  directions  were  noted  toward  the  outer  boundaries  of 
the  body-only  missile  configuration  run.  These  outer  boundaries  represented 
the  wind  tunnel  freestream  region,  where  crossflow  velocity  is  expected  to 
reach  zero.  Therefore,  to  duplicate  just  the  freestream  region,  this  prelimi- 
nary test  consisted  of  placing  the  pressure  probe  in  the  tunnel  with  no 
VLSAM  model  and  no  grid.  Thus,  the  expected  pitch  and  yaw  angles  should 
both  have  been  zero.  This  was  not  the  case  however. 

Results  of  the  preliminary  run  (R001A2)  are  listed  in  Appendix  D 
(Result  00.DAT  File).  The  yaw  angles  measured  ranged  from  -4.00°  to  -7.00°, 
with  an  average  of  -5.00°.  The  pitch  angles  (alpha)  ranged  from  approxi- 
mately +  17.4°  to  approximately  +18.4°,  with  the  average  +17.942°.  Thus,  to 
correct  for  these  errors,  +5.00°  was  added  to  the  yaw  angle  and  -17.942°  was 
subtracted  from  the  alpha  values  in  the  output  file  of  the  CONVERT  program 
(RESULT.DAT). 

Though  the  exact  cause  for  the  yaw  and  pitch  angle  errors  was  not 
known,  one  possibility  might  have  been  a  slight  bend  which  was  noted  in  the 
5-hole  pressure  probe.  Other  causes  may  have  been  improper  calibration  of 
the  probe  or  a  misalignment  of  the  traverser  assembly.  The  corrected  errors 
only  effect  the  velocity  vector  plots,  and  should  have  no  effect  on  the  pressure 
contour  plots. 
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III.  RESULTS 

The  following  sections  discuss  the  velocity  vector  plots  and  the  total  and 
static  pressure  coefficient  contour  plots  for  the  VLSAM  model  configurations 
A  (plus)  and  C  (cross),  both  with  and  without  turbulence.  All  plots  depict  the 
3"  by  5.5"  data  acquisition  field  and  the  position  of  the  missile  model  (nose 
aspect)  relative  to  the  field.  Vortex  sizes  and  vortex  distances  from  the  model 
surface  are  referenced  to  the  model  base  diameter  d  (1.75"). 

A.     CONFIGURATION  OA  ("PLUS"  WITHOUT  TURBULENCE) 

For  the  plus  configuration,  the  swirling  patterns  of  the  velocity  vector 
plot  (Figure  21)  clearly  denote  the  two  asymmetric  vortices,  which  form 
circles  that  rotate  in  opposite  directions.  Although  the  bottom  vortex  center  is 
evident,  this  plot  fails  to  show  the  center  of  the  top  vortex.  The  vortex 
strength  is  a  maximum  on  the  outer  edges  of  the  vortex  cores,  denoted  by  the 
large  vector  arrows,  where  the  velocities  flow  back  toward  the  missile  body. 
Towards  the  outer  boundaries  of  the  data  acquisition  field,  where  the  vortex 
strengths  are  minimal,  the  vectors  plot  as  points. 

The  total  pressure  coefficient  (Cpr)  contour  plot  (Figure  22)  shows  that 
the  extent  of  the  top  vortex  is  approximately  0.72d  at  a  distance  of  0.33d  from 
the  missile  surface.  The  bottom  vortex  is  0.83d  at  0.67d  from  the  body.  There 
are  more  changes  in  the  pressure  gradient  within  a  smaller  area  for  the  top 
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Figure  21.  Velocity  Vector  Plot  -  Configuration  0A 
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Figure  22.  Total  Pressure  Coefficient  -  Configuration  0A 
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vortex,  which  varies  in  Cpx  from  -0.14  to  -1.35  within  an  inch,  indicating 
greater  vortex  strength.  The  bottom  vortex  is  slightly  more  diffused.  The  Cpt 
plot  shows  that  the  distance  between  the  two  vortex  centers  is  approximately 
0.77d. 

The  static  pressure  coefficient  (Cps)  contour  plot  (Figure  23)  shows  that 
the  top  vortex  extends  approximately  l.Od  at  a  distance  of  0.49d  from  the 
body.  The  bottom  vortex  extends  about  the  same  (l.ld)  at  a  distance  of  0.55d. 
Again,  the  top  vortex  appears  to  be  stronger  than  the  bottom  vortex.  Cps  for 
the  top  vortex  varies  from  -0.55  to  -2.70  within  1.2  inches  while  the  bottom 
vortex  varies  from  -0.55  to  about  -2.18  in  1.6  inches.  The  distance  between  the 
vortices  is  approximately  0.70d  on  the  Cps  plot. 

B.     CONFIGURATION  3A  ("PLUS"  WITH  TURBULENCE) 

With  added  turbulence  (grid  3),  the  velocity  vector  plot  (Figure  24)  still 
indicates  vortex  asymmetry,  but  it  also  indicates  that  the  vortices  have  less 
strength  (smaller  vector  arrows)  than  for  the  no  grid  condition.  Again  the 
vortex  strength  is  maximized  on  the  outer  edges  of  the  vortex  cores. 

From  the  Cpr  contour  plot  (Figure  25),  the  top  vortex  extends  to  approx- 
imately 0.72d,  centered  at  about  0.44d  from  the  model.  The  bottom  vortex  is 
slightly  larger  (0.83d)  and  is  located  0.55d  from  the  model  body.  The  pressure 
gradient  of  the  top  vortex  is  steeper  than  for  the  bottom  vortex,  indicating 
greater  strength.  Cpx  varies  from  -0.14  to  -1.5  within  an  inch  for  the  top 
vortex  and  from  -0.14  to  -1.35  within  1.4  inches  for  the  bottom  vortex.  The 
distance  between  the  two  vortices  as  measured  on  the  Cpx  plot  is  approxi- 
mately 0.72d. 
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Figure  23.  Static  Pressure  Coefficient  -  Configuration  OA 
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Figure  24.  Velocity  Vector  Plot  -  Configuration  3 A 
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Figure  25.  Total  Pressure  Coefficient-  Configuration  3A 
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Examination  of  the  Cps  contour  plot  (Figure  26)  reveals  that  both 
vortices  are  roughly  0.83d  wide.  The  top  vortex  is  closer  to  the  missile  body 
(0.39d)  than  the  bottom  vortex  (0.55d).  More  contour  levels  in  a  smaller  area 
and,  consequently,  a  higher  gradient  and  stronger  vortex  exists  for  the  top 
vortex.  Cps  for  the  top  vortex  varies  from  -0.27  to  -1.60  within  0.8  inch  while 
the  bottom  vortex  Cps  ranges  from  -0.27  to  -1.36  within  1.0  inch.  On  the  Cps 
plot,  the  two  vortices  are  about  0.66d  apart. 

C.     CONFIGURATION  0C  ("CROSS"  WITHOUT  TURBULENCE) 

For  the  model  C  configuration  (cross),  the  velocity  vector  plot  (Figure  27) 
still  clearly  displays  the  two  asymmetric  vortices.  As  with  configuration  A, 
the  strengths  of  the  vortices  are  maximized  at  the  edges  of  the  circular  swirls, 
where  the  flow  is  back  toward  the  missile  body. 

Figure  28,  the  total  pressure  coefficient  plot,  displays  a  0.77d  wide  top 
vortex  located  at  about  0.5d  from  the  model  surface.  The  bottom  vortex 
extends  to  approximately  0.88d  and  is  0.66d  from  the  model.  The  bottom 
vortex  is  more  tightly  wrapped  (i.e.,  more  contour  levels  per  area)  closer  to  its 
core  than  is  the  top  vortex.  The  Cpx  range  for  the  bottom  vortex  (-0.14  to 
-1.23)  is  slightly  greater  than  the  range  for  the  top  vortex  (-0.14  to  -1.09).  The 
vortex  centers  are  roughly  0.66d  apart. 

From  the  static  pressure  coefficient  contour  (Figure  29),  both  vortices 
extend  to  0.95d  with  the  top  vortex  slightly  closer  to  the  model  body  (0.56d) 
than  the  bottom  vortex  (0.66d).  The  relative  strengths  of  the  vortices  is  diffi- 
cult to  interpret  from  the  Cps  plot,  but  the  top  vortex  appears  to  be  a  little 
less  diffused  than  the  bottom  vortex  and  therefore  stronger.  As  with  the  total 
pressure  coefficient  plot,  the  vortices  are  0.66d  apart. 
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Figure  26.  Static  Pressure  Coefficient  -  Configuration  3A 
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Figure  28.  Total  Pressure  Coefficient  -  Configuration  OC 
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Figure  29.  Static  Pressure  Coefficient  -  Configuration  OC 
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D.  CONFIGURATION  3C  ("CROSS"  WITH  TURBULENCE) 

The  velocity  vector  plot  (Figure  30)  shows  that  the  addition  of  turbulence 
(grid  3)  diffuses  the  asymmetric  vortices  for  model  configuration  C.  The  vector 
arrows  plot  smaller  for  this  turbulent  case  than  for  the  no  grid  case. 

The  Cpr  contour  plot  (Figure  31)  reveals  a  0.92d  bottom  vortex  at  a  dis- 
tance of  0.66d  from  the  model  body  and  a  0.83d  wide  top  vortex  at  0.55d  from 
the  model.  The  bottom  vortex  appears  to  be  slightly  more  tightly  wrapped 
towards  the  vortex  center  than  the  top  vortex.  Cpr  levels  for  both  vortices 
range  from  about  0.0  to  -1.23.  The  distance  between  the  vortices  is  approxi- 
mately 0.72d. 

Figure  32,  the  static  pressure  coefficient  plot,  like  the  Cpt  plot,  also 
shows  the  distance  between  the  vortices  to  be  0.72d.  The  top  vortex  (extent 
0.77d)  is  roughly  0.55d  from  the  model's  surface  while  the  bottom  vortex 
(extend  0.88d)  is  0.66d.  The  vortices  appear  to  have  the  same  strength  on  the 
Cps  plot,  in  which  the  coefficients  range  from  0.0  to  approximately  -1.1. 

E.  COMPARISONS 

1.      Between  Body  Configurations  (A  and  C) 

The  following  observations  were  noted  when  making  comparisons 

between  the  two  missile  configurations  for  both  the  turbulent  (grid  3)  and 

non-turbulent  (no  grid)  runs. 

(1)  From  the  total  pressure  coefficient  (Cpr)  contour  plots,  the  vortices 
are  slightly  larger  for  configuration  C.  This  is  more  pronounced  for 
the  run  conducted  with  added  turbulence  (grid  3). 
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Figure  30.  Velocity  Vector  Plot  -  Configuration  3C 
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Figure  31.  Total  Pressure  Coefficient  -  Configuration  3C 


59 


2.0 


1.0 


^      0.0 

< 


-1.0 


•2.0 


j '    '     ' i '    i     t    '    '     '     '    '    i    < i i    i    i    i    i    i    t 


Level 

CPS 

6 

-0.27 

5 

-0.82 

4 

-1.3S 

3 

-1.91 

2 

-2.45 

1 

-3.00 

MISSILE 


-2.0 


1.0  0.0  1.0 

X-AXIS  (IN) 


2.0 


Figure  32.  Static  Pressure  Coefficient  -  Configuration  3C 
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(2)  The  vortices  of  the  C  configuration  appear  to  be  centered  a  little  fur- 
ther away  from  the  model  surface  than  do  the  A  configuration 
vortices.  This  shift  is  on  the  order  of  0.2  inches,  and  is  apparent  in 
both  the  no-grid  and  grid  3  runs. 

(3)  The  relative  distance  between  vortex  centers  does  not  seem  to  be 
affected  when  the  configurations  are  changed. 

(4)  The  bottom  vortex  is  more  tightly  wrapped  (stronger)  for  configura- 
tion C  based  on  the  Cpt  plots.  This  is  more  noticeable  for  the  run 
conducted  without  turbulence. 


2.      Between  Turbulence  Levels  (0  and  3) 

The  following  observations  were  noted  when  making  comparisons 
between  the  two  turbulence  levels  tested  for  each  body  configuration. 
Turbulence  level  0  was  the  no-grid  condition  and  level  3  was  the  run  con- 
ducted with  a  grid  inserted. 

(1)  The  velocity  vector  plots  indicate  smaller  arrows,  therefore  weaker 
vortices  (less  crossflow),  when  turbulence  is  added.  This  trend  holds 
for  both  body  configurations. 

(2)  The  two  vortices  remain  in  approximately  the  same  positions  as  indi- 
cated on  both  the  vector  plots  and  the  static  pressure  coefficient  (Cps) 
plots). 

(3)  From  the  total  press  coefficient  (Cpt)  plots,  the  center  of  the  bottom 
vortex  for  configuration  A  shifts  slightly  closer  to  the  model  body 
when  the  turbulence  increases.  (-0.2  to  0.0  on  the  x-axis). 

(4)  Also  on  the  Cpr  plots,  with  added  turbulence,  the  top  vortex  core  for 
configuration  C  seems  to  shift  away  from  the  VLSAM  body. 

(5)  The  Cpr  plots  for  both  configurations  are  slightly  more  diffused  with 
added  turbulence.  Specifically,  the  vortex  centers  are  less  tightly 
wrapped,  thus  indicating  weaker  vortices  for  the  more  turbulent 
condition. 

(6)  There  are  no  noticeable  differences  in  the  static  pressure  coefficient 
contours  between  the  turbulent  and  non- turbulent  runs. 
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3.      With  Body-Only  Configuration  (B) 

The  following  observations  were  made  in  order  to  correlate  the 
results  of  this  study  with  the  previous  experiment  conducted  by  Lung  for  a 
body-only  VLSAM  model  configuration.  [Ref.  3]  The  purpose  is  to  note  how 
the  addition  of  wings  and  strakes  might  affect  the  vortex  flow  pattern  around 
the  missile.  The  velocity  vector,  total  pressure  and  static  pressure  coefficient 
plots  for  the  body-only  run  (OB)  are  displayed  in  Figures  33,  34,  and  35 
respectively.  These  comparisons  were  made  for  the  nominal  ambient  (no-grid) 
flowfield  condition. 

(1)  From  the  velocity  vector  plots,  the  vortex  pattern  for  configuration  B 
is  similar  to  the  one  for  body  C  in  that  the  vortex  asymmetry  is  more 
pronounced  for  these  cases.  Body  A  vortices  are  closer  to  the  model 
surface  and  asymmetry  is  less  pronounced. 

(2)  From  the  total  pressure  coefficient  contours,  the  vortices  appear  to  be 
located  much  closer  to  the  missile  body  for  body  A  than  for  either 
body  C  or  B.  Apparently  the  strake/wing-generated  vortices  for  this 
configuration  act  to  hold  the  nose-generated  vortices  near  the  body, 
resulting  in  the  higher  induced  side  forces  for  this  configuration 
observed  by  Rabang. 

(3)  Also  from  the  Cpr  plots,  the  strengths  of  the  bottom  vortices  appear 
to  be  roughly  the  same  for  bodies  A  and  B,  while  the  plot  for  body  C 
seems  to  be  more  tightly  wrapped  (stronger)  near  its  center. 

(4)  From  the  static  pressure  contours,  the  body-only  vortices  are  more 
diffused  and  larger  than  either  configuration  with  wings  (A  and  C). 

(5)  The  relative  distances  between  the  top  and  bottom  vortices  remains 
the  same  on  all  three  plots  for  all  body  configurations. 

In  general,  the  vortex  pattern  around  the  nose  of  the  missile  model  with 

wings  (in  either  configuration)  resembles  the  vortex  pattern  for  a  body-only 

configuration.  Though  there  are  subtle  differences  as  previously  noted,  the 
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relative  strengths,  sizes  and  positions  of  the  asymmetric  vortices  were  com- 
parable in  all  cases. 
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Figure  33.  Velocity  Vector  Plot  -  Configuration  0B 
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Figure  34.  Total  Pressure  Coefficient  -  Configuration  OB 
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Figure  35.  Static  Pressure  Coefficient  -  Configuration  OB 
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IV.  CONCLUSIONS  AND  RECOMMENDATIONS 

The  flowfield  about  a  Vertically-Launched  Surface-to-Air  Missile 
(VLSAM)  model  at  high  angle  of  attack  was  investigated  in  the  wind  tunnel 
of  the  Naval  Postgraduate  School.  Missile  "plus"  and  "cross"  body  configura- 
tions (A  and  C  respectively)  were  both  tested.  The  angle  of  attack  was  set  at 
50°  with  a  Reynolds  number  of  1.1x10s  for  all  data  runs.Two  flowfield  condi- 
tions were  treated:  the  nominal  ambient  wind  tunnel  condition  (no  grid)  and 
a  condition  with  a  grid-generated  turbulence  of  length  scale  1.08  inches  and 
1.88%  turbulence  intensity  (see  Table  2).  The  following  conclusions  were 
reached: 

(1)  An  increase  in  turbulence  intensity  tended  to  reduce  the  strength  of 
the  asymmetric  nose-generated  vortices. 

(2)  The  two  asymmetric  vortices  remained  in  approximately  the  same 
position  for  an  increase  in  turbulence. 

(3)  The  top  vortex  was  closer  to  the  model  surface  and  appeared  to  be 
stronger  for  both  body  configurations.  This  condition  was  more  pro- 
nounced for  configuration  A. 

(4)  The  wing/strake  arrangement  of  configuration  C  caused  the  vortices 
to  be  centered  further  away  from  the  model  surface  than  those  of  con- 
figuration A,  correlating  with  the  differences  in  induced  side  forces 
for  those  configurations  observed  by  Rabang. 

(5)  Configuration  C  vortices  are  more  diffused  and  larger.  This  was  more 
apparent  when  turbulence  was  added. 

(6)  Crossflow  velocity  vector  plots  agreed  with  the  behavior  denoted  by 
the  total  and  static  pressure  coefficient  contours  for  both  body  config- 
urations and  turbulence  levels. 
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(7)  Though  subtle  differences  exist,  the  addition  of  wings  and  tails  did 
not  greatly  alter  the  vortex  pattern  around  the  nose  of  the  missile 
model,  when  compared  to  the  body-only  configuration  tested  by  Lung. 

Recommendations  for  a  continued  study  of  the  behavior  of  asymmetric 
vortices  under  varying  flowfield  conditions  are  suggested  as  follows: 

(1)  Examine  the  vortices  at  positions  further  back  along  the  model  body, 
where  effects  from  the  wings  might  better  be  seen.  (Such  as  at  a 
length/diameter  ratio  of  9.) 

(2)  Investigate  asymmetric  vortex  behavior  for  just  a  body-only  configu- 
ration, at  the  position  described  above,  to  provide  comparisons. 

(3)  Continue  to  study  vortex  behavior  at  various  angles  of  attack  and 
turbulence  with  varying  intensities  and  length  scales  in  order  to  pro- 
vide a  large  data  base,  which  can  be  used  to  calculate  vorticity 
contours. 
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APPENDIX  A-  PPROBE  PROGRAM 


1  DEF  SEG: CLEAR  , &HFE00:GOTO  4  'Begin  PCIB  Program  Shell 

2  GOTO  1000  '  User  program 

3  GOTO  900  '  Error  handling 

4  I=&HFE00  '  Copyright  Hewlett-Packard  1984,1985 

5  PCIB.DIR$=ENVIRON$("PCIB" ) 

6  I$=PCIB.DIR$+"\PCIBILC.BLD" 

7  BLOAD  1  $ , I 

8  CALL  KPCIB.DIRS,  I%,J7„)  : PCIB . SEG= 1% 

9  IF  J%=0  THEN  GOTO  13 

10  PRINT  "Unable  to  load."; 

11  PRINT  "    (Error  #";J%;")" 

12  END 

13  ' 

l't    DEF  SEG=PCIB.SEG:O.S=5:C.S=in:I .V-15 

15  I .C=20:L.P=25:LD.FILE=30 

16  GET.MEM=35:L.S=40:PANELS=45:DEF.EPR=50 

17  PCIB.ERR$=STRING$< 64 , 32 )  :  PCIB . NAMES=STRING$ ( 16 , 32  ) 

18  CALL  DEF. ERR(PCIB. ERR, PCIB. ERRS, PCIB. NAME* , PCIB. GLBERR)  :  PCI U . BASERR= 2 55 

19  ON  ERROR  GOTO  3 

20  J=-l 

21  I$=PCJB.DIR$+"\PCIB.SYN" 

22  CALL  O.S( If ) 

23  IF    PC1B.ERRO0    THEN    ERROR    PC1B.BASERR 

24  1  =  0 

25  CALL  I . V( I .READ. REGISTER, READ. SELF1D. DEFINE. INITIALI ZE . SYSTEM ) 
2  6  IF  PCIB.ERRO0  THEN  ERROR  PCIB.BASERR 

27  CALL  I .V( I .ENABLE. SYSTEM. DISABLE. SYSTEM, INITIALIZE , POWER . ON ) 

28  IF  PCIB.ERRO0  THEN  ERROR  PCIB.BASERR 

2  9  CALL  I.V( I .MEASURE. OUTPUT. START, HALT)  i 

30  IF  PCIB.ERRO0  THEN  ERROR  PCIB.BASERR 

31  CALL  I . V( I .ENABLE. INT . TRIGGER , DISABLE . INT . TRIGGER. ENABLE . OUTPUT , DISABLE . OUTPU 
T) 

3  2  IF  PCIB.ERRO0  THEN  ERROR  PCIB.BASERR 

3  3  CALL  I . V( 1 .CHECK. DONE. GET. STATUS. SET. FUNCTION. SET. RANGE) 
34  IF  PCIB.ERRO0  THEN  ERROR  PCIB.BASERR 

3  5  CALL  I .V( I .SET. MODE, WRITE. CAL, READ. CAL, STORE. CAD 

36  IF  PC1B.ERRO0  THEN  ERROR  PCIB.BASERR 

37  CALL  I . V( I .DELAY. SAVE. SYSTEM, J, J ) 

38  IF  PCIB.ERRO0  THEN  ERROR  PCIB.BASERR 

39  1  =  1 

40  CALL  I  .V(  I,  SET.  GATF.TIME,  SET.  SAMPLES,  SET.  SLOPE,  SET.  SOURCE) 

41  IF  PCIB.ERRO0  THEN  ERROR  PCIB.BASERR 

4  2  CALL  I .C( I. FREQUENCY, AUTO. FREQ, PERIOD, AUTO. PER) 

43  IF  PCIB.ERRO0  THEN  ERROR  PCIB.BASERR 

44  CALL  I .C( I . INTERVAL. RATIO, TOTALIZE, R100MILLI ) 

45  IF  FCIB.ERRO0  THEN  ERROR  PCIB.BASERR 

46  CALL  I .C(I ,R1,R10,R100,R1KILO) 

47  IF  PCIB.ERRO0  THEN  ERROR  PCIB.BASERR 

48  CALL  I .C( I,R10MEGA,R100MEGA,CHAN.A,CHAN.B) 

49  IF  PCIB.ERRO0  THEN  ERROR  PCIB.BASERR 

50  CALL  I .C( I .POSITIVE, NEGATIVE.COMN, SEPARATE) 

51  IF  PCIB.ERRO0  THEN  ERROR  PCIB.BASERR 

52  1  =  2 

53  1  =  3 
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5k    CALL    I .V( I, ZERO. OHMS. SET. SPEED. J. J) 
5  5    IF    PCTB.ERROO    THEN    ERROR    PCIB.BASERR 

56  CALL    I .C( I ,DCVOLTS,ACVOLTS,OHMS,R200MILLI ) 

57  IF    PCIB.ERROO    THEN    ERROR    PCIB.BASERR 

58  CALL    I .C( I,R2,R20,R200.R2KILO) 

59  IF    PCIB.ERROO    THEN    ERROR    PCIB.BASERR 

60  CALL    I.C(I,R20KILO,R200KILO,R2MEGA,R20MEGA) 

61  IF    PCIB.ERROO    THEN    ERROR    PCIB.BASERR 

62  CALL    I.C(I,AUTOM.R2.5,R12.5,J) 

63  IF  PCIB.ERROO  THEN  ERROR  PCIB.BASERR 
6k    l=k 

6 5  CALL  I . V ( I , SET . COMPLEMENT , SET . DRIVER , OUTPUT . NO . WA IT , ENABLE . HANDSHAKE ) 

66  IF  PCIB.ERROO  THEN  ERROR  PCIB.BASERR 

67  CALL  I .V( I .DISABLE. HANDSHAKE, SET. THRESHOLD, SET. START. BIT. SET. NUM. BITS) 

68  IF    PCIB.ERROO    THEN    ERROR    PCIB.BASERR 

69  CALL  I .V( I , SET. LOGIC. SENSE, J, J, J) 

70  IF    PCIB.ERROO    THEN    ERROR    PCIB.BASERR 

71  CALL  I .C( I .POSITIVE, NEGATIVE, TWOS, UNSIGNED) 

72  IF    PCIB.ERROO    THEN    ERROR    PCIB.BASERR 
7  3    CALL    I .C(I .OC.TTL.RO.Rl) 

7k    IF    PCIB.ERROO    THEN    ERROR    PCIB.BASERR 
75    CALL    I .C(I ,R2,R3,Rd,R5) 

7  6    IF    PCIB.ERROO    THEN    ERROR    PCIB.BASERR 

77  CALL    I .C( I ,R6,R7.R8,R9) 

78  IF    PCIB.ERROO    THEN    ERROR    PCIB.BASERR 

79  CALL    I .C( I .R10.R11.R12.R13 ) 

80  IF  PCIB.ERROO  THEN  ERROR  .  PCI  B  .  BASERR 

81  CALL    I  .C(  I  ,R1<*,R15,R16,J  )  ' 

8  2    IF    PCIB.ERROO    THEN    ERROR    PCIB.BASERR 
83    1=6 

Bk    CALL  I .V( I , SET. FREQUENCY, SET. AMPLITUDE, SET. OFFSET, SET. SYMMETRY) 

85    IF    PCIB.ERROO    THEN    ERROR    PCIB.BASERR  i 

8  6  CALL  I .V( I .SET. BURST. COUNT, J, J, J) 

87  IF    PCIB.ERROO    THEN    ERROR    PCIB.BASERR 

88  CALL  I .C( I .SINE, SQUARE , TRIANGLE , CONTINUOUS) 

89  IF  PCIB.ERROO  THEN  ERROR  PCIB.BASERR 

90  CALL  I .C( I .GATED, BURST, J, J) 

91  IF    PCIB.ERROO    THEN    ERROR    PCIB.BASERR 

92  1  =  7 

93  CALL  I . V( I . AUTOSCALE, CALIBRATE, SET. SENSITIVITY, SET. VERT. OFFSET) 
9k    IF  PCIB.ERROO  THEN  ERROR  PCIB.BASERR 

95  CALL  I .V( I , SET. COUPLING, SET. POLARITY, SET. SWEEPSPEED, SET. DELAY) 

96  IF    PCIB.ERROO    THEN   ERROR   PCIB.BASERR 

97  CALL  I . V( I , SET. TRIG. SOURCE, SET. TRIG. SLOPE, SET. TRIG. LEVEL, SET. TRIG. MODE) 

98  IF    PCIB.ERROO    THEN    ERROR    PCIB.BASERR 

9  9  CALL  I .V( I ,GET.SINGLE.WF,GET.TWO.WF,GET.VERT. INFO, GET  TIMEBASE. INFO) 

100  IF    PCIB.ERROO    THEN    ERROR    PCIB.BASERR 

101  CALL  I .V( I .GET. TRIG. INFO , CALC .WFVOLT , CALC .WFTIME , CALC . WF . STATS ) 

102  IF  PCIB.ERROO  THEN  ERROR  PCIB.BASERR 

103  CALL  I .V( I . CALC. RISETI ME. CALC. FALLTI ME, CALC. PERIOD, CALC. FREQUENCY) 
10<»  IF  PCIB.ERROO  THEN  ERROR  PCIB.BASERR 

105  CALL  I .V( I , CALC. PLUSWIDTH, CALC. Ml NUSWIDTH. CALC. OVERSHOOT, CALC. PRESHOOT) 

106  IF    PCIB.ERROO    THEN    ERROR    PCIB.BASERR 

107  CALL  I .V( I , CALC. PK. TO. PK, SET. TIMEOUT, SCOPE. START, MEASURE. SINGLE. WF) 

108  IF    PCIB.ERROO    THEN    ERROR    PCIB.BASERR 

10  9  CALL  I .V( I , MEASURE. TWO. WF, J, J. J ) 
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110  IF  PCIB.ERRO0  THEN  ERROR  PCIB.BASERR 

111  CALL  I .C( I ,R10NANO,R100NANO,R1MICKO,R10MICRO) 

112  IF  PCIB.ERRO0  THEN  ERROR  PCIB.BASERR 

113  CALL  I .C( I .R100MICRO.R1MILLI , R10MILLI , R100MILLI ) 

114  IF  PCIB.ERRO0  THEN  ERROR  PCIB.BASERR 

115  CALL  I .C(I,R1.R3  0,R20NANO,R200NANO) 

116  IF  PCIB.ERRO0  THEN  ERROR  PCIB.BASERR 

117  CALL  I .C( I ,R2MICRO,R20MICRO,R200MICRO,R2MILLI ) 

118  IF  PCIB.ERRO0  THEN  ERROR  PCIB.BASERR 

119  CALL  I .C( I .R20MILLI .R200MILLI .R2.R20) 

120  IF  PCIB.ERRO0  THEN  ERROR  PCIB.BASERR 

121  CALL  I .C(I ,R50NANO,R500NANO,R5MICRO.R50MICRO) 

122  IF  PCIB.ERRO0  THEN  ERROR  PCIB.BASERR 

123  CALL  I .C( I .R500M1CRO.R5M1LLI .R50MILLI .R500MILLI ) 

124  IF  PCIB.ERRO0  THEN  ERROR  PCIB.BASERR 

125  CALL  I .C(I ,R5,R50,CHAN.A,CHAN.B) 

126  IF  PCIB.ERRO0  THEN  ERROR  PCIB.BASERR 

127  CALL  I .C( I, EXTERNAL, POSITIVE, NEGATIVE, AC) 

128  IF  PCIB.ERRO0  THEN  ERROR  PCIB.BASERR 

129  CALL  I . C( I ,DC, TRIGGERED, AUTO. TRIG, AUTO. LEVEL) 

130  IF  PCIB.ERRO0  THEN  ERROR  PCIB.BASERR 

131  CALL  I .C( I.X1.X10. STANDARD, AVERAGE) 

132  IF  PCIB.ERRO0  THEN  ERROR  PCIB.BASERR 

133  1=8 

3  34  CALL  I.  V(l,  OPEN.  CHANNEL,  CLOSE.  CHANNEL,  J,  J) 

135  IF    FCIB.ERRO0    THEN    ERROR    PCIB.BASERR 

136  CALL    C.S 

137  IF  PCIB.ERRO0  THEN  ERROR  PCIB.BASERR 

138  I$=PCIB.DIR$+"\PCIB.PLD" 

139  CALL  L.P(IS ) 

140  IF  PCIB.ERRO0  THEN  ERROR  PCIB.BASERR 

141  I$="DMM.01": 1-3 :J=0:K=0:L=1 

14  2  CALL  DEFINE(DMM.01, 1$, I, J.K.L) 

143  IF  PCIB.ERRO0  THEN  ERROR  PCIB.BASERR 

li*i*    I$="Func.Gen.01":I=6:J=0:K=l:L=l  ' 

li*5    CALL  DEFINE  (Func.  Gen.  01,  1$,  I  ,J,K,L) 

11*6    IF    PCIB.ERRO0    THEN    ERROR    PCIB.BASERR 

14  7    I$="Scope.01": 1  =  7  :  J=0 :K=2 :L=1 

148  CALL  DEFINE (Scope. 01, 1$, I ,J,K,L) 

149  IF  PCIB.ERRO0  THEN  ERROR  PCIB.BASERR 

150  1$= "Counter. 01": 1=1 : J=0 :K=3 :L=1 

151  CALL   DEFINE(Counter.01 , 1$, I ,J,K,L) 

152  IF    PCIB.ERRO0    THEN    ERROR    PCIB.BASERR 

15  3    I$="Dig. In.01":I=4:J=0:K=4:L=l 

154  CALL   DEFINE(Dig. In. 01, 1$, I ,J,K.L) 

155  IF    PCIB.ERRO0    THEN    ERROR    PCIB.BASERR 

156  I$="Dig.Out.01": 1=4 : J=1:K=4 :L=1 

157  CALL    DEFINE(Dig.Out.01 , If , I , J.K.L) 

158  IF    PCIB.ERRO0    THEN    ERROR    PCIB.BASERR 

159  I $= "Relay. Act. 01": I=8:J=0:K=5:L=1 

160  CALL   DEFINE(Relay.Act.01,I$, I.J.K.L) 

161  IF    PCIB.ERRO0    THEN    ERROR    PCIB.BASERR 

16  2    I $= "Relay. Mux. 01": I=2:J=0:K=6:L=1 

163  CALL    DEFlNE(Relay.Mux.01, 1$, I ,J,K,L) 

164  IF    PCIB.ERRO0    THEN    ERROR    PCIB.BASERR 

800  1$=ENVIR0N$( "PANELS" ) + "\PANELS . EXE" 

801  CALL    L.S(I$  ) 
8  99    GOTO    2 
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900  IF  ERR=PCIB.BASERR  THEN  GOTO  903 

901  PRINT  "BASIC  error  #";ERR;"  occurred  in  line  "  ;  ERL 
90  2  STOP 

903  TMPERR=PCIB.ERR: IF  TMPERR=0  THEN  TMPERR=PCIB .  G'.BERR 

904  PRINT  "PC  Instrument  error  l";TMPERR;"  detected  at  line  " ;  ERL 

905  PRINT  "Error:  ";PCIB.ERR$ 

906  IF  LEFTS (PCIB. NAMES, 1 )< >CHRS ( 32)  THEN  PRINT  "Instrument:  ";PCIB. NAMES 

907  STOP 

908  COMMON  PCIB. DIRS , PCIB. SEG 

909  COMMON  LD . FILE , GET . MEM , PANELS , DEF . ERR 

910  COMMON  PCIB. BASERR.PCIB. ERR, PCIB. ERRS, PCIB. NAMES, PCIB. GLBERR 

911  COMMON  READ. REGISTER, READ. SELFID, DEFINE, INITIALIZE . SYSTEM , ENADLE . SYSTEM, DISA 
BLE . SYSTEM .INITIALIZE. POWER . ON , MEASURE . OUTPUT , START , HALT , ENABLE . INT . TR1 GGER .DISA 
BLE . INT . TRIGGER , ENABLE . OUTPUT , DISABLE . OUTPUT , CHECK . DONE , GET . STATUS 

912  COMMON  SET . FUNCTION , SET . RANGE , SET . MODE , WRITE . CAL , READ . CAL , STORE . CAL , DELAY , SA 
VE . SYSTEM , SET . GATETIME , SET . SAMPLES , SET . SLOPE , SET . SOURCE , ZERO . OHMS , SET . SPEED , SET . 
COMPLEMENT , SET . DRIVER , OUTPUT . NO . WA I T , ENABLE . HANDSHAKE , DI SABLE . HANDSHAKE 

913  COMMON  SET . THRESHOLD , SET . START . BIT , SET . NUM . B ITS , SET . LOGI C . SENSE , SET . FREOUENC 
Y , SET . AMPLITUDE . SET . OFFSET , SET . SYMMETRY , SET . BURST . COUNT . AUTOSCALE , CALI BRATE , SET . 
SENS ITI V I TY , SET . VERT . OFFSET , SET . COUPL I NG , SET . POLAR I TY , SET . SWEEPSPEED 

914  COMMON  SET  .  DELAY  .  SET  .  TR I G  .  SOURCE  .  SET  .  TRI G  .  SLOPE  ,  SET  .  TR I G  .  LEVEL  ,  SET  .  TR I G  .  MODE 
, GET . S I NGLE . WF , GET . TWO . WF , GET . VERT . INFO , GET . TIMEBASE . INFO . GET . TRI G . INFO , CALC . WFV 
OLT , CALC . WFTIME , CALC . WF . STATS , CALC . RI SETIME , CALC . FALLTIME , CALC . PERIOD 

915  COMMON  CALC . FREQUENCY , CALC . PLUSWIDTH , CALC . MINUSWIDTH , CALC . OVERSHOOT , CALC . PRE 
SHOOT , CALC . PK . TO . PK , SET . TI MEOUT , SCOPE . START , MEASURE . S I NGLE . WF , MEASURE . TWO . WF , OPE 
N . CHANNEL , CLOSE . CHANNEL 

916  COMMON  FREQUENCY , AUTO . FREQ , PERIOD , AUTO . PER , INTERVAL , RATIO , TOTALI ZE , R100MI LLI 
, Rl , RIO , R100 , R1KILO , Rl OMEGA , Rl OOMEGA , CHAN . A , CHAN . B , POSITIVE , NEGATIVE , COMN , SEPARA 
TE , DCVOLTS , ACVOLTS , OHMS , R200MI LLI , R2 , R2  0 , R200 , R2KILO , R20KI LO , R200KI LO 

917  COMMON  R2MEGA , R20MEGA , AUTOM , R2 . 5 , Rl 2 . 5 , POSITIVE , NEGATI VE , TWOS , UNSIGNED , OC , TT 
L.RO.R1.R2.R3  ,R4,R5.R6,R7  ,R8  ,  R9  ,  RIO  ,  Rll .  R12  ,  R13  ,  R14 ,  Rl  5  ,  R16  ,  SINE  ,  SQUARE  ,  TRI  ANGLE 
,CONTINUOUS,GATED,BURST,R10NANO,R100NANO,R1MICRO.R10MICRO.R100MICRO 

918  COMMON  R1MILLI .R10MILLI , R100MILLI ,R1 , RIO , R20NANO , R200NANO , R2MICPO , R20MICRO , R 
200MICRO.R2MILLI , R20MILLI , R200MILLI , R2 , R20 . R50NANO , R500NANO , R5MICRO , R50MICRO , R50 
OMICRO , R5MILLI , R50MILLI .R500MILLI , R5 , R50 , CHAN. A . CHAN. B , EXTERNAL. POSITIVE 

919  COMMON  NEGATIVE, AC, DC, TRIGGERED, AUTO. TRIG, AUTO. LEVEL, X] . X10 , STANDARD . AVERAGE 

920  COMMON  DMM .  01 ,  Func  .  Gen  .  01 ,  Scope  .  01  .Counter.  01,  Dip, .  In  .  01 .  Dig  .Out .  01  .Re  lay.  Act 
.01, Re lay. Mux. 01 

999  'End  PCIB  Program  Shell 

1000  REM  This  step  initialzes  the  HP  system 
1010  CLS 

1020  OPTION  BASE  1 

1030  DIM  P(5),PA(50,5),PP(50,5) ,XPT(50) ,YPT(50) ,X(50) ,Y(50) , YAW (50) 

1040  REM 

1050  CALL  INITIALIZE. SYSTEM! PGMSHEL.HPC ) 

1060  REM 

1070  REM  SET  FUNCTIONON  THE  'DMM'  ,  'RELAY  MUX  ,  'RELAY  ACTUATOR' 

1080  REM 

1090  CALL  SET. FUNCTION (DMM. 01, DC VOLTS) 

1100  CALL  SET. RANGE (DMM. 01, AUTOM) 

1110  CALL  DISABLE. INT. TRI GGER (DMM. 01) 

1120  CALL  ENABLE. OUTPUT (RELAY. MUX. 01) 

1130  CALL  ENABLE. OUTPUT (RELAY. ACT. 01  ) 

1140  REM  *************  PROGRAM  TRAVERSE  ******************* 

1150  REM 

1160  REM      OPEN  THE  COM  PORT  AND  INITIALIZE  THE  MOTOR  SETTINGS 

1170  OPEN  "coml: 1200, n, 8, 1 , rs , cs , ds , cd"  AS  II 

1180  REM  SET  MOTOR  DEFAULT  VALUES 

1190  DATA  2000,2000,2000,2,2,2,0.000125,0.000125,0.000125  , 
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1200  READ  V1.V2.V3.R] ,R2,R3,C1,C2,C3 

1210  REM  DEFINE  CHARACTERS  FOR  DATA  REDUCTION  ALGORITHM 

1220  RN2$="RENAME  A : RAW . DAT  " 

1230  HEAD1S  =  "  #     X      Y       PI        P2        P3        P4        P5       VA 

W   " 

1240  FORMATS =  "II   II. II  11.11   HI. ill   I  *  *  .  1*1   1*1.1*1   H  I  .  R I «   111.111   III 

.11" 

1250  PRINT 

12  60  PRINT  "*******************************************" 

1270  PRINT  "**  USER  MUST  SELECT  'CAPS  LOCK'  FUNCTION  *  +  " 

12  80  PRINT  "*******************************************" 
1290  REM  DISPLAY  MOTOR  DEFAULT  SETTINGS 

1300  PRINT  "  ♦********************♦**♦**♦****+*+" 

1310  PRINT  "  INITIALIZED  VALUES  FOR  ALL  MOTOR  SETTINGS:" 

1320  PRINT  "  VELOCITY  =  1000  STEPS/SEC" 

1330  PRINT  "  RAMPfMOTOR  ACCELERATION)  =  2   (6000  STEPS/SEC" 2  ) " 

1340  PRINT  "  DEFAULT  INCREMENTAL  UNITS  ARE  INCHES" 

13  50  PRINT  "  ***********************************" 

1360  PRINT 

1370  PRINT  "NOTE!!  USE  MANUAL  CONTROL  TO  INITIALIZE  PROBE  POSITION  BEFORE" 

1380  PRINT  "  SELECTING  COMPUTER  CONTROLLED  MOVEMENT. 

1390  PRINT 

1400  INPUT  "MANUAL  CONTROL  OR  COMPUTER  CONTROL  ( ENTER  'MAN'  or  'CP')";CON$ 

1410  IF  CON$="CP"  THEN  3490 

1420  REM   OPTION  TO  CHANGE  DEFAULT  SETTINGS  OF  VELOCITY  OR  ACCELERATION  RAMP 

1430  PRINT 

1440  PRINT 

1450  PRINT  "  DO  YOU  WANT  TO  CHANGE  THE  VELOCITY  OR  ACCELERATION  RAMP" 

1460  PRINT  "         DEFAULT  SETTINGS?  (Y  or  N)" 

1470  PRINT 

1480  PRINT  "IF  'NO',  THIS  PROGRAM  WILL  THEN  LET  YOU  DEFINE  THE" 

1490  PRINT  "DISTANCE  YOU  WANT  TO  MOVE  (IN  INCHES).  IF  'YES'." 

1500  PRINT  "YOU  CAN  CHANGE  ANY  OR  ALL  OF  THE  DEFAULT  SETTINGS  FOR  ANY  MOTOR." 

1510  PRINT  . 

1520  PRINT 

153  0  PRINT 

1540  INPUT  "DO  YOU  WANT  TO  CHANGE  ANY  OF  THE  DEFAULT  SETTINGS?  (Y  or  N ) " ; D$ 

1550  IF  D$="Y"  THEN  1590 

1560  IF  D$="N"  THEN  2220 

1570  REM 

1580  REM    ****  OPERATOR  SELECTED  MOTOR  VARIABLES  ***** 

1590  PRINT 

1600  PRINT 

1610  INPUT  "WHICH  DEFAULT  VALUE?  ( ENTER  ' 1 ' FOR  VELOC  OR  '2'  FOR  ACCEL  RAMP ) " ; L 

1620  ON  L  GOTO  1690,1^30 

1630  PRINT  "DO  YOU  WANT  TO  CHANGE  THE  DEFAULT  VELOCITY?  (Y  OR  N)" 

1640  INPUT  V$ 

1650  IF  V$="Y"  THEN  1690 

1660  PRINT  "DO  YOU  WANT  TO  CHANGE  THE  DEFAULT  ACCELERATION  RAMP?  (Y  or  N)" 

1670  IF  R$="Y"  THEN  1990 

1680  IF  R$="fo"  THEN  1450 

1690  PRINT 

17  00  PRINT 

1710  INPUT  "WHICH  MOTOR  VELOCITY  DO  YOU  WISH  TO  CHANGE?  (1,2.  or  3)"; J 

1720  ON  J  GOTO  17  30,1830,1880 

1730  PRINT 

1740  PRINi 

1750  INPUT  "ENTER  DESIRED  VELOCITY  OF  MOTOR  II"; VI 

1760  PRINT 

1770  PRINT 

1780  PRINT 

1790  PRINT  "DO  YOU  WANT  TO  CHANGE  VELOCITY  OF  ANOTHER  MOTOR?  (Y  OR  N>" 
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1800  INPUT  V$ 

1810  IF  V$="Y"  THEM  1690 

1820  IF  V$="N"  THEN  1<i30 

1830  PRINT 

1840  PRINT 

1850  INPUT  "ENTER  DESIRED  VELOCITY  OF  MOTOR  2";V2 

1860  PRINT 

18  70  GOTO  17  80 
1880  PRINT 
1890  PRINT 

1900  INPUT  "ENTER  DESIRED  VELOCITY  OF  MOTOR  H3";V3 
1910  PRINT 
1920  GOTO  1780 
1930  PRINT 
191*0    PRINT 

19  50  INPUT  "WHICH  MOTOR  ACCEL  RAMP  DO  YOU  WANT  TO  CHANGE?  (1.  2.  or  3)";K 
1960  ON  K  GOTO  1970,2060,2120 

1970  PRINT 

1980  PRINT 

1990  INPUT  "ENTER  DESIRED  ACCELERATION  RAMP  OF  MOTOR  *1";R1 

2000  PRINT 

2010  PRINT 

2020  PRINT  "DO  YOU  WANT  TO  CHANGE  THE  ACCEL  RAMP  OF  ANOTHER  MOTOR?  (Y  or  N>?" 

2030  INPUT  RM$ 

2040  IF  RM$="Y"  THEN  1930 

2050  IF  RM$="N"  THEN  1450 

2060  PRINT 

2070  PRINT 

2080  INPUT  "ENTER  DESIRED  ACCELERATION  RAMP  OF  MOTOR  *2" ;R2 

2090  PRINT 

2100  PRINT 

2110  GOTO  2000 

2120  PRINT 

2130  PRINT 

2140  INPUT  "ENTER  DESIRED  ACCELERATION  RAMP  OF  MOTOR  «3";R3  i 

2150  PRINT 

2160  PRINT 

2170  GOTO  2000  , 

2180  REM 

2190  REM  DEFINE  DISTANCE  TO  MOVE  MOTOR 

2200  PRINT 

2210  PRINT 

2220  PRINT 

2  2  30  REM  INITIALIZE  MOTOR  INCREMENTS  TO  ZERO 

221*0    11  =  0 

2250  12=0 

2260  13=0 

2270  PRINT 

2  2  80  PRINT  "  ************************************************************* 

2290  PRINT  "  **       DEFINE  WHICH  MOTOR  YOU  WANT  TO  MOVE  ** 

2300  PRINT  "  **  ** 

2310  PRINT  "  **         NOTEIM   A  POSITIVE  ('+')  INCREMENT  TO  A  MOTOR    ** 

2320  PRINT  "  **       MOVES  TRAVERSER  AWAY  FROM  THAT  PARTICULAR  MOTOR     *♦ 

2330  PRINT  "  **  ** 

2340  PRINT  "  **   --  MOTOR  HI  MOVES  THE  PROBE  UPSTREAM  AGAINST  THE  FLOW   *  + 

2350  PRINT  "  **   --  MOTOR  H 2  MOVES  THE  PROBE  TOWARD  THE  ACCESS  WINDOW    ** 

2360  PRINT  "  **   --  MOTOR  *  3  MOVES  THE  PROBE  VERTICALLY  DOWNWARD         ** 

2  3  70  PRINT  "  ************************************************************* 

2380  PRINT 

2390  PRINT 

2400  INPUT  "WHICH  MOTOR  DO  YOU  WANT  TO  MOVE?  (1.2,  or  3)";L 

2410  ON  L  GOTO  2420,2680,2970 
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2420 
2430 
2440 
2450 
2460 
2470 
2480 
2490 
2500 
2510 
2520 
2530 
2540 
2550 
2560 
2570 
2580 
2590 
2600 
2610 
2620 
2630 
2640 
2650 
2660 
2670 
2680 
2690 
2700 
2710 
2720 
2730 
2740 
2750 
2760 
2770 
2780 
2790 
2800 
2810 
2820 
2830 
2840 
2850 
2860 
2870 
2880 
2890 
2900 
2910 
2920 
2930 
2940 
2950 
2960 
2970 
2980 
2990 
3000 
3010 
3020 
3030 


PRINT 

PRINT 

PRINT  "HOW  FAR  DO  YOU  WANT  TO  MOVE  MOTOR  #1?" 

PRINT  "  *********  (ENTER  DISTANCE  IN  INCHES)  *********** 

INPUT  II 

PRINT 

PR I NT "  *********************************" 

PRINT 

PRINT  "SUMMARY  OF  OPERATOR  INPUTS:" 

PRINT  "         MOTOR  #1    VELOCITY  =  ":V1 

PRINT  "  ACCELERATION  RAMP  =  ";R1 

PRINT  "  INCREMENTAL  DISTANCE  =  "; II ;" INCHES" 

PR I NT "  *********************************" 

PRINT  "DO  YOU  WANT  TO  CHANGE  ANY  OF  THESE  VALUES?  (Y  or  N),f 

PRINT 

PRINT  "ENTER  'N'  TO  START  MOTOR  MOVEMENT.   ENTER  'Y'  TO  RETURN" 

PRINT  "TO  VARIABLE  SELECTION  SUBROUTINE." 

INPUT  V$ 

IF  V$="Y"  THEN  1430 

GOSUB  3410 

PRINT 

PRINT  "DO  YOU  WANT  TO  MOVE  ANOTHER  MOTOR  ALSO?  (Y  or  N)?" 

INPUT  C$ 

IF  C$="Y"  THEN  2220 

IF  C$="N"  THEN  3260 

PRINT 

PRINT 

PRINT  "HOW  FAR  DO  YOU  WANT  TO  MOVE  MOTOR  #2?" 

PRINT  "  *********  (ENTER  DISTANCE  IN  INCHES)  *********" 

INPUT  12 

PRINT 

PRINT 

REM  DISPLAY  OPERATOR  SELECTED  MOTOR  VARIABLES 

PRINT"  *********************************" 

PRINT 

PRINT  "SUMMARY  OF  OPERATOR  INPUTS:" 

PRINT  "         MOTOR  #2    VELOCITY  =  ";V2 

PRINT  "  ACCELERATION  RAMP  =  ";R2' 

PRINT  "  INCREMENTAL  DISTANCE  =  "; 12 ; "INCHES" 

PR I NT "  *********************************" 

PRINT 

PRINT 

PRINT  "DO  YOU  WANT  TO  CHANGE  ANY  OF  THESE  VALUES?  (Y  or  N)" 

PRINT 

PRINT  "ENTER  'N'  TO  START  MOTOR  MOVEMENT.   ENTER  'Y'  TO  RETURN" 

PRINT  "TO  VARIABLE  SELECTION  SUBROUTINE." 

INPUT  V$ 

IF  V$="Y"  THEN  1430 

GOSUB  3  410 

PRINT 

PRINT  "DO  YOU  WANT  TO  MOVE  ANOTHER  MOTOR  ALSO?  (Y  or  N)?" 

INPUT  C$ 

IF  C$="Y"  THEN  2220 

IF  C$="N"  THEN  3260 

PRINT 

PRINT 

PRINT  "HOW  FAR  DO  YOU  WANT  TO  MOVE  MOTOR  ft  3  ? " 

PRINT  "  *********  (ENTER  DISTANCE  IN  INCHES)  ********** 

INPUT  13 

PRINT 

PRINT 

REM  DISPLAY  OPERATOR  SELECTED  MOTOR  VARIABLES 
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301*0 
3050 
3060 
3070 
3080 
3090 
3100 
3110 
3120 
3130 
3140 
3150 
3160 
3170 
3180 
3190 
3200 
3210 
3220 
3230 
32t*0 
3250 
3260 
3270 
3280 
3290 
3300 
3310 
3320 
3330 
3340 
3350 
3360 
3370 
3380 
3390 
31*00 
3kl0 
3U20 
3430 
3Ui*0 
3U50 
31*60 
31*70 
3480 
3490 
3500 
3510 
3520 
3530 
3540 
3550 
3560 
3570 
3580 
3590 
3600 
3610 
3620 
3630 
3640 
3650 
3660 
3670 
3680 
3690 
3700 


I 


PRINT"  *********************************" 

PRINT 

PRINT  "SUMMARY  OF  OPERATOR  INPUTS:" 

PRINT  "         MOTOR  It  3    VELOCITY  =  ";V3 

PRINT  "  ACCELERATION  RAMP  =  ";R3 

PRINT  "  INCREMENTAL  DISTANCE  =  " ; I  3 ; " INCHES" 

PR  I  ITT 

PR I NT "  *********************************" 

PRINT 

PRINT 

PRINT  "DO  YOU  WANT  TO  CHANGE  ANY  OF  THESE  VALUES?  (Y  or  N>" 

PRINT 

PRINT  "ENTER  'N'  TO  START  MOTOR  MOVEMENT.   ENTER  'Y'  TO  RETURN' 

PRINT  "TO  VARIABLE  SELECTION  SUBROUTINE." 

INPUT  V$ 

IF  V$="Y"  THEN  1430 

GOSUB  3  410 

PRINT 

PRINT 

INPUT  "DO  YOU  WANT  TO  INPUT  ANOTHER  MANUAL  MOTOR  MOVEMENT  (Y  or  N)";M$ 

IF  M$="Y"  THEN  2210 

PRINT 

PRINT  "DO  YOU  WANT  TO  INPUT  COMPUTER  CONTROLLED  MOTOR  MOVEMENT?" 

PRINT  "        *********  NOTE! !l  *********  " 

PRINT  "  ALL  PREVIOUS  MOTOR  INCREMENT  INPUTS  HAVE  BEEN  ZEROIZED." 

PRINT  "PROGAM  WILL  LET  YOU  CHOOSE  MANUAL  OR  CP-CONTROLLED  MOVEMENT." 

PRINT  "*****  (IF  'NO'.  THE  PROGRAM  WILL  END).  *****" 

PRINT 

INPUT  "DO  YOU  WANT  COMPUTER  CONTROLLED  MOTOR  MOVEMENT  (Y  or  N ) " ; N$ 

IF  N$="Y"  THEN  3500 

PRINT 

PRINT 

PRINT 

PRINT  "     ***************************** 

PRINT  "         THE  PROGRAM  HAS  ENDED." 
PRINT  "     ****************************" 

END 

REM    *******  MOTOR  MOVEMENT  SUBROUTINE  ********* 

PRINT  II,  "&"   :PRINT  #1,  "E" : "Cl=" ; CI ; " :C2=" ;C2 ; " : C3=" ; C3 

PRINT  HI,  "Il=": I1;":V1=";V1;":R1=";R1;  , 

PRINT  HI,  ":  I2=";  12;  ":V2  =  ";V2,-":R2  =  ";R2 

PRINT  HI,  "I3="; 13  ;  "  :  V3  =  " ; V3 ; ":R3=";R3; ":@" 

RETURN 

R  EM    ******************************************** 

REM    ******************************************** 

PR  I  ITT 

REM  *******  COMPUTER  CONTROLLED  MOVEMENT  ******* 

PRINT 

PRINT  "THE  PRESSURE  DATA  WILL  BE  WRITTEN  TO  FILES  ON  DRIVE  'A'  " 

PRINT 

PRINT  "YOU  WILL  BE  ASKED  TO  INPUT  FILE  NAMES  FOR  THESE." 

PRINT 

INPUT  "IS  A  FORMATTED  DISK  IN  DRIVE  'A'?   PRESS  'ENTER'  TO  CONTINUE" ;D$ 

PRINT 

PRINT 

PRINT 

PRINT  " 

PRINT  " 

PRINT  " 

PRINT  " 

PRINT  " 

PRINT  " 

PRINT  " 

PRINT  " 

PRINT 

REM  SET  INITIAL  MOVEMENT  DISTANCE  AND  NUMBER  OF  DATA  POINTS  TO  ZERO 

HT  =  0 


****************************************" 

**  NOTE  ! ! !                 **" 

**  COMPUTER  CONTROLLED  MOVEMENT        **" 

**  IS  PROGRAMMED  WITH  A           **" 

**  DEFAULTED  NEGATIVE  MOTOR  INCREMENT  **" 

**  (i.e.  MOTOR  H3  WILL  MOVE  UPWARD    **" 

**  BY  ENTERING  A  (+)   DISTANCE).      **" 
****************************************" 
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3710  WD=0 

3720  DIST=0 

3730  XPT=0 

37U0    YPT  =  0 

3750  N=0 

3760  PRINT 

3770  PRINT 

3780  INPUT  "WHAT  IS  THE  DIMENSION  (  X  ,  Y  )  ( IN  INCHES)  THAT  YOU  WANT  TO  MEASURE 

."  ;WD,HT 

3790  PRINT 

3800  INPUT  "WHAT  IS  THE  STEP  (IN  INCHES)  THAT  YOU  WANT  TO  MOVE.";DIST 

3810  YPT=INT(HT  /DIST)  +  1 

3820  XPT=INT(WD  /DIST)+  1 

3830  N=XPT*YPT 

3840  PRINT 

38  50  PRINT  "THERE  ARE   " ; XPT ; "  *  ";YPT;"  =  ";N;"  POINTS  TO  BE  MEASURED  " 

3  8  60  PRINT 

3870  INPUT  "ARE  THE  NUMBER  OF  POINTS  IS  OK.(Y  OR  N ) " ; C$ 

3880  IF  C$="N"  THEN  3780 

3890  CLS 

3900  N  =  XPT 

3910  IF  (N  <  1)  OR  (N  >  99)  GOTO  3780 

3920  REM  ***  GENERATING  STRING  STRING  SEGMENTS  FOR  DATA  FILE  NAMES 

3930  B$  =  MID$(STR$( 1)  ,  2):  REM  **  STRING  NUMBER  "1" 

3940  E$  =  MID$(STR$(N)  ,  2):  REM  **  ENDING  STRING  NUMBER  "N" 

3950  X$  =  "XXXXXX" 

3960  EX$  =  ".DAT" 

3970  CLS 

3980  PRINT  "DATA  FILES  WILL  BE  INCREMENTED  FROM:" 

3990  PRINT 

4000  PRINT  (X$  +  B$  +  EX$);  "   To   ";  (X$  +  E$  +  EX$ ) 

1*010    PRINT 

4020  PRINT 

4030  INPUT  "ENTER   DATA  FILE  NAME  (6  CHARACTERS  MAX  --  NO  EXTENSION )": F2$ 

4040  PRINT 

4050  PRINT 

4060  IF  LEN(F2$)  >  6  OR  LEN(F2$)  <  1  GOTO  4030 

4070  CLS 

4080  PRINT  N;  "DATA  FILES  WILL  BE  GENERATED  AND  INCREMENTED' AS  FOLLOWS:" 

4090  PRINT 

4100  PRINT 

4110  PRINT  (F2$  +  B$  +  EX$  )  ;  "   To   ";  (F2$  +  E$  +  EX$ ) 

4120  PRINT 

4130  PRINT 

4140  INPUT  "ARE  THE  NUMBER  OF  POINTS  AND  FILE  NAMES  OK . ( Y  OR  N)";  C$ 

4150  IF  C$  =  "N"  GOTO  3780 

4160  IF  C$  =  "Y"  GOTO  4180 

4170  GOTO  4140 

4180  CLS  j 

4190  PRINT 

4200  PRINT 

4210  REM  SET  INITIAL  POSITION  DATA 

4220  X(1)=-DIST 

4230  Y(1)=-DIST 

4240  FOR  IX=2  TO  XPT+1 

4250  X(IX)=0 

4260  NEXT  IX 

4270  FOR  JY=2  TO  YPT+1 

4280  Y(JY)=0 

4290  NEXT  JY 

4  300  FOR  1=1  TO  XPT 
4302  11=0 

4304  12=0 
4306  13=0 
4  310  FOR  J 


4  3  20  REM 


1  TO  YPT 
MOTOR   CP- CONTROLLED  MOTOR  MOVEMENT 
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4330  11=0 

k3i*0    12-0 

4350  13=0 

4360  REM  EACH  POINT  TAKE  10  TIMES  READINGS 

4370  X( I+1)=X( I )+DIST 

4380  XPT( J)=X( 1+1) 

4390  Y( J+l )=Y( J)+DIST 

4400  YPT(J)=Y(J+1 ) 

4405  INPUT  "  ADJUST  THE  WHEEL  TO  MAKE  THE  P2  =P3 . INPUT* THE  YAW  ANGLE" ; YAW( J ) 

4408  PRINT 

4410  INPUT  "  PRESS  'ENTER'  TO  START  THE  MEASUREMENT" ; MOVES 

4420  REM 

44  30  REM  READ  FIVE  CHANNELS  AND  DISPLAY  THE  DATA 

4440  REM 


4450 

STEPPER= 

4 

4460 

SWITCH  = 

3 

4470 

HOMER=8 

4480 

DELAY 1  = 

.1 

4490 

DELAY2  = 

1 

4500 

REM  SET 

THE 

S, 

V 

PORT 

TO 

14 

4510 

FOR  IL=1 

TO 

3 

4520  THYME  =  TIMER 

4  5  30  CALL  OUTPUT ( RELAY . ACT . 01 , STEPPER ) 

4540  CHKTIME  =  TIMER 

4550  IF  CHKTIME  <  (THYME  +  DELAY1)  GOTO  4540 

4560  CALL  OPEN . CHANNEL ( RELAY . ACT. 01 , SWITCH ) 

4570  CLS 

4580  NEXT  IL 

4590  PRINT 

4600  PRINT  "  NOW  IS  POINT  ";J 

4  610  REM  START  MEASURE  FROM  PORT  4  TO  PORT  8 

4620  FOR  JJ=1  TO  5 

4630  CALL  OUTPUT ( RELAY . ACT . 01 , STEPPER ) 

4640  CHKTIME  =  TIMER 

4650  IF  CHKTIME  <  (THYME  +  DELAY2 )  GOTO  4640 

4  660  REM  EACH  PORT  SAMPLE  10  TIMES 

4670  FOR  11=1  TO  10 

4  6  80  ROUT=l 

4690  CALL  OUTPUT ( RELAY . MUX . 01 , ROUT ) 

4700  CALL  MEASURE ( DMM . 01, VOLTS) 

4710  PA(II ,JJ)=VOLTS 

4  7  20  NEXT  II  I 

4730  CALL  OPEN . CHANNEL! RELAY . ACT . 01 , SWITCH ) 

4740  IF  JJ=5  THEN  4760 

"4  7  50  NEXT  JJ 

4  7  60  REM  HOME  THE  S.V  PORT  TO  0  4  8 

4770  CALL  OUTPUT ( RELAY . ACT . 01 , HOMER ) 

4780  CALL  OPEN . CHANNEL ( RELAY . ACT . 01 , HOMER ) 

4790  REM 

4800  REM  DISPLAY  THE  SAMPLE  DATA 

4  810  REM 

4820  PRINT  HEAD1S 

4830  FOR  IS=  1  TO  10 

4840  PRINT  USING  FORMAT*  ;  IS . XPT ( J  )  . YPT( J ) . PA ( IS . 1 ) , PA ( IS , 2 ) . PA ( IS . 3 ) . PA ( IS . 4 ) , PA 

( IS, 5)  ,YAW( J ) 

4  8  50  NEXT  IS 

4860  REM 

4870  REM  AVERAGE  THE  DATA 

4  8  80  REM 

4890  FOR  JA  =  1  TO  5 

4900  TOTAL  =  0 

4910  FOR  IA  =  1  TO  10 

4920  TOTAL  =  TOTAL  +  PA(IA.JA) 

49  30  NEXT  I A 

4940  AVERAGE  =  TOTAL  /10 

4950  P(JA)=AVERAGE 

4960  NEXT  JA 
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(.970  PRINT 

4980  PRINT  "THE  AVERAGES  ARE:  " 

5000  PRINT  HEAD1S 

5010  FOR  JD=1  TO  5 

5020  PP( J,JD)=P(JD) 

5030  NEXT  JD 

5040  PRINT  USING  FORMATS  ;  J  ,  XPT<  J  )  .  YPT  (  J  )  ,  PP(  J  .  1 )  .  PP(  J  .  2  )  ,  PP(  J  ,  3  >  ,  PP<  J  .  4  )  ,  PP(  J  ,  5  ) 

,YAW( J) 

5045  PRINT 

5048  PRINT  USING  "THE  NULLING  ERROR  IS  +| .  Ml  I"  ;PP(  J,  3  ) -PP(  J,  2  ) 

50(*9  PRINT 

5050  PRINT  "DO  YOU  WANT  RE-MEASURE  AGAIN  (Y  /  N)" 

5060  PRINT 

5062  PRINT  "IF  'Y'  WILL  RE-SAMPLE  AGAIN." 

5061*  PRINT 

5070  INPUT  "IF  'N'  WILL  MOVE  THE  TRAVERSER  STEP  UPWARD  (WAIT  7  SEC  )";C$ 

5075  PRINT 

5080  IF  C$  =  "Y"  THEN  4405 

5082  IF  C$="N"  THEN  5090 

5084  GO  TO  5070 

5090  IF  J=YPT  THEN  5160 

5100  REM 

5110  REM  MOVE  THE  TRAVERSER  STEP  UPWARD. 

5120  REM 

5130  I3=-DIST 

5140  GOSUB  3410 

5150  NEXT  J 

5160  REM***  STORE  DATA  BEFORE  NEXT  SAMPLE*** 

5170  OPEN  "A:\RAW.DAT"  FOR  OUTPUT  AS  12 

5180  PRINT  #2  ,HEAD1$ 

5190  FOR  ID=1  TO  YPT 

5200  PRINT  12  .USING  FORMATS ; ID , XPT( ID ) , YPT( ID ) . PP( ID, 1 ) , PP ( ID , 2 ) , PP( ID , 3 ) , PP( ID 

,4) ,PP( ID, 5) ,YAW( ID) 

5  210  NEXT  ID 

5220  CLOSE  12  '  ' 

5230  REM  ***  GENERATING  INCREMENTED  DATA  FILE  NAME 

5240  IF  (I  >  10)  OR  (I  =  10)  THEN  1$  =  MIDI  ( STR$  <  I  )  ,  2) 

5250  IF  (I  <  10)  THEN  1$  =  ( MID$ ( STR$ ( 0 ) .  2)  +  MID$ ( SIRS ( I ) .  2))  i 

5260  FI2$  =  (F2$  +  1$  +  EXS ) 

52  70  PRINT  ' 

5260  PRINT  "  WRITING  DATA  FILE  " ;  FI2S 

5290  DF2$=RN2$+FI2$ 

5300  REM  **  RENAME   DATA  FILE 

5310  SHELL  DF2$ 

5320  REM 

5  3  30  REM  MOVE  THE  TRAVERSER  TO  THE  NEXT  SAMPLE  POSITION 

5  34  0  REM 

5350  PRINT 

5360  IF  I=XPT  THEN  5430 

5370  INPUT  "THEN  PRESS  'ENTER'  FOR  NEXT  COLUMN  SAMPLE (  90  SEC)  ";MOVE$ 

5390  I2=-DIST 

5400  I3=HT 

5410  GOSUB  3410 

5420  NEXT  I 

5430  CLS 

5440  PRINT  "ALL  MOVEMENTS  COMPLETE" 

5450  PRINT 

5460  PRINT 

5470  PRINT  "YOU  WANT  TO  REPOSITION  TRAVERSER  FOR  ANOTHER  MOVEMENT  (Y  OR  N)?" 

5480  PRINT 

5490  PRINT  "IF  'Y',  THE  PROGRAM  WILL  TAKE  YOU  TO  MANUAL  CONTROL  SUBROUTINE." 

5500  PRINT  "IF  'N',  THE  PROGRAM  WILL  END." 

5510  PRINT 

5520  INPUT  "ANOTHER  MOVEMENT" ;R$ 

5530  IF  R$  =  "Y"  THEN  1370 

5540  IF  R$  =  "N"  THEN  3370 
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APPENDIX  B.  CALP  PROGRAM 


1  DEF  SEG:  CLEAR  ,  &I1FEOO  :  GOTO  4  "Begin  PCIB  Program  Shell 

2  GOTO  1000  '  User  program 

3  GOTO  900  '  Error  handling 

4  I=SHFE00  *  Copyright  Hewlett-Packard  1984,1985 

5  PCIB.DIRS=ENVIRON$("PCIB" ) 

6  1S=FCIB.DIR$+"\PCIBILC.BLD" 

7  BLOAD  I$,l 

8  CALL  I  (PCIB.  DIRS,  I7..J7.)  :  PCIB.  SEG=I% 

9  IF  J7„=0  THEN  GOTO  13 

10  PRINT  "Unable  to  load."; 

11  PRINT  "    (Error  l";J%;")" 

12  END 

13  ' 

14  DEF  SEG=PCIB.SEG:O.S=5:C.S=10: I .V=15 

15  I .C=20:L.P=25:LD.FILE=30 

16  GET.MEM=35:L.S=4Q:PANELS=45:DEF.ERR=50 

17  PCIB.ERRS=STRING$(64, 32)  :  PCIB . NAMES=STRINGS ( 16 , 32 ) 

18  CALL  DEF. ERR(PCIB. ERR, PCIB. ERRS, PCIB. NAMES, PCIB. GLBERR)  :  PCIB . BASERR= 255 

19  ON  ERROR  GOTO  3 

20  J=-l 

21  IS=FCIB.DIRS+"\PCIB.SYN" 
2  2  CALL  O.S( IS) 

23  IF    PCIB.ERRO0    THEN    ERROR    PCIB.BASERR 

24  1=0 

2  5  CALL  I .V( I , READ. REG I STER , READ. SELFID, DEFINE, INITIALIZE. SYSTEM) 
26  IF  PCIB.ERROO  THEN  ERROR  PCIB.BASERR 

2  7  CALL  I .V( I .ENABLE. SYSTEM, DISABLE. SYSTEM, INITIALIZE, POWER. ON) 

28  IF  PCIB.ERROO  THEN  ERROR  PCIB.BASERR 

29  CALL  I .V( I .MEASURE, OUTPUT, START, HALT) 

30  IF  PCIB.ERROO  THEN  ERROR  PCIB.BASERR  I 

31  CALL  I .V( I .ENABLE. INT. TRIGGER, DISABLE. INT . TRIGGER , ENABLE. OUTPUT , DISABLE . OUTPU 
T) 

32  IF  PCIB.ERROO  THEN  ERROR  PCIB.BASERR 

3  3  CALL  I .V( I .CHECK. DONE. GET. STATUS. SET. FUNCTION, SET. RANGE) 

34  IF  PCIB.ERROO  THEN  ERROR  PCIB.BASERR 

35  CALL  I .V(I , SET. MODE, WRITE. CAL, READ. CAL, STORE. CAD 

36  IF  PCIB.ERROO  THEN  ERROR  PCIB.BASERR 
3  7  CALL  I .V( I .DELAY, SAVE. SYSTEM, J, J) 

38  IF  PCIB.ERROO  THEN  ERROR  PCIB.BASERR 

3  9  1  =  1 

40  CALL  I .V( I , SET. GATETIME, SET. SAMFLES, SET. SLOPE, SET. SOURCE) 

41  IF  PCIB.ERROO  THEN  ERROR  PCIB.BASERR 

4  2  CALL  I . C( I .FREQUENCY, AUTO. FREO , PERIOD, AUTO. PER) 
4  3  IF  PCIB.ERROO  THEN  ERROR  PCIB.BASERR 

44  CALL  I .C( I, INTERVAL, RATIO, TOTALIZE, R100MILLI ) 

4  5    IF    PCIB.ERROO    THEN    ERROR    PCIB.BASERR 

46    CALL    I .C( I ,R1,R10,R100,R1KILO) 

4  7    IF    PCIB.ERROO    THEN    ERROR    PCIB.BASERR 

48  CALL    I .C( I .RIOMEGA.RIOOMEGA.CHAN.A.CHAN.B) 

49  IF  PCIB.ERROO  THEN  ERROR  PCIB.BASERR 

50  CALL  I .C( I .POSITIVE, NEGATIVE, COMN, SEPARATE) 

51  IF  PCIB.ERROO  THEN  ERROR  PCIB.BASERR 

52  1  =  2 

53  1  =  3 
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54  CALL  I  .V(  1  .ZERO.  OHMS,  SET.  SLEEP,  J  ,  J  ) 

55  IF  PCIB.ERROO  THEN  ERROR  PCID .  BASI5RR 

56  CALL  I .C( I , DCVOLTS.ACVOLTS.OHMS.R200MILLI ) 

57  IF  PCIB.ERROO  THEN  ERROR  PCIB.BASERR 

58  CALL  I .C( I ,R2,R20,R200,R2KILO) 

59  IF    PCIB.ERROO    THEN    ERROR    PCIB.BASERR 

60  CALL    I .C( I ,R20KILO,R200KILO,R2MEGA,R20MEGA) 

61  IF    PCIB.ERROO    THEN    ER.R0R    PCIB.BASERR 

62  CALL    I .C(I .AUTOM.R2. 5.R12.5.J) 

63  IF    PCIB.ERROO    THEN    ERROR    PCIB.BASERR 

64  1  =  4 

65  CALL  I .V( I .SET. COMPLEMENT, SET. DRIVER, OUTPUT. NO. WAIT, ENABLE. HANDSHAKE) 

66  IF  PCIB.ERROO  THEN  ERROR  PCIB.BASERR 

67  CALL  I . V ( I , DI SABLE . HANDSHAKE , SET . THRESHOLD , SET . START . BIT , SET . NUM . BI TS ) 

68  IF  PCIB.ERROO  THEN  ERROR  PCIB.BASERR 

69  CALL  I .V(I , SET. LOGIC. SENSE, J, J, J ) 

70  IF  PCIB.ERROO  THEN  ERROR  PCIB.BASERR 

71  CALL  I. C( I .POSITIVE, NEGATIVE, TWOS. UNSIGNED) 

72  IF  PCIB.ERROO  THEN  ERROR  PCIB.BAFERR 

7  3  CALL  I .C( I .OC.TTL.RO.Rl) 

74  IF  PCIB.ERROO  THEN  ERROR  PCIB.BASERR 

75  CALL  I .C(I ,R2,R3,R4,R5) 

76  IF  PCIB.ERROO  THEN  ERROR  PCIB.BASERR 

77  CALL  I.C( I .R6.R7.RR.R9) 

78  IF    PCIB.ERROO    THEN    ERROR    PCIB.BASERR 

79  CALL    I .C( I .R10.R11.R12.R13 ) 

80  IF    PCIB.ERROO    THEN    ERR.0R    PCIB.BASERR 

81  CALL    I .C(I .R14.R15.R16, J ) 

8  2  IF  PCIB.ERROO  THEN  ERROR  PCIB.BASERR 
83  1=6 

8  4  CALL  I .V(I , SET. FREQUENCY. SET. AMPLITUDE. SET. OFFSET. SET. SYMMETRY) 
85  IF  PCIB.ERROO  THEN  ERROR  PCIB.BASERR 
8  6  CALL  I .V( I .SET. BURST. COUNT, J, J.J ) 

8  7  IF  PCIB.ERROO  THEN  ERROR  PCIB.BASERR 

88  CALL  I .C( I .SINE, SQUARE. TRIANGLE, CONTINUOUS)  , 

89  IF  PCIB.ERROO  THEN  ERROR  PCIB.BASERR 

90  CALL  I .C( I .GATED, BURST. J , J) 

91  IF    PCIB.ERROO    THEN    ERROR    PCIB.BASERR 

92  1  =  7 

93  CALL  I .V( I , AUTOSCALE, CALIBRATE. SET. SENSITIVITY, SET. VERT. OFFSET) 

94  IF  PCIB.ERROO  THEN  ERROR  PCIB.BASERR 

9  5  CALL  I .V( I , SET. COUPLING, SET. POLARITY, SET. SWEEPSPEED, SET. DELAY) 
96  IF  PCIB.ERROO  THEN  ERROR  PCIB.BASERR 

9  7  CALL  I . V( I .SET. TRIG. SOURCE, SET. TRIG. SLOPE, SET. TRIG. LEVEL, SET. TRIG. MODE) 

98  IF  PCIB.ERROO  THEN  ERROR  PCIB.BASERR 

99  CALL  I .V( I , GET. S1NGLE.WF, GET. TWO. WF, GET. VERT. INFO, GET. TIMEBASE. INFO) 

100  IF  PCIB.ERROO  THEN  ERROR  PCIB.BASERR 

101  CALL  I . V( I .GET. TRIG. INFO . CALC . WFVOLT , CALC .WFTIME , CALC . WF . STATS ) 

102  IF  PCIB.ERROO  THEN  ERROR  PCIB.BASERR 

103  CALL  I .V( I , CALC. RISETIME, CALC. FALLTIME, CALC. PERIOD, CALC. FREQUENCY ) 

104  IF  PCIB.ERROO  THEN  ERROR  PCIB.BASERR 

105  CALL  I .V( I , CALC. PLUSWIDTH. CALC. MINUSWIDTH, CALC. OVERSHOOT, CALC. PRESHOOT) 

106  IF  PCIB.ERROO  THEN  ERROR  PCIB.BASERR 

10  7  CALL  I .V( I , CALC. PK. TO. PK, SET. TIMEOUT, SCOPE. START, MEASURE. S I NGLE.WF) 

108  IF  PCIB.ERROO  THEN  ERROR  PCIB.BASERR 

109  CALL  I .V( I, MEASURE. TWO. WF, J, J, J ) 
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110  IF    rCID.CTROO    THEN    ERROR    PCIB.BASERR 

111  CALL    1 .C( I .R10NANO.R100NANO.R1M1CRO.R10MJCRO) 

112  IF    PCIB.ERROO    THEN    ERROR    PCIB.PASFPR 

113  CALL    I .C( I .R100MICRO.R1MILLI.R10MILLI .R100MILLI ) 

114  IF    PCIB.ERROO    THEN    ERROR    PCIB.BASERR 

115  CALL    I .C(I ,R1,R10,R20NANO,R200NANO) 

116  IF    PCIB.ERROO    THEN    ERROR    PCIB.BASERR 

117  CALL    I .C( I .R2MICRO.R20MICRO.R200MICRO.R2MILLI ) 

118  IF    PCIB.ERROO    THEN    ERROR    PCIB.BASERR 

119  CALL    I .C( I,R2OMILLI,R200MILLI,R2.R20) 

120  IF    PCIB.ERROO    THEN    ERROR    PCIB.BASERR 

121  CALL    I .C(I,R50NANO.R500NANO,R5MICRO,R50MICRO) 

122  IF    PCIB.ERROO    THEN    ERROR    PCIB.BASERR 

123  CALL    I.Cd .R500MICRO.R5MILLI .R50MILLI .R500MILLI ) 

124  IF    PCIB.ERROO    THEN    ERROR    PCIB.BASERR 

125  CALL    I  .Cd  .R5.R50.CHAN.A.CHAN.B) 

126  IF  PCIB.ERROO  THEM  ERROR  PCIB.BASERR 

127  CALL  I .C( I, EXTERNAL, POSITIVE. NEGATIVE. AC) 

128  IF  PCIB.ERROO  THEN  ERROR  PCIB.BASERR 

129  CALL  I .C( I .DC. TRIGGERED, AUTO. TRIG. AUTO. LEVEL) 

130  IF  PCIB.ERROO  THEN  ERROR  PCIB.BASERR 

131  CALL  I .C( I .XI ,X10, STANDARD, AVERAGE) 

132  IF  PCIB.ERROO  THEN  ERROR  PCIB.BASERR 

133  1=8 

134  CALL  I. V( I, OPEN. CHANNEL, CLOSE. CHANNEL, J, J) 

135  IF  PCIB.ERROO  THEN  ERROR  PCIB.BASERR 

136  CALL  C.S 

137  IF  PCIB.ERROO  THEN  EPROR  PCIB.BASERR 

138  I$=PCIB.DIRf+"\PCIB.PLD" 

139  CALL  L.P( If ) 

140  IF  PCIB.ERROO  THEN  ERROR  PCIB.BASERR 
14  1  I$="DMM.01": I=3:J=0:K=0:L=1 

142  CALL  DEFINE(DMM. 01.1$, I, J.K.L) 

143  IF  PCIB.ERROO  THEN  ERROR  PCIB.BASERR 

144  If ="Func.Gen.01 " : I = 6 : J=0 : K= 1 : L= 1 

145  CALL    DEFINE(FUNC.GEN.01 , If , I , J.K.L) 

146  IF    PCIB.ERROO    THEN    EPROR    PCIB.BASEPR 
14  7  If = "Scope. 01" :I=7:J=0:K=2:L=1 

148  CALL    DEFINE(SCOPE.01, If . I .J.K.L) 

149  IF    PCIB.ERROO    THEN    ERROR    PCIB.BASERR 

150  If ="Counter.01": 1=1: J=0:K=3 :L=1 

151  CALL  DEFlNEf COUNTER. 01 , If , I . J.K.L) 

152  IF  PCIB.ERROO  THEN  ERROR  PCIB.BASERR 

153  If ="Dig. In. 01": 1=4: J=0:K=4:L=1 

154  CALL  DEFINEfDIG. IN.01, If , I , J.K.L) 

155  IF  PCIB.ERROO  THEN  ERROR  PCIB.BASERR 
1.5  6  If  =  "Dig.  Out.  01":  1  =  4  :  J=l:K=<i  :L=1 

157  CALL  DEFINEfDIG. OUT. 01 , If , I . J. K.L) 

158  IF  PCIB.ERROO  THEN  ERROR  PCIB.BASERR 

159  If ="Relay.Act.01":I=8:J=0:K=5:L=l 

160  CALL    DEFINE(RELAY. ACT. 01,1$ . I . J.K.L) 

161  IF    PCIB.ERROO    THEN    ERROR    PCIB.BASERR 

162  I $= "Relay. Mux. 01": I =2 : J=0 :K=6 : L= 1 

163  CALL  DEFINE(RELAY.MUX.01, 1$. I .J.K.L) 

164  IF  PCIB.ERROO  THEN  ERROR  PCIB.BASERR 

800  If =ENVIRON$ ( "PANELS" ) +"\PANELS . EXE" 

801  CALL  L.S(1$) 
899  GOTO  2 
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900  IF  ERR=PCIB.BASERR  THEN  GOTO  903 

901  PRINT  "BASIC  error  »";ERR;"  occurred  in  line  " ; ERL 

902  STOP 

90  3  TMPERR=PCIB.ERR: IF  TMPERR  =  0  THEN  TMPER£=PCIB.  GLBERP 

904  PRINT  "PC  Instrument  error  »";TMPERR;"  detected  at  line  " ; ERL 

905  PRINT  "Error:  ";PC1B.ERP$ 

906  IF  LEFTSIPCIB. NAMES, 1 )<>CHR$(32 )  THEN  PRINT  "Instrument:  ";PCIB. NAMES 

907  STOP 

908  COMMON  PCIB . DIRS . PCIB . SEG 

909  COMMON  LD . FILE , GET . MEM , PANELS , DEF . ERR 

910  COMMON  PC  I  B  .  BASERP. ,  PC  I B  .  ERR  .  PC  I  B  .  ERR$  .  PC  I B  .  NAMES  ,  PCI  B .  GLBERR 

911  COMMON  READ .REGISTER , READ . SELFID , DEFINE . INITI ALI ZE . SYSTEM , ENABLE . SYSTEM . DISA 
BLE . SYSTEM , INITI ALI ZE , POWER . ON , MEASURE , OUTPUT , START , HALT . ENABLE . INT . TRIGGER , DISA 
BLE . I NT . TRI GGER , ENABLE . OUTPUT , D I SABLE . OUTPUT , CHECK . DONE , GET . STATUS 

912  COMMON  SET . FUNCTION , SET . RANGE , SET . MODE .WRITE . CAL , READ . CAL , STORE . CAL , DELAY . SA 
VE . SYSTEM , SET . GATETIME , SET . SAMPLES , SET . SLOPE , SET . SOURCE , ZERO . OHMS , SET . SPEED , SET . 
COMPLEMENT  ,  SET  .  DP. I VER  .  OUTPUT  .NO. WAIT,  ENABLE  .  HANDSHAKE  ,  DI  SABLE  .  HANDSHAKE 

913  COMMON  SET . THRESHOLD , SET . START .BIT. SET . NUM .BITS, SET . LOGI C . SENSE , SET . FREQUENC 
Y . SET . AMPLITUDE . SET . OFFSET , SET . SYMMETPY . SET . BURST . COUNT , AUTOSCALE , CALI BRATE , SET . 
SENSITIVITY , SET . VERT . OFFSET , SET . COUPLING , SET. POLARITY , SET. SWEEPSPEED 

9 1 4  COMMON  SET . DELAY . SET . TRI G . SOURCE , SET . TRI G . SLOPE , SET . TRI G . LEVEL . SET . TRIG .MODE 
, GET . S I NGLE . WF . GET . TWO . WF , GET . VERT . I NFO . GET . TI MEBASE . INFO . GET . TRI G . INFO . CALC . WFV 
OLT.CALC.WFTI ME. CALC.WF. STATS. CALC.R I SETIME.CALC.FALLTIME.CALC. PERIOD 

915  COMMON  CALC . FREQUENCY . CALC . PLUSWI DTH , CALC . MI NUSWI DTH . CALC . OVERSHOOT . CALC . PRE 
SHOOT , CALC . PK . TO . PK , SET .TIMEOUT , SCOPE . START . MEASURE .SINGLE . WF , MEASURE . TWO . WF , OPE 
N . CHANNEL . CLOSE . CHANNEL 

916  COMMON  FREQUENCY. AUTO. FREQ, PERIOD, AUTO. PEP, INTERVAL , RATIO , TOTALI ZE , R100MILLI 
, Rl . P 10, Rl 00. R1K I LO.P IOMEGA. R100MEGA, CHAN. A, CHAN. B. POSITIVE, NEGATIVE. COMN.SEPARA 
TE.DCVOLTS.ACVOLTS.OHMS.P200MILLI , R2 , R2  0 . R200 , R2KI LO , R20KILO , R200KILO 

917  COMMON  R2MEGA , R20MEGA , AUTOM .R2.5.P12.5, POS I TI VE , NEGATI VE , TWOS , UNS I GNED . OC , TT 
L,R0,R1.P2,R3,R4,R5.P6,P7  ,  P8  ,  R9  ,  RIO  ,  Rl  1 ,  R12  ,  Rl  3  ,  R  14  .  R15  .  R16  ,  SINE  .  SQUARE  .  TRI  ANGLE 
, CONTINUOUS, GATED, BURST, R10NANO,  P.  100NANO.R1MICRO.P10MICRO.R100MICRO 

918  COMMON  R1MILLI .R10MILLI .R100M1LLI , P 1 , RIO . R20NANO , R200NANO , R2MICRO . R20MICRO . R 
200MICRO.R2MILLI .P20MILLI , R200MILLI . R2 , R20 . P50NANO , R500NANO . R5MICRO , R50MICRO , R50 
OMI CRO . R5MILLI , R50MILLI , R500M1LLI , R5 , R50 , CHAN . A , CHAN . B , EXTERNAL , POSITIVE 

919  COMMON  NEGATI VE , AC , DC , TRI GGERED , AUTO . TRI G , AUTO . LEVEL , X 1 , X 1 0 , STANDARD , AVEP AGE 

920  COMMON  DMM .01, FUNC . GEN .01, SCOPE .01, COUNTER .01, DIG. IN. 01. DIG. OUT .01. RELAY . ACT 
.01  .PELAY.M1JX.01 

999  'End  PCIB  Program  Shell 

1000  P.EM 

1010  REM  This  step  initial zes  the  HP  system 

1020  CLS 

1030  OPTION  BASE  1 

1040  DIM  P(  10)  ,PA(50,6)  ,  PP  (  50,6)  .XPTUO)  ,CAL(40) 

1050  CALL  INITIALIZE. SYSTEM ( PGMSHEL . HPC ) 

1060  REM 

1070  PEM  All  PC  devices  now  have  an  initial  state 

1080  REM  Set  function  on  the  DMM  and  Relay  MUX 

1090  REM 

1100  CALL  SET. FUNCTION (DMM. 01, DC VOLTS) 

1110  CALL  SET. RANGE (DMM. 01, AUTOM) 

1120  CALL  DISABLE.  INT.  TRIGGER  ( r>MM  .  01  ) 

1130  CALL  ENABLE. OUTPUT (RELAY. MUX. 01 ) 

1140  FORMATS-"**   t  *  .  M  H  #   *»  .  t  HI*   14. fill   Iff  .•«#!   If. Mil   H.tftfft" 

1200  FOR  1=1  TO  10 

1210  CAL( I  )=0.0 
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1220  NEXT  I 

3510  REM 

1520  REM  READ  THE  VOLTAGE  OF  48TH   CHANNEL  AND  DISPLAY  THE  DATA 

1530  REM 

151*0  PRINT  "  CHOOSE  6  POINTS" 

1550  PRINT 

1550  PRINT  "THE  CALIBRATION  WILL  BE  STORES  IN  'CAL.DAT'" 

1560  REM 

1570  REM  Begin  sampling  loop 

1580  REM 

1600  FOR  J=l  TO  1 

1610  PRINT 

1630  FOR  JJ=1  TO  6 

1631  INPUT  "INPUT  THE  CALIBRATION  PRESSURE" ; CAL( JJ  ) 
163  2  INPUT  "PRESS  'ENTER'  TO  START   MEASUREMENT" ; MOVES 
1640  FOR  11  =  1  TO  10 

1650  ROUT=l 

1660  CALL  OUTPUT (RELAY. MUX. 01, ROUT) 

1670  CALL  MEASURE (DMM. 01, VOLTS) 

1680  PA(II,JJ)=VOLTS 

1690  NEXT  II 

1700  IF  JJ  =  6  THEN  1740 

17  30  NEXT  J J 

1740  REM 

1750  REM  DISPLAY  THE  SAMPLE  DATA 

1760  REM 

17  80  FOR  IS=  1  TO  10 

1790  PRINT  USING  FORMATS  ;  IS  ,  PA  (  IS  ,  1  )  .  PA(  IS  ,  2  )  ,  PA  (  IS  ,  3  )  ,  PA  (  IS  ,  4  )  ,  PA  (  IS  ,  5  )  ,  PA  (  IS  .  6 

) 

1800  NEXT  IS 

1810  REM 

18  20  REM  AVERAGE  THE  DATA  , 
18  30  REM 

1840  FOR  JA  =  1  TO  6 

1850  TOTAL  =  0 

1860  FOR  IA  =  1  TO  10 

1870  TOTAL  =  TOTAL  +  PA(IA.JA) 

18  80  NEXT  I A 

1890  AVERAGE  =  TOTAL  /10 

1900  P(JA)=AVERAGE 

1920  NEXT  JA 

1930  PRINT 

1940  PRINT  "THE  AVERAGE  ARE:  " 

2000  FOR  JD=1  TO  6 

2010  PP(J,JD)=P(JD) 

2020  NEXT  JD 

2055  PRINT  USING  FORMATS ; J , PP( J , 1 ) , PP ( J , 2 ) , PP( J , 3  )  , PP ( J , 4  )  , PP ( J  ,  5  )  , PP ( J , 6  ) 

2070  PRINT 

2080  INPUT  "DO  YOU  WANT  RE-MEASURE  AGAIN  7  (Y  /  N ) " ; C$ 

2090  IF  CS="Y"  THEN  1580 

2101  REM***  STORE  DATA  BEFORE  NEXT  SAMPLE*** 

2102  OPEN  "A:\CAL.DAT"  FOR  OUTPUT  AS  *  2 

2106  FOR  ID=1  TO  6 

2107  PRINT  H2, USING  FORMATS ; ID. PP( J , ID ) , CAL( ID ) 

2108  NEXT  ID 

2109  CLOSE  «2 
2210  NEXT  J 
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APPENDIX  C.  CONVERT  PROGRAM 


*  THIS  PROGRAM  CONVERTS  THE  VOLTAGE  OF  TRANSDUCER  INTO  PHYSICAL  * 

*  PRESSURE, VELOCITY, YAW  ANGLE  AND  PITCH  ANGLE.  THOSE  DATA  ARE    * 

*  USED   FOR  PLOT  PROGRAM  LATER.  * 
t^t********************************************* ***************** 

CHARACTER* 12  FNAME 
CHARACTER* 12  NAME 
CHARACTER* 2  A (50) 
CHARACTER* 80  ST 


REAL  K.INTR 
INTEGER  COLS.RWS 


DATA  A/ 


01 
10 
19 
28 

37 

46 


'02 

1 11 

'20 
'29 
1  38 
'47 


WRITE  (* 
READ  (*, 
WRITE  (* 
READ  (*, 
WRITE  (* 
READ  (*, 
WRITE  (* 
READ  (*, 
WRITE  (* 
READ  (*, 
WRITE  (* 
READ  (*, 
WRITE  (* 
READ  (*, 
WRITE  (* 
READ  ( * , 
WRITE  (* 
READ  ( 
WRITE 
READ  ( 
WRITE 
READ  ( 
WRITE 
READ  ( 
WRITE 
READ  ( 
CONVERT  THE 


(A\)  ' 
15)  '  ) 
(A\)  ' 
15)  '  ) 
<A\)  ' 
A12)-' 
(A\)  ' 
F4.2) 
(A\)  ' 
F4.  2  ) 
(A\)  ' 
F3.  1) 
(A\)  ' 
F3.  1  ) 
(A\)  * 
F6.4) 
(A\)  ' 
F9.  6) 
(A\)  ' 
F9.  6) 
(  A\  )  ' 
F4  .  2  ) 
(A\  )  ' 
F5.  2) 
(A\)  ' 
' (F5. 2 ) 
PRESSURE 


09 
18 
27 
36 
'45 


OF  DATA  PTS  IN  A  COL  (UP/DOWN)  = 
FILE  NAME?  (IE  R001A2XX.DAT) 
(F4.2)  =  ' 


DTPTS 

'03'  ,  '01*'  ,  '05'  ,  '06'  ,  '07'  ,  »08' 

'  12'  ,  '  13'  ,  'Ik'   ,  '15'  ,  '16'  ,  '17' 

'  21'  .  '  22'  ,  '  23'  ,  '  24  '  ,  '  25'  ,  '26' 

'  30'  ,  '  31'  ,  '32'  ,  '33'  ,  '  34'  ,  '35' 

'  39*  ,  '40'  ,  'Ul'   ,  'k2'   ,  '1*3'  ,   '44' 

'48' , '49' , ' 50'/ 

'  I  OF  COLS  (AWAY  FROM  MSL )  = 
COLS 

'  0 
RWS 

'  DATA 

NAME 

*  PI 
)  PI 

'  PF  (F4.2)  =  ' 
)  PF 

'  TI  (F3. 1 )  =  ' 
)  TI 

'  TF  (F3.1)  =  ' 
)  TF 

1  K  (F6.4)  =  ' 
)  K 

'  SLOPE  FOR  DELTAP  ( F9 . 6 )  =  ' 
)  SLOPE 

'  INTERCEPT  FOR  DELTAP  ( F9 . 6 ) 
)  INTR 

'  QM1  FACTOR  ( F4 . 2 )  =  ' 
)  OM1FAC 

1  X  OFFSET  =  ' 
)  XOFF 

'  Y  OFFSET  =  ' 
)  YOFF 

UNIT  FROM  inHg  TO  psf 


PATM=(PI+PF)*35.  3631 

R=1716. 5 

E=0.0123 

T=(TI+TF)/2.+460 

RO=PATM/(R*T) 

DTPTS=RWS*COLS 
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*  OPEN  A  NEW  FILE  TO  STORE  THE  REDUCED  DATA 

OPEN(2,FILE= 'RESULT.DAT' ,STATUS= 'OLD' ) 
WR1TE(2,222)  DTPTS 
222    FORMAT  (15) 

*  OPEN  A  SEQUENTIAL  OF  DATA  FILE 

DO  20  1=1, COLS 

NAME(7:8 )=A( I ) 
FNAME=NAME 
OPEN( 1,FILE=FNAME) 
READ( 1 , 100,END=20 )ST 
100     FORMAT (A6 5) 

15       READ(1,  1000,END=30)NO,X,Y,V1,V2,V3,V4,V5.BETA 
1000     FORMAT ( I2.F7.2.F6. 2 , 5F9 . 3 . F8 . 2 ) 

*  CONVERT  THE  VOLTAGE  TO  PRESSURE  IN  LBF/FT**2 

P1=DELTAP(V1,SL0PE, INTR ) ♦ 2 . 04  7 5+PATM 
P2  =  DELTAP ( V2 . SLOPE , I NTR ) *  2 . 04  7  5  +  PATM 
P3=DELTAP(V3, SLOPE, I NTR ) *2 . 047 5+PATM 
P4=DELTAP(V4, SLOPE, 1NTR ) *2 . 047 5+PATM 
P5=DELTAP(V5, SLOPE, I NTR ) *2 . 0475+PATM 

*  CALCULATE  THE  PITCH  ANGLE  IN  DEGREES 

F=(P4-P5 )/(Pl-P2 ) 
ALPHA=FPITCH(P) 

*  CALCULATE  THE  VELOCITY  IN  FT/SEC 

YSLOP=FYSLOP( ALPHA) 

VELM-SQRT( ( 2*YSL0P*(P1-P2) )/(RO*K) ) 

VEL=VELM+( 1+E) 

*  CALCULATE  THE  LOCAL  DYNAMIC  PRESSURE 

QM1=QM1FAC*2 . 04  7  5/K 
QM=RO+VEL**2/2 . 
Q1=QM1*( 1+2+E) 
Q=QM+( 1+2+E) 

*  CALCULATE  THE  YAW  ANGLE  IN  DEGREES 

YAW=FYAW(BETA+  5.0) 

*  CALCULATE  THE  VELOCITY  COMPONENTS 

BETAR=YAW* .017453 

ALPHAR= (ALPHA- 17 .942)*. 017453 

VELY=VEL*SIN(ALPHAR) 

VELX=VEL*COS ( ALPHAR ) +SIN( BETAR ) 

*  CALCULATE  THE  TOTAL  PRESSURE  IN  LBF/IN**2 

PTC=FPT( ALPHA) 
PT1=P1-0*PTC 
PT=PT1/144 . 
CPT=(PT1-PATM-Q1)/01 

*  CALCULATE  THE  STATIC  PRESSURE  IN  LBF/IN**2 

PSl=PTl-0 

PS=PS1/144 

CPS=(PS1-PATM)/Q) 

*  WRITE  ZEPO  VALUES  IF  VELOCITY  IS  TOO  HIGH  (BAD) 

IF(VEL.GT. 200.0)  THEN 
Z  -  0  .  0 

WRITE  (2.2000)-X+XOFF,Y+YOFF,Z.Z,Z,Z,Z,Z.Z,Z,Z 
ELSE 

WRITE ( 2 , 2000 ) - X+XOFF , Y+YOFF , VEL . VELX . VELY , YAW , 
C      ALPHA-17. 942, PT.CPT.PS, CPS 
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2000       FORMAT ( 11F10. 3 ) 
ENDIF 
GO  TO  15 
30     CLOSE(l) 
20     CONTINUE 
CLOSE(2) 
STOP 

END 

************************************************************* 

*  THIS  FUNCTION  CONVERTS  THE  VOLTAGE  TO  PHYSICAL  PRESSURE 

FUNCTION  DELTAP (X, SLOPE, INTR) 
REAL  INTR 
DELTAP=X*SLOPE+INTR 

END 

************************************************************* 

*  THIS  FUNCTION   CALCULATES  THE  PITCH  ANGLE 

FUNCTION  FPITCH(X  ) 

FPlTCH=3.7  59  +  5  3.  7.568*X-1 .  3085*X**2-1 .  6  583*X**3 

*         -0.8061*X**4+16.5115*X**5 

END 

************************************************************ 

*  THIS  FUNCTION  CALCULATES  THE  VELOCITY  PRESSURE  COEFFICENT 

FUNCTION  FYSLOP(X) 
IF(X.LT. -10)THEN 

FYSLOP  =  0.981-0.010  2*X-3.000E-«»  +  X+*2-2.500E-6*X**3 
ELSE  IF( (X.GB. -10) .AND. (X.LE. 10) )THEN 

FYSLOP= 0.98-0. 006 +X+ 2  .  OOOE-4*X**2 
ELSE 

FYSLOP=0.  9801-0.  003  5*X-1.  143E-(»*X**2  +  5  .  8  33E- 6*X**  3 
END  IF 

END 

************************************************************* 

*  THIS  FUNCTION  CALCULATES  THE  YAW  ANGLE 

FUNCTION  FYAW(X)  ' 

IF( (X.GE.O) .AND. (X.LE. 180) )  THEN 

FYAW=-X 
ELSE 

FYAW=360-X 

END  IF 

END 
************************************************************* 

*  THIS  FUNCTION  CALCULATES  THE  TOTAL  PRESSURE  COEFFICIENT 

FUNCTION  FPT(X) 
IFIX.LE.-30)  THEN 

FPT=-0.01 
ELSE  IF( (X.GT. -30) .AND. (X.LT. -20) )  THEN 

FPT=0.02+1.00E-3*X 
ELSE  IF( (X.GE. -20) .AND. (X.LE.30) )  THEN 

FPT  =  0 
ELSE 

FPT  =  0.03-1 .  r>0E-3*X 
END  IF 
END 
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APPENDIX  D.  RESULTOO.DAT 
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1 

.000 

.000 

122.209 

-7.000 

18. 214 

15.002 

.017     1 

4.875 

-.060 

2 

.000 

.500 

122. 155 

-7.000 

18.388 

15.002 

.015     1 

4.875 

-.061 

3 

.000 

1.000 

122.267 

-7.000 

18.254 

15.002 

.017     1 

4.875 

-  .062 

4 

.000 

1.500 

122.199 

-7.000 

18.297 

15.002 

.016     1 

4.875 

-.061 

5 

.000 

2.000 

122.310 

-7.000 

18.433 

15.002 

.018     1 

4.875 

-.061 

6 

.000 

2.500 

122.155 

-7.000 

18. 388 

15.002 

.016     1 

4.875 

-.061 

7 

.000 

3.000 

122. 293 

-7.000 

18. 302 

15.002 

.018     1 

4.875 

-.061 

8 

.000 

3.  500 

122.359 

-7.000 

18.287 

15.002 

.018     1 

4.875 

-  .062 

9 

.000 

4  .000 

122.287 

-7.000 

18.357 

15.002 

.017     1 

4.875 

-.062 

10 

.000 

4.  500 

125.471 

-5.000 

17. 651 

15.002 

.017     1 

4.  868 

-  .  118 

11 

.000 

5.000 

125. 733 

-5.000 

17.773 

15.002 

.018     1 

4.868 

-.122 

12 

.000 

5.  500 

125.881 

-5.000 

17.614 

15.002 

.019     1 

4.867 

-.  124 

1 

.500 

.000 

125.918 

-5.000 

17.784 

15.002 

.017     1 

4.867 

-.127 

2 

.500 

.500 

125.733 

-5.000 

17.773 

15.002 

.018     1 

4.868 

-.122 

' 

3 

.500 

1.000 

126. 131 

-5.000 

17.814 

15.002 

.015     1 

4.866 

-.132 

4 

.500 

1.500 

125.773 

-5.000 

17.714 

15.002 

.016     1 

4.867 

-.125 

5 

.500 

2.000 

126.015 

-5.000 

17.764 

15.002 

.015     1 

4.867 

-.131 

6 

.500 

2.500 

125.812 

-5.000 

17.655 

15.002 

.013     1 

4.867 

-.128 

7 

.500 

3.000 

125. 624 

-5.000 

17.669 

15.001 

.012     1 

4.867 

-.126 

8 

.500 

3.500 

125.731 

-5.000 

17. 570 

15.002 

.015     1 

4.867 

-.126 

9 

.500 

4.000 

125.777 

-5.000 

17.688 

15.002 

.014     1 

4.867 

-.127 

10 

.500 

4.500 

125.944 

-5.000 

17.601 

15.002 

.013     1 

4.867 

-.131 

11 

.500 

5.000 

126.524 

-5.000 

17.454 

15.002 

.015     1 

4.865 

-.140 

12 

.500 

5.500 

125.918 

-5.000 

17.784 

15.002 

.015     1 

4.  867 

-.129 

1 

1.000 

.000 

126 

.  500 

-5 

.000 

17 

.408 

15 

.001 

4 

.  865 

-.  142 

2 

1.000 

.500 

126 

.  304 

-5 

.000 

17 

.  474 

15 

.002 

k 

.866 

-.136 

3 

1.000 

1 

.000 

126 

.189 

-5 

.000 

17 

.625 

15 

.002 

1* 

.866 

-.134 

It 

1.000 

1 

.500 

125 

.976 

-5 

.000 

17 

.822 

15 

.001 

4 

.  866 

-.  13  it 

5 

1.000 

2 

.000 

125 

.656 

-5 

.000 

17 

.662 

15 

.001 

4 

.867 

-.  126 

6 

1.000 

2 

.500 

125 

.659 

-5 

.000 

17 

.636 

15 

.001 

4 

.867 

-.127 

7 

1  000 

3 

.000 

125 

.571 

-5 

.000 

17 

.604 

15 

.001 

4 

.867 

-  126 

8 

1.000 

3 

.500 

125 

.821 

-5 

.000 

17 

.602 

15 

.001 

k 

.867 

-.131 

9 

1.000 

it 

.000 

125 

.881 

-5 

.000 

17 

.614 

15 

.001 

4 

.867 

-.130 

10 

1.000 

it 

.500 

125 

.617 

-5 

.000 

17 

.722 

15 

.001 

It 

.867 

-.129 

11 

1.000 

5 

.000 

125 

.492 

-5 

.000 

17 

.723 

15 

.001 

4 

.867 

-.125 

12 

1.000 

5 

.  500 

126 

.317 

-5 

.000 

17 

.396 

15 

.002 

4 

.866 

-.137 

1 

1.500 

.000 

126 

.253 

-5 

.000 

17 

.611 

15 

.001 

It 

.  866 

-.  138 

2 

1.  500 

.500 

125 

.896 

-5 

.000 

17 

.  509 

15 

.001 

It, 

,867 

-.131 

3 

1.500 

1, 

.000 

125 

.926 

-5 

.000 

17 

.732 

15 

.001 

k . 

866 

-.133 

4 

1.500 

1, 

.500 

125, 

.860 

-5 

.000 

17 

.  543 

15 

.001 

1* . 

866 

-  .  132 

5 

1.500 

2, 

.000 

125, 

,405 

-5, 

.000 

17 

.895 

15 

.001 

It. 

867 

-.  126 

6 

1.500 

2. 

,500 

125. 

,941 

-5 

,000 

17, 

.627 

15, 

.001 

It. 

866 

-.135 

I 

7 

1.  500 

3. 

000 

126. 

016 

-5. 

,000 

17. 

,536 

15, 

,001 

It. 

866 

-.139 

8 

1.500 

3. 

500 

126. 

063 

-5. 

,000 

17. 

424 

15. 

,001 

It. 

866 

-.  136 

9 

1.500 

4. 

000 

125. 

824 

-5. 

000 

17. 

576 

15. 

,001 

u  . 

866 

-.  135 

10 

1.500 

it . 

500 

125. 

881 

-5. 

000 

17. 

614 

15. 

001 

4. 

866 

-.  133 

11 

1.500 

5. 

000 

125. 

853 

-5. 

000 

17. 

594 

15. 

001 

4. 

866 

-.134 

12 

1.500 

5. 

500 

125. 

784 

-5. 

000 

17. 

635 

15. 

001 

4 . 

866 

-.  133 

1 

2.000 

000 

125. 

709 

-5. 

000 

17. 

727 

15. 

001 

It . 

866 

-.133 

2 

2.000 

500 

125. 

242 

-5. 

000 

17. 

725 

15. 

001 

4 . 

867 

-.  125 

3 

2.000 

1. 

000 

125. 

934 

-5. 

000 

17. 

679 

15. 

001 

It. 

866 

-.  138 

It 

2.000 

1. 

500 

125. 

925 

-5. 

000 

17. 

529 

15. 

001 

4. 

866 

-.137 

5 

2.000 

2. 

000 

125. 

833 

-5. 

000 

17. 

726 

15. 

001 

1* . 

866 

-.136 

6 

2.000 

2. 

500 

125. 

844 

-5. 

000 

17. 

876 

.15. 

001 

It. 

866 

-.136 

.012 

1 

.015 

1 

.015 

1 

.011 

1 

.012 

1 

,012 

1 

.011 

1 

.011 

1 

.013 

1 

.009 

1 

,010 

1 

.014 

1 

,011 

1 

,012 

1 

,011 

1 

Oil 

1 

008 

1 

009 

1 

007 

1 

010 

1 

007 

1 

010 

1 

008 

1 

008 

1 

007 

1 

007 

1 

006 

1 

006 

1 

006 

1 

006 

1 
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7 

2.000 

3.000 

125.890 

-5.000 

17 

.765 

15.001 

.008 

1 

4,-866 

-.  135 

8 

2  .000 

3.  500 

126.253 

-5.000 

17 

.611 

15.001 

.009 

1 

4.865 

-  .  141 

9 

2  .000 

4.000 

125.720 

-5.000 

17 

.649 

15.001 

.008 

1 

4.866 

-.132 

10 

2.000 

4.500 

125.289 

-5.000 

17 

.613 

15.001 

.006 

1 

4.867 

-.127 

11 

2.000 

5.000 

125.982 

-5.000 

17 

.770 

15.001 

.009 

1 

4.866 

-.136 

12 

2.000 

5.500 

125.549 

-5.000 

17 

.762 

15.001 

.007 

1 

4.867 

-.  130 

1 

2.500 

.000 

125.971 

-5.000 

17 

.848 

15.001 

.005 

1 

4.866 

-.139 

2 

2.500 

.500 

125.873 

-5.000 

17 

.667 

15.001 

.006 

1 

4.866 

-.13  7 

3 

'2.500 

1.000 

125.760 

-5.000 

17 

.589 

15.001 

.005 

1 

4.866 

-  .  136 

4 

2.500 

1.500 

125.886 

-5.000 

17 

.791 

15.001 

.006 

1 

4.866 

-.137 

5 

2.  500 

2.000 

125.235 

-4.000 

17 

.778 

15.000 

.005 

1 

4.867 

-.127 

6 

2.500 

2.  500 

125.696 

-4.000 

17 

.603 

15.000 

.005 

1 

4.866 

-.135 

7 

2.500 

3.000 

125.117 

-4.000 

17 

.726 

15.000 

.005 

1 

4.867 

-.  124 

8 

2.500 

3.500 

125.394 

-4.000 

17 

,539 

15.000 

.004 

1 

4.867 

-  .  130 

9 

2.500 

4.000 

125. 169 

-4.000 

17 

.793 

15.001 

.006 

1 

4.867 

-.124 

10 

2.500 

4.500 

125.871 

-4.000 
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1  .000  .000    150.248    -16.000     15.148     14.700      -.041 
4.516  -.604 

2  .000  .250    155.911    -15.000     14.475     14.701      -.039 
4.502  -.722 

3  .000  .500    160.315    -15.000     13.866     14.701      -.037 
4.491  -  .817 

4  .000  .750    164.509    -14.000     13.393     14.701      -.037 
4. 480  -.911 

5  .000  1.000    169.716    -14.000     12.697     14.701      -.039 
4.465  -1.033 

6  .000  1.250    175.393    -13.000     10.631     14.700      -.047 
4. 448  -1. 177 

7  .000  1.500    180.980     -8.000      2.980     14.688      -.147 
4.420  -1.414 

8  .000  1.750    185.168      9.000     -5.242     14.637      -.581 
4.357  -1.955 

9  .000  2.000    179.544     39.000    -23.961     14.556     -1.262 
4.293  -2.493 

10  .000  2.250    181.551     73.000    -13.886     14.587      -.999 
4. 318  -2. 281 

11  .000  2.500    199.011     84.000     -1.061     14.673      -.274 
4. 349  -2.016 

12  .000  2.750    192.122     81.000      -.585     14.649      -.480 
4.347  -2.035 

13  .000  3.000    190.920     61.000    -17.689     14.604      -.856 
4 . 307  -2. 379 

14  .000  3.250    186.191     24.000     -3.723     14.533     -1.456 
4 . 250  -2 . 857 

15  .000  3.500    190.574      4.000      3.291     14.621      -.713 
4. 325  -2. 227                                                                j 

16  .000  3.750    188.810     -8.000      4.141     14.(679      -.219 
4 . 388  -1.687 

17  .000  4.000    178.573    -11.000      3.863     14.693      -.102 
4. 433  -1. 310 

18  .000  4.250    171.008    -13.000      4.608     14.696      -.077 
4.4  57  -1.101 

19  .000  4.500    162.822    -33.000      6.181     14.696      -.075 

4.480  - . 911 

20  .000  4.750    155.653    -13.000      8.050     14.697      -.074 
4.499  -.752 

21  .000  5.000    151.143    -12.000      9.417     14.696      -.075 
4. 510  - .656 

22  .000  5.250    150.027    -12.000     10.473     14.697      -.073 
4. 513  - .631 

23  .000  5.500    147.773    -12.000     11.461     14.697      -.073 
4.518  -.585 

1  .250  .000    147.058    -12.000     16.164     14.692      -.112 
4.515  -.609 

2  .250  .250    152.564    -16.000     15.207     14.696      -.077 
4.506  -.688 

3  .250  .500    157.323    -16.000     14.221     14.696      -.076 
4 . 494  -  .789 

4  .250  .750    159.896    -17.000     14.424     14.690      -.131 

4.481  -.901 

5  .250  1.000    169.024    -16.000     12.905    .14.696      -.078 
4.463  -1.056 
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4.613  .018 
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4.609  -.034 

3  .000 

4.604  -.092 
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4.598  -.152 

5  .000 

4.596  -.180 
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4.586  -.288 
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4.566  -.522 
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4.526  -.970 
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4.483  -1. 465 
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4.501  -1.258 

11  .000 
4. 531  - .913 
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4.519  -1.055 
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14  .000 
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10 

2.500 

2.250 

84.224 

4  .000 

22. 530 

14.692 

-.085     1 

4.635 

.  262 

11 

2.500 

2.500 

84.007 

5.000 

21.010 

14.692 

-.086     1 

4.635 

.264 

12 
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2.750 

84 .061 

6.000 

19. 271 

14.693 

-.079     1 

d.636 

.  271 

13 

2.500 

3.000 

84.715 

5.000 

16.808 

14.692 

-.088     1 

4.634 

.251 

14 

2.500 

3.250 

84.589 

4.000 

15.321 

14.692 

-.095     1 

k  .63<* 

.246 

15 

2.500 

3.500 

86.336 

2.000 

13.493 

14.692 

-.094     1 

4.631 

.219 

16 

2.500 

3.750 

87. 139 

1.000 

12.460 

14.691 

-.096     1 

4.630 

.205 

17 

2.500 

4.000 

88.613 

.000 

12.040 

14.692 

-.095     1 

4.628 

.  182 

18 

2.500 

4.250 

8  9.726 

-2.000 

11.616 

14.692 

-.094     1 

4  .  626 

.164 

19 

2.  500 

4  .  500 

90.245 

-3.000 

11.725 

14.692 

-.090     1 

4.626 

.  160 

20 

2.  500 

4  .750 

90.971 

-4.000 

12.015 

14.692 

-.093     1 

4.  625 

.  145 

21 

2.  500 

5.000 

91.  511 

-5.000 

12.219 

14.692 

-.094     1 

4.  624 

.  135 

- 

22 

2.500 

5.250 

91. 362 

-6.000 

12.652 

14.692 

-.095     1 

4.  624 

.  137 

23 

2.500 

5.  500 

.  91.328 

-6.000 

12.983 

14.691 

-.096     1 

4  .  624 

.136 

1 

2.  750 

.000 

91 . 392 

-7.000 

19. 116 

14.693 

-.077     1 

4.  625 

.  155 

2 

2.750 

.250 

91.334 

-6.000 

19.617 

14.693 

-.075     1 

4.626 

.  158 

i 

3 

2.750 

.500 

91 .078 

-6.000 

19.948 

14  .69,3 

-.073     1 

4.626 

.  163 

4 

2.  750 

.750 

89.948 

-4.000 

20.765 

14.694 

-.070     1 

4.  628 

.  185 

5 

2.750 

1.000 

90. 279 

-3.000 

20.865 

14.694 

-.069     1 

4.628 

.  181 

6 

2.750 

1.250 

89. 358 

-2.000 

21. 300 

14.694 

-.071     1 

4.  629 

.  194 

7 

2.750 

1.500 

88. 115 

-2.000 

21.736 

14.693 

-.077     1 

4.  630 

.209 

8 

2.750 

1.750 

87.428 

-1.000 

21.959 

14.693 

-.077     1 

4.  631 

.219 

9 

2.750 

2.000 

86. 165 

1.000 

21.878 

14.693 

-.078     1 

4  .  633 

.239 

10 

2.750 

2.250 

85.369 

4.000 

21.404 

14.693 

-.081     1 

4  .634 

.248 

11 

2.750 

2.  500 

85. 128 

5.000 

20.256 

14.693 

-.080     1 

4.634 

.253 

12 

2.750 

2.750 

85.040 

5.000 

18.772 

14.693 

-.079     1 

4.634 

.255 

13 

2.750 

3.000 

85. 167 

3.000 

17.104 

14.692 

-.090     1 

4.633 

.242 

14 

2.750 

3.250 

85.256 

2.000 

15.700 

14.692 

-.095     1 

4.633 

.236 

15 

2.750 

3.500 

86.455 

2.000 

14.368 

14.691 

-.096     1 

4.631 

.  216 
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2.   /bU 

J  .  7bU 

u;  .240 

.UUU 

1J . 3bU 

14 .oyj 

-.oyi    i 

4 
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.208 

17 

2.  7  50 

4.000 

88.755 

-1.000 

12.647 

14.692 

-.095     1 

U 

628 

.  180 

18 

2.750 
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88.990 

-2  .000 

12.  573 

14.691 

-.096     1 

4 

627 

.  175 

19 

2.  750 
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89.990 

-3.000 

12.328 

14.692 

-.094     1 

4 

626 

.  161 

20 

2.750 

4.750 

90.463 

-4.000 

12.409 

14.692 

-.094     1 

4 

625 

.  152 

21 

2.750 

5.000 

90.749 

-5.000 

12.439 

14.692 

-.095     1 

4 

625 

.  147 

22 

2.750 

5.250 

91.233 

-5.000 

12.568 

14.691 

-.098     1 

4 

624 

.136 

23 

2.750 

5.  500 

90.735 

-5.000 

13.073 

14.691 

-.101     1 

4. 

624 

.  142 

1 

3.000 

.000 

91.045 

-6.000 

19.015 

14.693 

-.078     1 

4 

626 

.  159 

2 

3.000 

.250 

90.844 

-5.000 

19.513 

14.694  ' 

-.071     1 

4 

627 

.  169 

3 

3.000 

.  500 

90.071 

-5.000 

20. 110 

14 .693 

-.077     1 

It  . 

627 

.  177 

4 

3  .000 

.750 

89.  285 

-4.000 

20. 535 

14.693 

-.075     1 

It 

629 

.  191 

5 

3  .000 

1  .000 

89.378 

-3.000 

20.613 

14.694 

-.071     1 

k 

629 

.  194 

6 

3.000 

1.250 

88.  641 

-2.000 

20.975 

14.693 

-  - .075     1 

It 

630 

.  201 

7 

3.000 

1  .  500 

87  .780 

-1.000 

21.311 

14.693 

-.076     1 

it 

631 

.  215 

8 

3.000 

1.750 

86.667 

.000 

21. 519 

14.693 

-.081     1 

4 

632 

.  227 

9 

3.000 

2  .000 

86. 102 

1  .000 

21.138 

14.693 

- .079     1 

U 

633 

.  238 

10 

3  .000 

2.  250 

85.416 

3  .000 

20. 609 

14 . 693 

-.079     1 

4 

634 

.  250 

11 

3.000 

2  .  500 

85.428 

2.000 

19.804 

14.693 

-.079     1 

It 

634 

.  249 

( 

12 

3  .000 

2.750 

85. 133 

3.000 

18.556 

14.693 

-.083     1 

k 

634 

.  250 

13 

3.000 

3.000 

85.059 

2.000 

17.413 

14.692 

-.092     1 

It 

633 

.24  2 

14 

3  .000 

3.  250 

85.748 

1.000 

16.023 

14.692 

-.094     1 

4 

632 

.  229 

15 

3.000 

3.  500 

85.929 

1.000 

15.21-9 

14. 691 

-.099     1 

u 

631 

.221 

16 

3.000 

3.750 

87 . 250 

.000 

14.095 

14.692 

-.094     1 

It 

.630 

.  205 

17 

3.000 

4.000 

88 .061 

.000 

13.  537 

14.692 

-.093     1 

It 

.629 

.  193 

18 

3.000 

4  .  250 

88.435 

-2 .000 

13 . 227 

14.691 

-.097     1 

u 

.628 

.  183 

19 

3.000 

4.500 

89.064 

-2.000 

13.298 

14.691 

-.097     1 

It 

.627 

.  173 

20 

3.000 

4.750 

88.967 

-3.000 

13.290 

14.691 

- . 101     1 

4 

.627 

.  171 

21 

3.000 

5.000 

89.440 

-3.000 

13.249 

14.691 

-.099     1 

4 

.626 

.  164 

22 

3.000 

5.250 

89.756 

-4.000 

13.409 

14.691 

-.101     1 

4 

.626 

.  157 

23 

3.000 

5.500 

89.789 

-4.000 

13.653 

14.691 

-.101     1 

It 

.626 

.  157 
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APPENDIX  G.  RESULTOC.DAT 


299 

1  .000 
4.585  -.225 

2  .000 
4 . 581  - .262 

3  .000 
4  . 576  - .301 

4  .000 
4.  576  -  .  304 

5  .000 
4.  570  -.351 

6  .000 
4.554  -.<*88 

7  .000 
4.529  -.699 

8  .000 
4.501  -.942 

9  .000 

4 .476  -1. 148 

10  .000 

4 .477  -1. 145 

11  .000 
4.481  -1.112 

12  .000 
4.445  -1 .418 

13  .000 
4.411  -1.699 

14  .000 
4.423  -1.597 

15  .000 
4.468  -1.220 

16  .000 
4.491  -1.026 

17  .000 
4.525  -.736 

18  .000 
4.545  -.564 

19  .000 
4.558  -.452 

20  .000 

4. 561  -.427 

21  .000 
4. 567  - .384 

22  .000 
4.572  -.336 

23  .000 
4. 574  - . 321 

1  .250 
4.585  -.228 

2  .250 
4.576  -.305 

3  .250 
4.570  -.355 

4  .250 

4.562  -.424 

5  .250 
4.557  -.462 


000 

135. 

701 

-  12 

000 

14  . 

346 

14. 

737 

250 

140. 

986 

-13 

000 

12. 

7  34 

14. 

745 

500 

143 

144 

-14 

000 

12. 

003 

14. 

745 

750 

143 

325 

-15 

000 

10. 

993 

14. 

745 

1 

000 

145 

356 

-17 

000 

9 

953 

14 

745 

1 

250 

143 

712 

-15 

000 

4 

871 

14 

725 

1 

500 

149 

850 

4 

000 

-11 

861 

14 

715 

1 

750 

153 

629 

15 

000 

-19 

567 

14 

695 

2 

000 

151 

278 

44 

000 

-29 

743 

14 

665 

2 

2  50 

162 

101 

70 

000 

-12 

751 

14 

694 

2 

500 

167 

065 

77 

000 

2 

944 

14 

711 

2 

750 

175 

726 

73 

000 

9 

831 

14 

699 

3 

000 

172 

112 

55 

000 

16 

570 

14 

656 

3 

250 

168 

349 

26 

000 

21 

499 

14 

657 

3 

500 

168 

201 

8 

000 

21 

068 

14 

701 

3 

750 

149 

545 

-6 

000 

8 

895 

14 

675 

4 

000 

157 

828 

-17 

000 

8 

079 

14 

.731 

4 

.250 

154 

562 

-16 

.000 

12 

.284 

14 

.  742 

4 

.  500 

149 

529 

-14 

000 

12 

413 

14 

.  743 

4 

.750 

148 

561 

-14 

000 

12 

446 

14 

.  744 

5 

.000 

146 

364 

-12 

.000 

12 

679 

14 

.  743 

5 

.250 

143 

843 

11 

.000 

12 

963 

14 

.743 

5 

.500 

143 

.081 

-11 

.000 

13 

.261 

14 

743 

.000 

132 

.867 

-12 

.000 

15 

.605 

14 

731 

.250 

142 

.775 

-13 

.000 

13 

.  341 

14 

.744 

.500 

145 

.758 

-15 

.000 

12 

.454 

14 

.745 

.750 

148 

.977 

-16 

.000 

11 

.239 

14 

745 

1 

.000 

150 

.799 

-18 

.000 

9 

.884 

14 

.745 

.063     1 

.  128     1 

.132     1 

.  132     1 

.  127     1 

-.043     1 

-.128     1 

-.291     1 

- . 547     1 

-.307     1 

-.160     1 

- .258     1 

-  .627     1 

-.615     1 

-.241     1 

-.462     1 

.007     1 

.  107     1 

.  112     1 

.  116     1 

.  114     1 

.  Ill     1 

.  Ill     1 

.007     1 

.121     1 

.131     1 

.  128     1 

.128     1 
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6 

.250 

1  . 

250 

154. 

569 

-16.000 

6.363 

14.734 

.038      1 

It 

537 

-  .633 

7 

.250 

1  . 

500 

162. 

4  36 

-2.000 

-5.  137 

14  .  7  20 

-.084      1 

It 

502 

-  .930 

8 

.250 

1. 

750 

162. 

148 

12.000 

-16.843 

14.673 

-  .483     1 

u 

456 

-1.  321 

9 

.  2  50 

2. 

000 

162. 

590 

42.000 

-27  .952 

14. 622 

-.911     1 

4 

404 

-1.760 

10 

.250 

2. 

250 

164. 

872 

74.000 

-13  .  870 

14.667 

- .533     1 

1* 

kit  3 

-1.435 

11 

.250 

2. 

500 

173. 

461 

82.000 

.397 

14.706 

-.201     1 

4 

458 

-1.  306 

12 

.250 

2. 

750 
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977 
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2.980 

14.688 

-  .355     1 

It 

432 

-1.521 

13 

.250 

3. 

000 

149. 

982 

60.000 

2.154 

14.594 

-1.149     1 

4 

t*09 

-1.722 

14 

.250 

3. 

250 

168. 

175 

19.000 

7.057 

14.611  • 

-1.011    1 

It. 

377 

-1.989 

15 

.250 

3 

500 

175. 

525 

3.000 

10.704 

14.681 

-.410     1 

4 

427 

-1.  565 

16 

.250 

3 

750 

167. 

877 

-10.000 

6.064 

14.709 

-.172     1 

It 

477 

-1.  143 

17 

.  250 

4 

000 

159. 

261 

-17.000 

9.054 

14.732 

.023     1 

It 

523 

-  .751 

18 

.250 

4 

250 

154 

163 

-16.000 

11. 100 

14.742 

.106     1 

4 

546 

-  .556 

19 

.  250 

4 

500 

14  8 

679 

-15.000 

11 .766 

14.744 

.  118     1 

It 
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-  .428 

20 

.250 

4 
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"14  6 

031 

-14.000 

12.087 

14.744 

.  118     1 

It 

568 

-  .  373 

21 

.  250 

5 
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145 

914 

-12.000 

12. 199 

14.744 

.  121     1 

It 

568 

-  .  368 

22 

.  250 

5 

250 

143 

882 
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12. 613 

14.744 

.118     1 

It 
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-  .  330 

23 

.250 

5 
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481 

-10.000 

12. 946 

14. 744 

.  118     1 

It 

576 
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1 
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14.892 
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.103     1 

It 

582 
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.250 
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098 
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13.863 
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It 
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3 
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454 
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13. 243 

14.  744 

.119     1 

k 
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U 
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14.744 

.  120     1 

It 
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.454 

-18.000 

10.963 
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.117     1 

It 

.  542 

-  .  595 
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.  500 
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.  250 

160 

.952 
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7.847 

14.730 

.072     1 
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.  524 

-  .740 

7 

.  500 
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.533 

-9.000 

1.  301 

14.718 
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It 

.484 

-1.085 
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.500 
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.098 

5.000 

-3.847 

14.667 
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It 

.447 

-1.  395 
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.480 

3  5.000 

-26.721 

14.566 
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It 

.396 

-1.826 
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.707 
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-4. 380 

14.669 
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It 

.426 

-1. 575 
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.794 
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14.723 
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.447 
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1? 

.500 

2 

.750 

103 

.419 

78.000 

-2.722 
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10 
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175 

792 
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7.052 

14.680 
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-1.  586 
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5.  303 

14.719 
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.459 

-1.294 

12 
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180 
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-1.937 

14.719 
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-1.  371 
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14 
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.750 
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-1.039 

17 

.750 

4 

.000 

160 

.593 

-15.000 

4.599 

14.741 

.092     1 

It 

.528 

-.711 

18 

.750 

4 

.250 

151 

.566 

-14 .000 

7.923 

14.742 

.101     1 

It 

.  552 

-  .506 

19 

.750 

4 

.500 

144 

.754 

-13.000 

9.982 

14.742 

.102     1 

It 

.  569 

-.  363 

112 


20 

.750 

4  . 

750 

142. 

311 

-  12 .000 

11. 

115 

14  . 

74  2 

.  1U4     1 

4 

575 

-  .  312 

21 

.750 

5. 

000 

14  2. 

236 

-11 .000 

1  1  . 

374 

14  . 

742 

.  102     1 

It 

575 

-  .  313 

22 

.750 

5. 

250 

141. 

997 

-10.000 

11  . 

797 

14. 

742 

.  104     1 

It 

576 

-  .306 

23 

.750 

5. 

500 

140. 

999 

-9.000 

12. 

260 

14. 

742 

.  105     1 

I* 

578 

-  .  286 

1 

1.000 

000 

139. 

824 

-12.000 

15. 

877 

14. 

741 

.094      1 

It 

580 

-.274 

2 

1  .000 

250 

141. 

010 

-13.000 

15. 

974 

14. 

742 

.  102     1 

4 

578 

-  .  288 

3 

1.000 

500 

145. 

331 

-15.000 

15. 

286 

14. 

742 

.  104     1 

It 

568 

-.373 

4 

1.000 

750 

149. 

707 

-16.000 

15. 

314 

14. 

742 

.  103     1 

It. 

557 

-  .465 

5 

1.000 

1 

000 

155. 

589 

-17 .000 

14. 

970 

14. 

742 

.  103     1 

1* . 

542 

-.590 

6 

1.000 

1 

250 

162 

647 

-19.000 

14. 

998 

14. 

741 

.091     1 

ft. 

522 

-  .  759 

7 

1  .000 

1 

500 

167 

301 

-15.000 

14 

658 

14. 

725 

-.043     1 

It 

494 

-1  .000 

8 

1  .000 

1 

750 

167 

898 

-4.000 

16 

258 

14 

706 

-.203     1 

It 

473 

-1.  175 

9 

1  .000 

2 

000 

160 

472 

20.000 

19 

773 

14 

680 

-.425     1 

k 

467 

-1 . 226 

«■ 

10 

1  .000 

2 

250 

165 

125 

50.000 

14 

600 

14 

694 

-.302     1 

it 

469 

-1.  209 

11 

1.000 

2 

500 

.166 

154 

72.000 

8 

456 

14 

729 

-.009     1 

k 

501 

-  .  940 

12 

1  .000 

2 

750 

167 

262 

7  2.000 

1 

24  6 

14 

733 

.029     1 

k 

502 

-  .  927 

13 

1  .000 

3 

000 

164 

428 

51 .000 

-19 

543 

14 

696 

-.283     1 

k 

473 

-1.174 

i 

14 

1  .000 

3 

250 

164 

282 

20.000 

-25 

.  503 

14 

.693 

-  .  316     1 

it 

470 

-1 . 204 

15 

1  .000 

3 

.  500 

160 

203 

-3 .000 

-15 

.672 

14 

.715 

-  .  123     1 

it 

482 

-1.102 

16 

1.000 

3 

.750 

162 

036 

-12.000 

-5 

.  622 

14 

.  735 

.046     1 

it 

519 

-.791 

17 

1.000 

4 

.000 

154 

477 

-14.000 

2 

.555 

14 

.741 

.095     1 

u 

544 

-  .  574 

18 

1  .000 

4 

.  250 

14  0 

.954 

-13 .000 

6 

.621 

14 

.743 

.  114     1 

It 

560 

-.438 

• 

19 

1  .000 

4 

.500 

144 

.  130 

-12.000 

8 

.945 

14 

.  744 

.119     1 

4 

.572 

-  .334 

20 

1  .000 

4 

.  750 

141 

.491 

-11 .000 

10 

.368 

14 

.744 

.  119     1 

k 

.579 

-  .281 

21 

1  .000 

5 

.000 

141 

.625 

-  10.000 

11 

.289 

14 

.  744 

.118     1 

It 

.578 

-  .  285 

22 

1.000 

5 

.250 

139 

.  809 

-9.000 

11 

.967 

14 

.74  3 

.  115     1 

It 

.  582 

-  .  254 

23 

1.000 

5 

.500 

141 

.326 

-9.000 

12 

.  335 

14 

.  744 

.119     1 

It 

.579 

-  .  278 

1 

1  .250 

.000 

130 

.251 

-16.000 

18 

.744 

14 

.7  40 

.087     1 

It 

.600 

-  .099 

2 

1  .250 

.250 

132 

.031 

-17 .000 

18 

.427 

14 

.  743 

.  109     1 

It 

.599 

-  .  110 

3 

1  .250 

.  500 

133 

.449 

-10.000 

19 

.196 

14 

.743 

.  113     1 

It 

.596 

-.132 

» 

113 


4 

1  .250 

750 

135. 

252 

-20.000 

19. 

900 

14  . 

742 

.  107     1 

k  . 

592 

-  .  172 

5 

1.250 

1  . 

000 

137  . 

600 

-23 .000 

20. 

789 

14  . 

742 

.104     1 

u . 

586 

-  .  220 

6 

1.250 

1. 

250 

138. 

908 

-25.000 

22. 

934 

14. 

739 

.079     1 

4 . 

580 

-.  271 

7 

1.250 

1. 

500 

139. 

116 

-26.000 

26. 

961 

14. 

7  30 

.003     1 

4 . 

571 

-  .  350 

8 

1.250 

1. 

750 

131  . 

898 

-17.000 

31. 

812 

14. 

705 

-.210     1 

4 . 

562 

-.426 

9 

1.250 

2. 

000 

127. 

652 

2.000 

38. 

120 

14. 

692 

-.321     1 

4. 

557 

-.461 

10 

1.250 

2 

250 

117. 

633 

29.000 

31. 

219 

14. 

685 

-.378     1 

4. 

571 

-.346 

11 

1.250 

2. 

500 

129. 

667 

56.000 

20. 

458 

14. 

723 

-.059     1 

4. 

58d 

-.235 

12 

1.250 

2. 

750 

134. 

967 

58.000 

6. 

096 

14. 

738 

.074     1 

4. 

588 

-.200 

13 

1  .250 

3 

000 

133. 

208 

43.000 

-18. 

586 

14. 

715 

-.127     1 

4 . 

568 

-.368 

Id 

1.250 

3 

250 

130. 

398 

8.000 

-39. 

174 

14 

699 

-.264     1 

4 

558 

-.454 

15 

1.  250 

3 

500 

139 

705 

-8  .000 

-26 

762 

14 

728 

-.018     1 

4 

567 

-  .383 

16 

1  .250 

3 

750 

141 

158 

-12  .000 

-10 

015 

14 

737 

.061     1 

4 

573 

-.333 

- 

17 

1.250 

4 

000 

138 

039 

-14.000 

859 

14 

742 

.  106     1 

4 

585 

-.226 

18 

1.  250 

4 

250 

.135 

955 

-13 .000 

6 

206 

14 

743 

.  110     1 

4 

590 

-.  183 

19 

1.  250 

4 

500 

133 

573 

-13. 000 

8 

855 

14 

743 

.111     1 

'. 

596 

-  .  137 

20 

1  .  250 

4 

7  50 

133 

391 

-13 .000 

10 

640 

14 

.743 

.  113     1 

4 

596 

-.132 

i 

21 

1  .250 

5 

000 

133 

472 

-13.000 

11 

799 

14 

.743 

.  115     1 

4 

596 

-  .  131 

22 

1.250 

5 

2  50 

132 

834 

-12.000 

12 

.232 

14 

.74*3 

.113     1 

4 

597 

-.  122 

2  3 

1.250 

5 

.500 

132 

.974 

-12.000 

12 

.880 

14 

.743 

.  115     1 

4 

597 

-.122 

1 

1.  500 

.000 

134 

.633 

-10.000 

18 

.295 

14 

.735 

.042     1 

4 

585 

-  .226 

2 

1.500 

.250 

139 

.641 

-11.000 

17 

.47  5 

14 

.74  2 

.107     1 

<t 

582 

-  .  257 

3 

1.  500 

.  500 

142 

.828 

-12.000 

17 

.  458 

14 

.74  3 

.  112     1 

4 

575 

-  .315 

It 

1.500 

.7  50 

144 

.611 

-13.000 

18 

.004 

14 

.74  3 

.113     1 

4 

.571 

-  .350 

5 

1.500 

1 

.000 

148 

.444 

-14 .000 

10 

.  541 

14 

.744 

.116     1 

4 

.562 

-.425 

6 

1.500 

1 

.  250 

153 

.  197 

-15.000 

20 

.  159 

14 

.74  3 

.  113     1 

4 

.51*9 

-  .528 

7 

1  .500 

1 

.  500 

152 

.768 

-13.000 

24 

.  145 

14 

.742 

.  103     1 

4 

.549 

-  .529 

8 

1.  500 

1 

.750 

149 

.813 

-7 .000 

29 

.800 

14 

.735 

.045     1 

4 

.  550 

-.525 

9 

1.500 

2 

.000 

140 

.858 

10.000 

33 

.559 

14 

.717 

-.108     1 

4 

.553 

-  .496 

10 

1.500 

2 

.250 

13  2 

.801 

35.000 

30 

.724 

14 

.712 

-.154     1 

4 

.566 

-.  388 

» 

114 


11 

1.500 

2. 

500 

128. 

548 

53  .000 

25.001 

14.  733 

.024     1 

k 

596 

-.132 

12 

1.  500 

2. 

750 

125. 

862 

51 .000 

13. 554 

14.742 

.  100     1 

k 

611 

-.008 

13 

1.500 

3. 

000 

130. 

583 

41.000 

-5. 156 

14.740 

.089     1 

It 

600 

-  .  104 

14 

1.500 

3. 

250 

137. 

768 

21 .000 

-20. 669 

14.732 

.022     1 

u 

576 

-.306 

15 

1.500 

3. 

500 

141. 

283 

3.000 

17.513 

14.736 

.055     1 

It 

572 

-.341 

16 

1.500 

3. 

750 

143. 

024 

-5.000 

-7  .419 

14.741 

.099     1 

4 

573 

-.  332 

17 

1.  500 

4. 

000 

141. 

164 

-8.000 

.351 

14.743 

.  114     1 

It 

579 

-.280 

IB 

1.500 

4 

250 

139. 

051 

-9.000 

4.993 

14.743 

.114     1 

4. 

564 

-.239 

19 

1.500 

4. 

500 

138. 

066 

-8.000 

7  .706 

14.743  . 

.  114     1 

4 

586 

-  .220 

20 

1.500 

4 

750 

138. 

196 

-9.000 

9.  576 

14.743 

.116     1 

It . 

586 

-  .220 

21 

1.500 

5 

000 

139 

079 

-8.000 

11.097 

14.745 

.  132     1 

It . 

586 

-  .  221 

22 

1.  500 

5 

250 

139 

367 

-9.000 

11. 635 

14.  743 

.  115     1 

It . 

583 

-  .  243 

23 

1  .  500 

5 

500 

139 

201 

-8.000 

11  .  915 

14.743 

.114     1 

It . 

583 

-  .  242 

1 

1.750 

000 

136 

003 

-10.000 

18.036 

14.741 

.095     1 

4 . 

5P8 

-  .  199 

2 

1.750 

250 

•138 

038 

-10.000 

18.  130 

14. 743 

.111     1 

4 . 

586 

-  .  222 

3 

1.750 

500 

140 

745 

-11 .000 

18. 312 

14. 743 

.112     1 

4 . 

500 

-  .  274 

it 

1.750 

750 

141 

745 

-11.000 

18.  896 

14.743 

.  116     1 

4 

578 

-  .290 

i 

5 

1.750 

1 

000 

142 

073 

-11.000 

20. 125 

14.743 

.  113     1 

4 

577 

-  .299 

6 

1.750 

1 

250 

144 

335 

-12.000 

2  2.251 

14.743 

.  112     1 

It 

571 

-  .  345 

7 

1.750 

1 

500 

142 

761 

-10.000 

26.418 

14.743 

.  108     1 

u 

57'» 

-.318 

8 

1.750 

1 

750 

138 

.788 

-6.000 

33.601 

14.742 

.  102     1 

4 

583 

-  .245 

9 

1.750 

2 

.000 

135 

.251 

7.000 

3  9.858 

14. 737 

.063     1 

u 

586 

-.217 

10 

1.750 

2 

.250 

126 

.872 

28.000 

38. 345 

14.735 

.048     1 

It 

603 

-  .078 

11 

1.7  50 

2 

.500 

119 

.  810 

38.000 

28.628 

14.741 

.095     1 

fa 

.623 

.091 

12 

1.750 

2 

.750 

114 

.731 

30 .000 

19.354 

14.743 

.113     1 

It 

.631* 

.  192 

13 

1.750 

3 

.000 

118 

.630 

30.000 

5.7  04 

14.743 

.  Ill     1 

It 

.627 

.  127 

14 

1.750 

3 

.250 

124 

.230 

19.000 

-5. 190 

14.742 

.  102     1 

4 

.615 

.023 

15 

1.750 

3 

.500 

131 

.934 

6.000 

-7  .685 

14.742 

.  104     1 

It 

.598 

-.114 

16 

1.750 

3 

.750 

134 

.980 

-1  .000 

-3. 566 

14.743 

.  110     1 

It 

.592 

-  .  164 

17 

1.750 

4 

.000 

135 

.951 

-4 .000 

1.752 

14.744 

.117     1 

It 

.591 

-.175 

115 


18 

1.750 

4  .  250 

136. 167 

-6.000 

5.624 

14.743 

.  113     1 

4 . 

590 

-  .  184 

19 

1  .750 

4  .  500 

137 .488 

-7 .000 

7  .878 

14  .743 

.135     1 

4. 

587 

-  .207 

20 

1.750 

4.7  50 

136.808 

-7.000 

9.622 

14.743 

.113     1 

k . 

589 

-  .  196 

21 

1.750 

5.000 

138.312 

-8 .000 

10. 576 

14.743 

.  113     1 

It . 

585 

-.225 

22 

1.750 

5.250 

138.281 

-8.000 

11 .512 

14.743 

.115     1 

k . 

586 

-  .222 

23 

1.750 

5.500 

137.  170 

-8.000 

12.306 

14.743 

.113     1 

It . 

588 

-  .203 

1 

2  .000 

.000 

135. 553 

-9.000 

17 .879 

14.742 

.  105     1 

u . 

591 

-  .  180 

2 

2.000 

.250 

135.580 

-8.000 

18.774 

14.743 

.  113     1 

4. 

591 

-.  173 

3 

2.000 

.  500 

137.043 

-9.000 

19.083 

14.743 

.111     1 

4  . 

588 

-  .  203 

k 

2.000 

.750 

136. 293 

-8.000 

20.233 

14.743 

.113     1 

4. 

590 

-.  186 

5 

2.000 

1.000 

137. 823 

-9.000 

20.977 

14.74  3 

.109     1 

4 

586 

-  .220 

6 

2.000 

1.250 

135.653 

-8.000 

23.599 

14.743 

.  Ill     1 

4 

591 

-.  177 

7 

2.000 

1  .500 

134. 506 

-6.000 

27 .254 

14.743 

_  .111     1 

4 

594 

-.154 

3 

2.000 

1.750 

130.342 

.000 

32. 902 

14  .743 

.  Ill     1 

4 

603 

-  .077 

9 

2.000 

2  .000 

127 . 358 

8.000 

37.165 

14.743 

.112     1 

4 

609 

-  .022 

10 

2.000 

2.  250 

121.656 

19.000 

35.322 

14.743 

.  109     1 

It 

621 

.074 

11 

2  .000 

2  .  500 

117 .405 

27.000 

28 .040 

14.742 

.  106     1 

It 

629 

.  142 

12 

2.000 

2.750 

115.629 

28.000 

20.396 

14.743 

.  112     1 

It 

633 

.  176 

1 

13 

2.000 

3.000 

116.073 

14.000 

11.884 

14.743 

.  109     1 

It 

632 

.  167 

Id 

2.000 

3.  250 

119.315 

17.000 

3.472 

14.743 

.108     1 

It 

625 

.  112 

15 

2.000 

3.  500 

124 .807 

8.000 

.529 

14.743 

.108     1 

It 

614 

.038 

16 

2.000 

3.750 

131.573 

1  .000 

1.534 

14.743 

.  108     1 

It 

600 

-  .  102 

37 

2.000 

4  .000 

132.475 

-2.000 

4.  119 

14.743 

.  110     1 

4 

.  598 

-  .  117 

18 

2.000 

4.250 

134 .709 

-4.000 

6.642 

14.743 

.108     1 

u 

.593 

-  .  161 

19 

2.000 

4.500 

136.228 

-5.000 

8.425 

14.742 

.107     1 

4 

.  589 

-.  191 

20 

2.000 

4.750 

136.170 

-6.000 

9.841 

14.743 

.111     1 

4 

.590 

-  .  186 

21 

2.000 

5.000 

138.069 

-7.000 

10.800 

14.743 

.  109     1 

It 

.  585 

-  .224 

22 

2.000 

5.250 

137.440 

-7.000 

11.649 

14.743 

.113     1 

It 

.587 

-.208 

23 

2.000 

5.500 

136.740 

-7.000 

12.  133 

14.743 

.  110     1 

k 

.589 

-.  198 

1 

2.250 

.000 

135. 362 

-8.000 

18.392 

14.742 

.  107     1 

It 

.  591 

-.  175 

116 


2 

2.250 

250 

136. 

262 

-0 .000 

10.  740 

]4  .  742 

.  106     1 

k . 

589 

-  .  193 

3 

2.250 

500 

135. 

637 

-7 .000 

19. 325 

14  .  742 

.  10b     1 

u . 

590 

-.  182 

k 

2.250 

750 

135. 

]07 

-6.000 

20. 334 

14.743 

.  109     1 

4 . 

592 

-.167 

5 

2.250 

1. 

000 

134. 

338 

-6.000 

21.686 

14.742 

.  106     1 

u . 

593 

-  .156 

6 

2.250 

1. 

250 

132. 

534 

-6.000 

23.747 

14.742 

.106     1 

u  . 

597 

-.122 

7 

2.250 

1. 

500 

131. 

471 

-4.000 

26. 193 

14.742 

.106     1 

4 . 

600 

-  .  103 

8 

2.250 

1. 

750 

128. 

279 

-2.000 

29.454 

14.742 

.107     1 

4  . 

607 

-  .out* 

9 

2.250 

2. 

000 

124. 

259 

7  .000 

30.897 

14.742 

.  105     1 

4 . 

615 

.025 

10 

2.250 

2. 

250 

120. 

398 

13.000 

30.254 

14.742  . 

.105     1 

4 . 

623 

.091 

11 

2.250 

2. 

500 

117. 

111 

18.000 

25.819 

14.742 

.106     1 

4 . 

629 

.147 

12 

2.250 

2 

750 

115. 

876 

18.000 

21.013 

14.742 

.107     1 

4 . 

632 

.  168 

13 

2.250 

3 

000 

115 

611 

18.000 

15.443 

14  .  742 

.  105     1 

4 

632 

.  170 

14 

2.250 

3 

250 

120 

879 

12.000 

9.  187 

14.742 

.  107     1 

4 

622 

.085 

v 

15 

2.250 

3 

500 

123 

448 

8.000 

6.053 

14  .74  2 

.  104     1 

4 

616 

.038 

16 

2.250 

3 

750 

•12  7 

522 

3.000 

5.719 

14.742 

.106     1 

4 

608 

-  .031 

17 

2.250 

4 

000 

131 

820 

.000 

6.437 

14.743 

.108     1 

4 

599 

-  .  107 

18 

2.  250 

4 

250 

132 

895 

-3 .000 

7  .990 

14.742 

.  106     1 

4 

597 

-  .  129 

i 

19 

2.  250 

4 

500 

134 

755 

-4.000 

9.  152 

14.  742 

.106     1 

4 

593 

-  .  164 

20 

2.250 

4 

.750 

135 

773 

-6.000 

10.432 

14.7  4'3 

.  108     1 

4 

590 

-.181 

21 

2.250 

5 

.000 

135 

463 

-6.000 

11 .203 

14.743 

.108     1 

4 

591 

-  .175 

22 

2.  250 

5 

.250 

136 

525 

-7.000 

12.045 

14.743 

.  108     1 

4 

589 

-.195 

23 

2.250 

5 

.  500 

136 

.245 

-7.000 

12.418 

14.743 

.  109     1 

4 

.590 

-.  189 

* 

1 

2.500 

.000 

133 

.212 

-7 .000 

18.899 

14.  742 

.  107     1 

4 

.  596 

-.134 

2 

2.500 

.250 

134 

.589 

-7.000 

19. 150 

14. 743 

.  114     1 

4 

.  594 

-  .  153 

3 

2.500 

.500 

135 

.  737 

-8.000 

19. 506 

14.744 

.118     1 

4 

.  592 

-  .  171 

4 

2.500 

.750 

134 

.1B6 

-7.000 

20.348 

14.743 

.116     1 

4 

.595 

-  .144 

5 

2.500 

1 

.000 

132 

.  180 

-6.000 

21.764 

14.743 

.112     1 

4 

.599 

-.110 

6 

2.500 

1 

.250 

130 

.867 

-6.000 

23.229 

14.743 

.112     1 

4 

.602 

-  .086 

7 

2.500 

1 

.500 

129 

.322 

-3.000 

24. 708 

14.743 

.113     1 

4 

.605 

-  .057 

8 

2.500 

1 

.750 

128 

.325 

.000 

26.047 

14.743 

.112     1 

4 

.607 

-.040 

• 

117 


9 

2  .  500 

2  .000 

124  . 

946 

3  .000 

26. 810 

14.743 

.115     1 

61  5 

.023 

10 

2.500 

2  .  250 

121  . 

048 

0  .000 

2  6.057 

14.743 

.  115     1 

621 

.077 

11 

2.500 

2.500 

118. 

980 

12.000 

23.961 

14.743 

.115     1 

627 

.  125 

12 

2.500 

2.750 

118. 

270 

13.000 

20.871 

14.744 

.116     1 

62B 

.  138 

13 

2.500 

3.000 

119. 

701 

12.000 

16.567 

14.744 

.116     1 

625 

.  1  14 

14 

2.500 

3.250 

121. 

593 

10.000 

12.813 

14.743 

.113     1 

621 

.079 

15 

2.500 

3.500 

124. 

659 

5.000 

10. 150 

14.743 

.  114     1 

615 

.027 

16 

2.500 

3.7  50 

126. 

673 

2.000 

8.682 

14.743 

.115     1 

611 

-.007 

17 

2.  500 

4.000 

130. 

092 

-1  .000 

8.662 

14.744 

.116     1 

60<i 

-  .067 

18 

2.500 

4.250 

132. 

067 

-3.000 

9.495 

14.743 

.116     1 

600 

-  .  104 

19 

2.500 

4  .500 

133. 

770 

-4.000 

10.  173 

14.743 

.116     1 

596 

-  .136 

20 

2.500 

4  .  750 

134 

656 

-5.000 

10.  938 

14  .  743 

.  113     1 

59*. 

-.155 

21 

2.500 

5.000 

134 

945 

-0.000 

11.654 

14  .743 

.  112     1 

593 

-  .  162 

- 

22 

2.  500 

5.250 

135 

749 

-7.000 

12. 251 

14.743 

.  113     1 

591 

-  .176 

23 

2.  500 

5  .  500 

,13  5 

960 

-7 .000 

12.703 

14.743 

.114     1 

591 

-.179 

1 

2.750 

.000 

134 

168 

-8 .000 

18 . 668 

14.743 

.  10B     1 

594 

-.  151 

2 

2.750 

.250 

135 

286 

-8.000 

18.770 

14  .  743 

.111     1 

592 

-  .  169 

i 

3 

2.750 

.  500 

133 

362 

-7 .000 

19. 8  00 

14.743 

.  112     1 

596 

-  .  132 

k 

2.7  50 

.750 

133 

330 

-7.000 

20.206 

14.  7*43 

.112     1 

596 

-.  131 

5 

2.750 

1.000 

132 

317 

-6.000 

21  .  125 

14.743 

.113     1 

599 

-  .  Ill 

6 

2.750 

1  .  250 

128 

822 

-4 .000 

22.820 

14.743 

.113     1 

606 

-.048 

7 

2.750 

1  .500 

127 

.548 

-2.000 

23.608 

14.743 

.112     1 

609 

-  .026 

• 

8 

2.750 

1.750 

126 

.496 

.000 

24.278 

14.743 

.111     1 

611 

-  .008 

9 

2.750 

2.000 

124 

.958 

3.000 

2  4.495 

14.743 

.  114     1 

614 

.022 

10 

2.  750 

2.250 

121 

.745 

7.000 

24.064 

14.743 

.111     1 

621 

.074 

11 

2.750 

2.500 

120 

.189 

10.000 

22. 501 

14.743 

.114     1 

624 

.103 

12 

2.750 

2.750 

119 

.673 

10.000 

20. 362 

14.743 

.112     1 

625 

.  110 

13 

2.750 

3.000 

120 

.237 

9.000 

17.457 

14.743 

.114     1 

624 

.103 

14 

2.750 

3.250 

122 

.532 

7.000 

14.598 

14.743 

.112     1 

.619 

.062 

15 

2.750 

3.500 

124 

.683 

4.000 

12.629 

14.743 

.  109     1 

.614 

.022 

• 

118 


16 

2.7  50 

3.750 

127. 

4  32 

1  .000 

11 . 250 

14.743 

.111     1 

It . 

609 

-  .025 

17 

2.750 

4  .000 

129. 

310 

-1.000 

10. 685 

14.743 

.  Ill     1 

4. 

605 

-  .059 

18 

2.  750 

4  .  250 

130. 

544 

-2 .000 

10. 684 

14.743 

.  109     1 

4 . 

602 

-  .082 

19 

2.750 

4  .  500 

132. 

113 

-3.000 

11.232 

14.743 

.  112     1 

4 . 

599 

-.109 

20 

2.750 

4.750 

134. 

199 

-4.000 

11.598 

14  .743 

.  Ill     1 

4. 

594 

-.149 

21 

2.750 

5.000 

134. 

250 

-5.000 

12.084 

14.743 

.  Ill     1 

k . 

594 

-  .  150 

22 

2.750 

5.250 

134. 

819 

-6.000 

12.688 

14.743 

.109     1 

4 . 

593 

-  .  162 

23 

2.750 

5.  500 

135. 

451 

-6.000 

12.928 

14.743 

.113     1 

u  . 

592 

-  .  171 

1 

3.000 

.000 

133. 

856 

-6.000 

10.782 

14.74  3 

.114     1 

4 . 

595 

-.  140 

2 

3  .000 

.250 

131. 

943 

-5.000 

19.493 

14.743 

.111     1 

u . 

599 

-.107 

3 

3.000 

.  500 

132. 

3  24 

-5.000 

19.619 

14.743 

.  Ill     1 

4. 

598 

-  .  114 

k 

3  .000 

.  750 

132. 

137 

-6.000 

20.044 

14.743 

.  109     1 

4 

599 

-  .  112 

5 

3  .000 

1.000 

130 

158 

-4.000 

21.078 

14.743 

.  110     1 

4 

603 

-  .075 

6 

3.000 

1  .250 

128 

937 

-3.000 

21  . 836 

14  .74? 

.107     1 

4 

605 

-.055 

7 

3  .000 

1  .  500 

126 

693 

-1.000 

22.696 

14.  743 

.113     1 

4 

611 

-  .010 

• 

8 

3.000 

1.  750 

125 

627 

1  .000 

23.022 

14.743 

.  Ill     1 

4 

613 

.007 

9 

3  .000 

2.000 

123 

270 

3.000 

23.437 

14. 743 

.  108     1 

4 

617 

.046 

| 

10 

3  .000 

2  .250 

121 

979 

5.000 

22.902 

14. 743 

.  109     1 

4 

620 

.069 

11 

3.000 

2.  500 

120 

487 

7  .000 

21.763 

14.743 

.  Ill     1 

4 

623 

.095 

1 

12 

3.000 

2.750 

120 

773 

8.000 

19.838 

14.743 

.  Ill     1 

4 

.623 

.091 

13 

3.000 

3.000 

122 

215 

6.000 

17.752 

14.743 

.  112     1 

4 

.620 

.067 

14 

3.000 

3.  250 

124 

.387 

4.000 

15. 308 

14.743 

.113     1 

4 

.615 

.030 

. 

15 

3.000 

3.  500 

125 

.427 

2.000 

13.911 

14.743 

.  114     1 

4 

.613 

.014 

16 

3.000 

3.  750 

127 

.  281 

1.000 

12.724 

14. 743 

.  113     1 

4 

.609 

-  .020 

17 

3.000 

4  .000 

128 

.661 

.000 

12.058 

14.743 

.  109     1 

4 

.606 

-  .048 

18 

3.000 

4  .250 

130 

.582 

-2.000 

11.948 

14.743 

.  Ill     1 

4 

.602 

-  .082 

19 

3.000 

4  .  500 

132 

.  377 

-3.000 

12.096 

14.743 

.  110     1 

4 

.  598 

-.116 

20 

3.000 

4.750 

132 

.726 

-4.000 

12. 532 

14.743 

.111     1 

4 

.598 

-.  121 

21 

3.000 

5.000 

133 

.751 

-5.000 

12.591 

14.743 

.112     1 

4 

.595 

-  .  139 

22 

3.000 

5.  250 

133 

.585 

-5.000 

13.246 

14.743 

.  108     1 

4 

.595 

-  .  140 

23 

3.000 

5.  500 

133 

.965 

-6.000 

13.427 

,14.743 

.  110     1 

4 

.  595 

-  .  145 

119 


APPENDIX  H.  RESULT3C.DAT 


299 

1 

.000 

000 

96 

737 

-11 

000 

13. 

372 

14 

690 

.172     1 

4 

eiit 

.303 

2 

.000 

250 

97 

806 

-12 

000 

12. 

936 

14 

690 

.168     1 

k 

612 

.279 

3 

.000 

500 

99 

054 

-13 

000 

12. 

116 

14 

690 

.170     1 

4 

610 

.259 

It 

.000 

750 

99 

880 

-14 

000 

11 

094 

14 

690 

.  169     1 

k 

609 

.243 

5 

.000 

1 

000 

100 

846 

-16 

000 

10 

170 

14 

690 

.  169     1 

It 

607 

.225 

6 

.000 

1 

250 

100 

238 

-15 

000 

8 

460 

14 

681 

.065     1 

It 

599 

.133 

7 

.000 

1 

500 

98 

974 

-3 

000 

-4 

648 

14 

653 

-.250     1 

4 

573 

-.  160 

8 

.000 

1 

750 

98 

912 

15 

000 

-13 

092 

14 

633 

-.478     1 

it 

553 

-.386 

9 

.000 

2 

000 

94 

907 

37 

000 

-22 

132 

14 

602 

-  .837     1 

k 

528 

-  .  673 

10 

.000 

2 

250 

99 

870 

65 

000 

-13 

838 

14 

619 

-.639     1 

It 

538 

-  .  565 

11 

.000 

2 

500 

104 

963 

73 

000 

2 

199 

14 

640 

-.403     1 

It 

550 

-  .426 

12 

.000 

2 

750 

108 

612 

71 

000 

10 

560 

14 

636 

*  - .452     1 

It 

539 

-  .  547 

13 

.000 

3 

000 

103 

958 

59 

000 

15 

672 

14 

602 

-  .837     1 

It 

513 

-.840 

14 

.000 

3 

250 

102 

772 

30 

000 

19 

279 

14 

595 

-.907     1 

It 

509 

-  .  888 

15 

.000 

3 

500 

103 

466 

11 

000 

21 

003 

14 

625 

-.567     1 

It 

538 

-.  561 

1 

16 

.000 

3 

750 

101 

685 

-4 

000 

12 

623 

14 

,6  3  5 

-.457     1 

It 

551 

-  .417 

17 

.000 

4 

000 

107 

602 

-15 

000 

10 

328 

14 

676 

.011     1 

k 

582 

-.064 

18 

.000 

4 

250 

104 

468 

-16 

.000 

13 

372 

14 

683 

.083     1 

It 

593 

.070 

19 

.000 

4 

500 

103 

733 

-14 

000 

12 

835 

14 

688 

.  143     1 

k 

600 

.  144 

20 

.000 

4 

750 

101 

574 

-13 

000 

12 

751 

14 

687 

.  138     1 

It 

603 

.  180 

21 

.000 

5 

000 

99 

819 

-12 

000 

12 

832 

14 

687 

.  135     1 

4 

606 

.  210 

22 

.000 

5 

.250 

98 

659 

-11 

.000 

12 

899 

14 

687 

.  130     1 

4 

607 

.226 

23 

.000 

5 

.500 

96 

.949 

-10 

.000 

13 

331 

14 

.686 

.125     1 

It 

609 

.252 

1 

.250 

.000 

96 

.502 

-11 

.000 

14 

210 

14 

.688 

.  148     1 

It 

.612 

.284 

2 

.250 

.250 

98 

.  184 

-13 

.000 

13 

427 

14 

689 

.  157     1 

k 

.610 

.261 

3 

.250 

.500 

99 

.601 

-14 

.000 

32 

825 

14 

.690 

.161     1 

u 

.608 

.240 

it 

.250 

.750 

102 

.050 

-16 

.000 

11 

671 

14 

.690 

.  168     1 

4 

.605 

.201 

5 

.250 

1 

.000 

103 

.  829 

-17 

.000 

10 

458 

•  14 

.689 

.  161     1 

It 

.601 

.  160 

120 


6 

.  250 

1  .250 

104. 

395 

-17 .000 

8  .  573 

14.  682 

.074      1 

It 

593 

.063 

7 

.250 

1  .  500 

103. 

435 

-7 .000 

-  .  823 

14  .653 

-.249     1 

4 

566 

-  .  242 

8 

.  250 

1  .  750 

103. 

086 

-1.000 

-7  .  788 

14  .  624 

- . 507      1 

It 

537 

-  .  574 

9 

.250 

2.000 

94. 

559 

39.000 

-22. 304 

14.579 

-  1  .093     1 

4 

506 

-.923 

10 

.250 

2.250 

100. 

821 

68.000 

-13. 193 

14.601 

-.842     1 

4 

518 

-.786 

1] 

.250 

2.500 

110. 

243 

77.000 

.516 

14.640 

-.396     1 

4 

541 

-.524 

12 

.250 

2.750 

113. 

049 

74.000 

5.440 

14.633 

- .475     1 

4 

529 

-.662 

13 

.250 

3.000 

101. 

824 

61 .000 

3.994 

14.575 

-1.139     1 

4 

(.90 

-1. 102 

14 

.250 

3.250 

101. 

485 

26.000 

4.513 

14. 558 

-1.329     1 

4 

474 

-1.  285 

15 

.250 

3.500 

106. 

355 

6.000 

11.286 

14.606 

-.792     1 

4. 

513 

-.843 

16 

.250 

3  .750 

107. 

919 

-9.000 

8.  127 

14.639 

-.408     1 

4. 

544 

-.4  90 

17 

.  250 

4.000 

106. 

988 

-17 .000 

10. 581 

14.673 

-.023     1 

It 

580 

-.086 

18 

.250 

4.250 

104 

264 

-18.000 

12.  185 

14.685 

.  .110     1 

4 

596 

.101 

19 

.  250 

4.  500 

102 

121 

-16.000 

12.377 

14.687 

.  138     1 

It 

602 

.  170 

20 

.250 

4.750 

101 

333 

-14 .000 

12. Ill 

14  .  688 

.  141     1 

It 

604 

.  188 

21 

.250 

5  .000 

99 

732 

-12.000 

12. 258 

14.687 

.  137     1 

It 

606 

.214 

22 

.250 

5.  250 

97 

907 

-11.000 

12. 518 

14. 687 

.  132     1 

4 

609 

.  242 

23 

.250 

5.500 

96 

575 

-10.000 

12.897 

14.687 

.132     1 

4 

611 

.266 

1 

.500 

.000 

91 

565 

-] 1.000 

16. 192 

14.6,74 

-.010     1 

It 

606 

.  211 

2 

.  500 

.250 

97 

493 

-32.000 

14.454 

14.686 

.  126     1 

4 

609 

.  244 

3 

.500 

.500 

100 

384 

-14.000 

13.425 

14.688 

.  147     1 

It 

606 

.211 

4 

.  500 

.750 

103 

031 

-16.000 

12.52? 

14.688 

.  148     1 

4 

602 

.163 

5 

.500 

1  .000 

105 

.289 

-18.000 

11.701 

14.688 

.147     1 

It 

598 

.  118 

6 

.500 

1  .  250 

107 

.485 

-18.000 

10.302 

14.683 

.083     1 

It 

.588 

.010 

7 

.500 

1.  500 

108 

.  919 

-10.000 

4  .025 

14.661 

-.164     1 

It 

.  564 

-  .266 

8 

.500 

1.750 

99 

.871 

5.000 

-1 .323 

14.612 

-.719     1 

4 

.531 

-.645 

9 

.500 

2.000 

92 

.828 

33.000 

-15.353 

14. 558 

-1.331     1 

u 

.488 

-1. 131 

10 

.  500 

2.250 

105 

.376 

77.000 

-8.410 

14.605 

-.804     1 

It 

.514 

-  .835 

11 

.500 

2.500 

114 

.442 

78.000 

1.004 

14.651 

-.277     1 

It 

.  544 

-.493 

12 

.500 

2.750 

115 

.768 

75.000 

.447 

14.642 

-.375     1 

It 

.533 

-.619 

121 


13 

.500 

3  .000 

106. G14 

61.000 

-8.569 

14.  579 

-1.088     1 

It 

487 

-1  .  144 

Id 

.500 

3  .  250 

105.033 

23.000 

-14.855 

14.559 

-1 .317     1 

It 

469 

-1  .  341 

15 

.  500 

3.500 

110.667 

1.000 

-1.026 

14.609 

-.753     1 

4 

509 

-.890 

16 

.500 

3.750 

111. 552 

-12.000 

3.680 

14.649 

-.302     1 

It 

bkl 

-.458 

17 

.500 

4  .000 

107.804 

-18.000 

8.075 

14.672 

-.038     1 

It 

577 

-.  117 

18 

.500 

4.250 

105.384 

-18.000 

10.495 

14.685 

.110     1 

It 

594 

.079 

19 

.500 

4  .500 

100.400 

-16.000 

11.451 

14.687 

.130     1 

4 

604 

.  194 

20 

.  500 

4.750 

99. 509 

-14.000 

11.722 

14.687 

.127     1 

4 

606 

.207 

21 

.500 

5.000 

98.404 

-12.000 

11.801 

14.686  . 

.124     1 

4 

607 

.225 

22 

.500 

5.250 

97.393 

-11.000 

12. 154 

14.686 

.121     1 

It 

608 

.  241 

23 

.  500 

5.500 

96.005 

-10.000 

12.782 

14.686 

.118     1 

4 

610 

.262 

1 

.750 

.000 

97.198 

-12.000 

15. 176 

14.688 

.145     1 

k 

611 

.  268 

2 

.750 

.250 

98.757 

-13.000 

14.856 

14.689 

.151     1 

It 

609 

.245 

*" 

3 

.750 

.500 

100.395 

-14.000 

14.467 

14.689 

.154     1 

It 

607 

.219 

4 

.750 

.  750 

•103.847 

-15.000 

13.774 

14.689 

.156     1 

It 

601 

.  155 

5 

.750 

1.000 

106.623 

-16.000 

13.358 

14.689 

.151     1 

It 

596 

.095 

6 

.750 

1.250 

109. 165 

-19.000 

12.820 

14.683 

.082     1 

It 

585 

-  .024 

1 

7 

.750 

1.  500 

109.984 

-12.000 

9.  363 

14.663 

-.142     1 

It 

564 

-  .  265 

8 

.  750 

1  .750 

106.534 

.000 

7.824 

14i630 

-.512     1 

It 

537 

-  .  566 

9 

.750 

2.000 

9  7.337 

26.000 

1  .073 

14.581 

-1.075     1 

4 

.503 

-  .  954 

10 

.  750 

2.250 

103. 220 

61.000 

.406 

14.60  3 

-.817     1 

It 

.  516 

-  .  806 

11 

.750 

2.  500 

117 .028 

74 .000 

3.171 

14.662 

-.148     1 

It 

.  550 

-.  420 

12 

.750 

2.750 

117.409 

74.000 

-.  160 

14.658 

-.198     1 

It 

.545 

-  .478 

13 

.750 

3.000 

112.016 

56.000 

-15. 666 

14.605 

-.799     1 

It 

.  502 

-  .964 

14 

.750 

3.250 

111.196 

23.000 

-22. 328 

14.597 

-.884     1 

It 

.496 

-1.032 

15 

.750 

3.  500 

115.  114 

-1 .000 

-10. 267 

14.626 

-.564     1 

It 

.517 

-  .794 

16 

.750 

3.750 

112.534 

-12.000 

-1.891 

14.659 

-.184     1 

It 

.  556 

-.360 

17 

.750 

4.000 

109. 175 

-16.000 

4.855 

14.680 

.055     1 

4 

.  583 

-.052 

18 

.750 

4.250 

104. 321 

-15.000 

8  .493 

14.685 

.111     1 

It 

.  596 

.101 

19 

.750 

4  .  500 

99.889 

-13.000 

10. 106 

14.687 

.130     1 

It 

.  605 

.  204 

122 


20 

.  750 

4 

750 

98. 

034 

-12.000 

11  . 

0/9 

14  . 

687 

.134     1 

4 

609 

.  241 

21 

.  750 

5 

000 

97  . 

562 

-11.000 

11  . 

411 

14  . 

686 

.  126     1 

u 

609 

.242 

22 

.750 

5 

250 

97. 

200 

-10.000 

11. 

755 

14. 

687 

.  132     1 

4 

610 

.  255 

23 

.750 

5 

500 

96 

015 

-9.000 

12. 

364 

14. 

686 

.125     1 

4 

611 

.270 

1 

1.000 

000 

92. 

973 

-11.000 

16. 

732 

14. 

679 

.045     1 

u 

609 

.2<»3 

2 

1.000 

250 

97. 

219 

-12.000 

15. 

783 

14. 

686 

.  125     1 

It 

609 

.247 

3 

1.000 

500 

99. 

888 

-13.000 

15. 

442 

14. 

688 

.145     1 

4 

607 

.218 

It 

1  .000 

750 

102 

472 

-15.000 

15. 

370 

14. 

688 

.148     1 

4 

602 

.  173 

5 

1.000 

1 

000 

105 

913 

-16.000 

15. 

095 

14. 

688 

.141     1 

It. 

596 

.099 

6 

1  .000 

1 

250 

109 

237 

-18.000 

15. 

245 

14 

685 

.115     1 

4. 

588 

.007 

7 

1.000 

1 

500 

110 

946 

-14 .000 

14. 

219 

14 

675 

-.008     1 

4 . 

574 

-.  151 

8 

1  .000 

1 

750 

108 

157 

-3.000 

13 

826 

14 

646 

-  .  332     1 

It 

550 

-  .418 

9 

1  .000 

2 

000 

99 

929 

21 .000 

13 

533 

14 

612 

-.721     1 

1* . 

530 

-  .  648 

- 

10 

1.000 

2 

250 

101 

992 

51.000 

6 

347 

14 

616 

-.669     1 

Ik . 

531 

-  .635 

11 

1.000 

2 

500 

105 

643 

69.000 

6 

770 

14 

654 

-.239     1 

u. 

563 

-.275 

12 

1  .000 

2 

750 

111 

612 

71.000 

1 

544 

14 

671 

-.046     1 

It. 

569 

-  .203 

13 

1.000 

3 

000 

109 

155 

54.000 

-15 

730 

14 

.629 

-.521     1 

It 

532 

-  .627 

1 

14 

1.000 

3 

250 

107 

114 

22.000 

-25 

062 

14 

.619 

-  .635     1 

4 

526 

-.701 

i 

15 

1.000 

3 

500 

108 

876 

-1.000 

-18 

298 

14 

.637 

-.433     1 

It 

540 

-.  534 

16 

1  .000 

3 

.750 

108 

996 

-12.000 

-6 

883 

14 

.666 

-.101     1 

4 

569 

-.  204 

17 

1  .000 

4 

.000 

105 

.868 

-13.000 

2 

.44  7 

14 

.682 

.072     1 

It 

590 

.032 

18 

1.000 

4 

.250 

102 

.  374 

-14 .000 

6 

.782 

14 

.686 

.  117     1 

It 

600 

.  143 

• 

19 

1.000 

4 

.500 

98 

.241 

-12.000 

9 

.104 

14 

.686 

.  121     1 

It 

607 

.225 

20 

1.000 

4 

.750 

96 

.  145 

-11.000 

10 

.421 

14 

.686 

.  122     1 

It 

611 

.264 

21 

1.000 

5 

.000 

96 

.077 

-10.000 

11 

.098 

14 

.686 

.  123     1 

4 

.611 

.266 

22 

1.000 

5 

.250 

96 

.339 

-10.000 

11 

.624 

14 

.685 

.115     1 

It 

.  610 

.253 

23 

1.000 

5 

.500 

94 

.905 

-9.000 

12 

.  157 

14 

.685 

.114     1 

4 

.612 

.278 

1 

1.250 

.000 

96 

.335 

-11.000 

15 

.709 

14 

.686 

.  121     1 

4 

.610 

.  260 

2 

1.  250 

.250 

97 

.455 

-13.000 

15 

.844 

14 

.686 

.123     1 

It 

.609 

.24  2 

3 

1.  250 

.  500 

98 

.732 

-13.000 

16 

.025 

14 

.686 

.  126     1 

4 

.607 

.  221 

123 


4 

1.250 

.750 

101. 

681 

-14 .000 

15.953 

14.686 

.  125     1 

(. . 

602 

.  165 

5 

1  .250 

1  .000 

104. 

803 

-16.000 

16. 256 

14. 686 

.  121     1 

4 . 

596 

.  102 

6 

1.  250 

1.250 

107. 

173 

-17 .000 

17.238 

14.685 

.  114     1 

<4 . 

592 

.048 

7 

1.250 

1.500 

108. 

121 

-16.000 

18.968 

14.681 

.059     1 

u  . 

585 

-  .026 

8 

1.  250 

1.750 

101. 

937 

-7.000 

21. 148 

14.653 

-.252     1 

4 . 

568 

-.217 

9 

1.250 

2.000 

95. 

570 

3  2.000 

21. 238 

14.629 

-.531     1 

it . 

554 

-.379 

10 

1.250 

2.250 

93. 

882 

41.000 

15.812 

14.626 

-.563     1 

4 . 

554 

-.  381 

11 

1.250 

2.500 

100. 

797 

61.000 

11.438 

14.672 

-.041     1 

k 

589 

.015 

12 

1.250 

2.750 

102. 

990 

62.000 

4.491 

14.687  . 

.129     1 

4 

600 

.144 

13 

1.250 

3.000 

106. 

615 

50.000 

-12.738 

14.665 

-.119     1 

4 . 

572 

-.175 

14 

1.250 

3.250 

101. 

650 

20.000 

-28. 102 

14.640 

-.399     1 

u . 

556 

-  .  359 

15 

1  .250 

3  .  500 

105 

706 

.000 

-22.013 

14  .  663 

-.137     1 

4 

572 

-.  174 

16 

1.250 

3.  750 

105 

923 

-10.000 

-9.544 

14.677 

.023     1 

k . 

586 

-  .018 

*" 

17 

1.250 

4  .000 

102 

964 

-13.000 

.221 

14.683 

.092     1 

It 

597 

.  107 

18 

1.  250 

4.250 

.  99 

508 

-12.000 

5.  272 

14  .685 

.110     1 

4 

604 

.  191 

19 

1.250 

4.500 

96 

297 

-11.000 

7  .984 

14.685 

.  106     1 

'i 

609 

.24  4 

20 

1.  2  50 

4.750 

96 

357 

-11.000 

9.495 

14.685 

.  104     1 

4 

609 

.  242 

1 

21 

1.  250 

5  .000 

95 

758 

-10.000 

10. 397 

14 . 684 

.103     1 

4 

609 

.  251 

22 

1.250 

5.250 

96 

329 

-10.000 

10. 988 

14 .684 

.102     1 

4 

608 

.240 

23 

1.250 

5.  ' 

95 

217 

-10.000 

11  .729 

14.684 

.097     1 

4 

610 

.255 

1 

1.  500 

.000 

95 

.007 

-10.000 

16.130 

14.686 

.123     1 

4 

612 

.285 

2 

1.  500 

.250 

96 

.037 

-11.000 

16.911 

14.686 

.  123     1 

4 

611 

.  267 

* 

3 

1.500 

.500 

97 

.264 

-12.000 

17  .  ■  -0 

14.686 

.125     1 

4 

.609 

.247 

4 

1.500 

.750 

99 

.451 

-13.000 

17 . 358 

14.687 

.  128     1 

4 

.606 

.210 

5 

1.  500 

1.000 

100 

.871 

-14.000 

18.  126 

14.686 

.124     1 

4 

.603 

.  179 

6 

1.  500 

1.250 

102 

.806 

-14.000 

19.671 

14.686 

.121     1 

4 

.600 

.  140 

7 

1.500 

1.500 

102 

.248 

-13.000 

22.922 

14.684 

.094     1 

4 

.  598 

.124 

8 

1.500 

1  .750 

100 

.171 

-7.000 

26. 166 

14.675 

-.004     1 

4 

.  593 

.064 

9 

1  .500 

2.000 

93 

.062 

9.000 

28.769 

14.654 

-.247     1 

u 

.  583 

-  .051 

10 

1.  500 

2.250 

87 

.  502 

34.000 

26.085 

14.648 

-  .312     1 

4 

.  585 

-  .023 

► 

124 


11 

1.  500 

2.  500 

88. 

434 

51.000 

20. 

498 

14.673 

-.028     l 

It . 

609 

.246 

12 

1  .  500 

2.750 

87  . 

630 

52.000 

11  . 

521 

14  .683 

.091     i 

4 . 

621 

.  378 

13 

1.  500 

3.000 

92. 

393 

41 .000 

-6. 

264 

14  .  682 

.070     1 

k . 

612 

.  278 

14 

1.  500 

3  .  250 

95. 

325 

21.000 

-19. 

830 

14 . 672 

-.034     1 

4. 

598 

.122 

15 

1.500 

3.  500 

97. 

392 

4.000 

-17. 

915 

14.676 

.003     1 

4. 

598 

.  123 

16 

1.500 

3.750 

98. 

584 

-5.000 

-8. 

052 

14. 681 

.066     1 

4. 

602 

.  164 

17 

1.500 

4  .000 

97. 

081 

-8.000 

147 

14.684 

.099     1 

It 

607 

.  224 

18 

1.500 

4.250 

95 

480 

-9.000 

5. 

111 

14.685 

.  105     1 

4. 

610 

.258 

19 

1.500 

4.500 

94 

482 

-10.000 

7. 

935 

14.684  . 

.103     1 

it 

611 

.275 

20 

1.500 

4.750 

94 

502 

-9.000 

9 

457 

14.685 

.  104     1 

4 

611 

.275 

21 

1.  500 

5.000 

94 

006 

-9.000 

10 

302 

14.684 

.099     1 

It 

612 

.  278 

22 

1.  500 

5.  250 

94 

448 

-9.000 

10 

910 

14.684 

.098     1 

It 

611 

.  270 

23 

1.500 

5.500 

93 

488 

-8.000 

11 

756 

14.684 

.094     1 

4 

612 

.282 

*■ 

1 

]  .750 

.000 

92 

894 

-9.000 

17 

627 

14.683 

.088     1 

4 

612 

.286 

2 

1.  750 

.  250 

•  94 

517 

-10.000 

17 

520 

14.685 

.  112     1 

4 

612 

.283 

3 

1.750 

.500 

95 

581 

-11.000 

17 

955 

14.686 

.  116     1 

4 

611 

.  268 

it 

1.750 

.  750 

96 

595 

-11 .000 

18 

597 

14  .  68  5 

.115     1 

4 

609 

.  24  9 

i 

5 

1.750 

1.000 

97 

845 

-11.000 

19 

581 

14.  686 

.  125     1 

it 

608 

.236 

6 

1.750 

1.250 

97 

916 

-11 .000 

21 

410 

14.686 

.120     1 

4 

608 

.  230 

7 

1.750 

1.  500 

95 

953 

-10.000 

25 

.644 

14.685 

.  Ill     1 

4 

610 

.  256 

8 

1.750 

1.750 

94 

.386 

-6.000 

30 

.707 

14.682 

.074     1 

4 

609 

.247 

9 

1.750 

2.000 

90 

.  723 

8.000 

35 

.600 

14.673 

-.024     1 

it 

606 

.  212 

10 

1.750 

2.  250 

86 

.  236 

26.000 

32 

.641 

14.672 

-.042     1 

4 

.611 

.268 

11 

1.750 

2.500 

83 

.219 

39.000 

25 

.691 

14.680 

.054     1 

it 

.624 

.411 

12 

1.750 

2.750 

81 

.682 

40.000 

16 

.  561 

14.686 

.  121     1 

4 

.631 

.  502 

13 

1.750 

3.000 

83 

.031 

32.000 

3 

.477 

14.685 

.111     1 

It 

.629 

.471 

lit 

1.750 

3.250 

87 

.326 

18.000 

-7 

.288 

14.682 

.081     1 

It 

.620 

.  373 

15 

1.750 

3.  500 

91 

.390 

7.000 

-8 

.551 

14.683 

.086     1 

It 

.615 

.  310 

16 

1.750 

3.750 

92 

.676 

-1.000 

-3 

.705 

14.683 

.091     1 

It 

.613 

.294 

17 

1.750 

4  .000 

92 

.762 

-5.000 

1 

.497 

14. 684 

.099     1 

It 

.614 

.  300 

125 


18 

1.750 

4.250 

93. 135 

-6.000 

5.401 

14.684 

.101     1 

(..613 

.  295 

19 

1.750 

4.  500 

93. 124 

-8 .000 

7  .974 

14 .684 

.098     1 

4.  613 

.293 

20 

1.750 

4.750 

93.009 

-8.000 

9.260 

14.684 

.095     1 

4.613 

.292 

21 

1.750 

5.000 

93.389 

-8.000 

10.438 

14.683 

.092     1 

4.612 

.  282 

22 

1.750 

5.250 

93.368 

-8.000 

11. 150 

14.684 

.097     1 

4.613 

.288 

23 

1.750 

5.500 

93.181 

-8.000 

11.886 

14.683 

.092     1 

4  .612 

.286 

1 

2.000 

.000 

93.  140 

-8.000 

17 .710 

14.683 

.086     1 

4.  612 

.280 

2 

2.000 

.250 

94.026 

-8.000 

17.834 

14.685 

.110     1 

4.613 

.289 

3 

2.000 

.500 

94.754 

-9.000 

18.071 

14.685  • 

.114     1 

4.612 

.281 

4 

2.000 

.750 

94.910 

-9.000 

18.865 

14.685 

.114     1 

4.612 

.277 

5 

2.000 

1  .000 

95.039 

-9.000 

19.911 

14.605 

.112     1 

4.611 

.  274 

6 

2.000 

1.250 

93. 587 

-8 .000 

22.007 

14.685 

.  115     1 

4.  614 

.  302 

7 

2.000 

1  .  500 

91.878 

-7.000 

25.237 

14.685 

.108     1 

4.  616 

.  324 

8 

2.000 

1.750 

90.  137 

-2.000 

29.506 

14  .684 

.099     1 

4.618 

.  345 

9 

2  .000 

2.000 

•  86.717 

7.000 

34.  178 

14.682 

.077     1 

4.621 

.  378 

10 

2.000 

2.250 

83.605 

18.000 

33.  172 

14.681 

.068     1 

4  .  624 

.419 

11 

2  .000 

2.500 

80.851 

28.000 

26. 148 

14.683 

.085     1 

4.  629 

.478 

12 

2.000 

2.750 

80.075 

30.000 

18.652 

14.686 

.  123     1 

4.  634 

.527 

i 

13 

2.000 

3.000 

80.632 

24.000 

9.609 

14.684 

.097     1 

4.  631 

.493 

14 

2.000 

3.  250 

83.079 

15.000 

1.474 

14.683 

.087     1 

4  .  627 

.446 

15 

2.000 

3.500 

86.533 

9.000 

-1.162 

14.683 

.084     1 

4.  622 

.  389 

16 

2.000 

3.  750 

89.123 

1.000 

.64P 

14.683 

.084     1 

4.618 

.  347 

17 

2.000 

4.000 

90.442 

-3.000 

3.640 

14.683 

.086     1 

4.616 

.326 

18 

2.000 

4.250 

91.362 

-5.000 

6.351 

14.683 

.086     1 

4.615 

.  311 

19 

2.000 

4.500 

92.002 

-6.000 

8.271 

14.683 

.085     1 

4.614 

.299 

20 

2.000 

4.750 

92.287 

-7.000 

9.510 

14.682 

.079     1 

4.613 

.288 

21 

2.000 

5.000 

93.240 

-7.000 

10.354 

14.682 

.081     1 

4.611 

.27  3 

22 

2.000 

5.250 

92.851 

-7.000 

11.256 

14.682 

.077     1 

4.612 

.277 

23 

2  .000 

5.500 

92.639 

-7 .000 

11.871 

14.682 

.076     1 

4.  612 

.279 

1 

2.250 

.000 

92.403 

-8.000 

18.107 

J4.683 

.091     1 

4.614 

.298 

126 


2 

2.250 

250 

93.  692 

-8.000 

18. 184 

14  .686 

.  121     1 

4.  614 

.  306 

3 

2.250 

500 

93.746 

-8.000 

18.877 

14.686 

.121     1 

4  .  6 1  4 

.  305 

4 

2.250 

750 

93.966 

-8.000 

19.615 

14.686 

.122     1 

4.  614 

.  302 

5 

2.250 

1. 

000 

93. 590 

-7.000 

20. 713 

14.686 

.  121     1 

k.eiu 

.  308 

6 

2.250 

1. 

250 

91.805 

-6.000 

22.343 

14.686 

.125     1 

4.617 

.  343 

7 

2.250 

1. 

500 

90.291 

-4.000 

24.914 

14.686 

.122     1 

4.  619 

.  365 

8 

2.250 

1. 

750 

88.769 

-1.000 

27 . 282 

14.686 

.123     1 

4.622 

.  392 

9 

2.  250 

2 

000 

85.854 

1.000 

30.087 

14.685 

.105     1 

4.624 

.420 

10 

2.250 

2 

250 

82.768 

13.000 

28.943 

14.685  . 

.107     1 

4.629 

.471 

11 

2.  250 

2 

500 

80.886 

20.000 

25.  116 

14.685 

.  107     1 

4.  631 

.499 

12 

2.  250 

2 

750 

80.979 

21.000 

19. 520 

14. 687 

.137     1 

4.634 

.  528 

13 

2.250 

3 

000 

81 .587 

18.000 

13.444 

14 .686 

.  120     1 

4.  631 

.  502 

14 

2  .250 

3 

250 

82.561 

13.000 

7  .897 

14 .684 

.  100     1 

4.628 

.467 

15 

2.  250 

3 

500 

85.327 

6.000 

4.955 

14.684 

.098     1 

4.624 

.  422 

16 

2.250 

3 

750 

•  88.056 

2.000 

4  .944 

14.684 

.  101     1 

4.621 

.381 

17 

2.  250 

4 

000 

90.051 

-2.000 

6.074 

14. 684 

.098     1 

4.  618 

.  345 

18 

2  .  250 

4 

250 

91.463 

.000 

7  .  592 

14.684 

.099     1 

4.616 

.  323 

i 

19 

2.  250 

4 

.500 

92.  141 

-5.000 

8.893 

14  .684 

.098     1 

4.615 

.310 

20 

2.  250 

4 

.750 

92. 377 

-6.000 

10.016 

14.68  4 

.096     1 

4.  614 

.  303 

21 

2.250 

5 

.000 

92.780 

-6.000 

10. 764 

14. 684 

.094     1 

4.613 

.295 

22 

2.  250 

5 

.  250 

92.515 

-6.000 

11.521 

14  .684 

.097     1 

4.  614 

.  303 

23 

2  .  250 

5 

.  500 

92.687 

-7.000 

12. 127 

14.684 

.097     1 

4.614 

.  300 

1 

2  .  500 

.000 

90.  377 

-8.000 

18 . 355 

14.679 

.045     1 

4  .  612 

.  286 

2 

2.  500 

.250 

92.261 

-8.000 

18. 384 

14.684 

.094     1 

4  .  614 

.303 

3 

2.  500 

.500 

92.515 

-8.000 

18.957 

14.684 

.099     1 

4.614 

.  304 

4 

2.500 

.750 

92.688 

-7.000 

19.447 

14.684 

.099     1 

4.  614 

.  302 

5 

2.500 

1 

.000 

91.432 

-6.000 

20.727 

14.684 

.100     1 

4.616 

.  324 

6 

2.500 

1 

.  250 

90.060 

-5.000 

22. 253 

14.684 

.096     1 

4.617 

.343 

7 

2.  500 

1 

.500 

88.704 

-3.000 

23.742 

14.684 

.097     1 

4.620 

.  367 

8 

2.500 

1 

.750 

86.956 

.000 

25.346 

•14.683 

.090     1 

4.621 

.388 

127 


9 

4.  624 

10 

4  .627 

11 

'..629 

12 

4  .  630 

13 

4.628 

14 

(..626 

15 

4.  623 

16 

4.620 

17 

4.618 

18 

4.616 

19 

4.614 

20 

4.613 

21 

4.613 

22 

4.  613 

23 

4.613 

1 

4  .  615 

2 

4.616 

3 

4  .  616 

4 

4.  618 

5 

4.617 

6 

4.  619 

7 

4.621 

8 

4.623 

9 

4  .  625 

10 

4.  628 

11 

4.629 

12 

4.630 

13 

4.628 

14 

4  .625 

15 

4.624 


500 

2. 

000 

85. 

017 

3  .000 

26. 

416 

14  . 

683 

414 

500 

2  . 

250 

82  . 

750 

8  .000 

25. 

923 

14. 

683 

448 

500 

2. 

500 

80. 

820 

12.000 

23. 

480 

14. 

683 

479 

500 

2. 

750 

81. 

486 

14.000 

19. 

431 

14. 

685 

490 

500 

3. 

000 

81. 

651 

12.000 

15. 

271 

14. 

683 

468 

500 

3. 

250 

82. 

508 

10.000 

11. 

210 

14. 

682 

440 

500 

3. 

500 

85. 

286 

5.000 

8. 

610 

14. 

683 

407 

500 

3. 

750 

87  . 

315 

1.000 

7  . 

650 

14. 

683 

376 

500 

4. 

000 

89. 

155 

-1.000 

8. 

057 

14. 

603 

344 

500 

4. 

250 

90. 

679 

-3.000 

8. 

715 

14. 

683 

322 

500 

4. 

500 

91. 

148 

-5.000 

9. 

540 

14 

682 

307 

500 

4 

750 

91 

837 

-6.000 

10 

196 

14 

682 

295 

500 

5 

000 

92 

040 

-7.000 

11 

007 

14 

682 

290 

500 

5 

250 

92 

076 

-7.000 

11 

834 

14 

682 

289 

500 

5 

500 

.  92 

143 

-8.000 

12 

116 

14 

682 

288 

750 

000 

91 

217 

-7 .000 

18 

257 

14 

683 

310 

750 

250 

91 

071 

-6.000 

IB 

570 

14 

.684 

324 

750 

.500 

91 

.242 

-6.000 

19 

.077 

14 

.  684 

324 

750 

.750 

90 

.064 

-6.000 

20 

.103 

14 

.  684 

343 

750 

1 

.000 

90 

.  326 

-6.000 

20 

.643 

14 

.684 

338 

750 

1 

.  250 

88 

.862 

-4.000 

21 

.  857 

14 

.684 

363 

750 

1 

.500 

87 

.337 

-2.000 

22 

.948 

14 

.684 

386 

• 

750 

1 

.750 

86 

.183 

.000 

23 

.884 

14 

.683 

401 

750 

2 

.000 

84 

.029 

3.000 

24 

.597 

14 

.683 

431 

750 

2 

.250 

82 

.292 

6.000 

24 

.046 

14 

.683 

457 

750 

2 

.500 

81 

.353 

9.000 

22 

.450 

14 

.683 

473 

750 

2 

.750 

81 

.458 

10.000 

19 

.475 

14 

.684 

.487 

.750 

3 

.000 

82 

.045 

10.000 

16 

.387 

14 

.683 

.460 

.750 

3 

.250 

83 

.289 

7.000 

13 

.230 

14 

.682 

.428 

.750 

3 

.500 

84 

.357 

4  .000 

11 

.507 

14 

.682 

.412 

• 

.086 

.084 

.086 

.  106 

.087 

.072 

.082 

.083 

.082 

.086 

.079 

.078 

.077 

.076 

.077 

.083 

.094 

i 

.097 

.097 

.095 

.097 

.094 

.090 

.087 

.086 

.088 

.103 

.086 

.072 

.073 


128 


16 

2.750 

3.750 

86. 106 

1.000 

10.404 

14.682 

.07  2     1 

621 

.  304 

17 

2.  750 

4  .000 

86.862 

.000 

10.048 

14.681 

.069     1 

620 

.  368 

18 

2.  750 

4  .250 

88 . 866 

-3.000 

10. 242 

14.  681 

.066     1 

617 

.  333 

19 

2.  750 

4.500 

89.  571 

-4.000 

10. 572 

14.682 

.072     1 

616 

.3  27 

20 

2.7  50 

4.750 

90. 345 

-5.000 

11 .021 

14 .681 

.068      1 

615 

.  310 

21 

2.750 

5.000 

90.  176 

-6.000 

11  .  803 

14.681 

.069     1 

615 

.314 

22 

2.7  50 

5.250 

90.795 

-6.000 

12. 137 

14. 681 

.064     1 

614 

.  299 

23 

2.  750 

5.500 

90.414 

-6.000 

12. 509 

14.681 

.064     1 

614 

.  305 

1 

3.000 

.000 

90.729 

-6.000 

18.408 

14.682 

.080     1 

615 

.316 

2 

3  .000 

.250 

91.000 

-6.000 

18.565 

14.683 

.087     1 

615 

.318 

3 

3.000 

.500 

90.7  94 

-5.000 

19.062 

14.683 

.086     1 

615 

.  320 

4 

3.000 

.750 

89.922 

-5.000 

19.655 

14.683 

.086     1 

617 

.  335 

5 

3.000 

1.000 

89.061 

-4 .000 

20.550 

14.683 

.087     1 

618 

.  350 

6 

3.000 

1.250 

88.460 

-3 .000 

21.208 

14.683 

*  .088     1 

619 

.  362 

7 

3  .000 

1  .  500 

87 .082 

-1.000 

22.031 

14.603 

.083     1 

621 

.37  9 

8 

3  .000 

1  .750 

85. 296 

1  .000 

22.807 

14 . 682 

.081     1 

623 

.405 

9 

3.000 

2.000 

84 .083 

3  .000 

22 . 941 

14.682 

.079     1 

625 

.423 

10 

3  .000 

2.  250 

82.996 

5.000 

22.115 

14 . 682 

.079'     1 

626 

.  439 

11 

3.000 

2.  500 

82.411 

7  .000 

21.  137 

14.682 

.081     1 

627 

.450 

i 

12 

3  .000 

2.750 

02.629 

8.000 

19.  176 

14.683 

.084      1 

627 

.450 

13 

3.000 

3.000 

83. 322 

6.000 

16.658 

14.682 

.074     1 

625 

.429 

14 

3  .000 

3  .  250 

83.675 

4.000 

14.871 

14.681 

.065     1 

624 

.415 

15 

3.000 

3.  500 

85.012 

3  .000 

13.020 

14.681 

.061     1 

622 

.  390 

16 

3.000 

3.750 

86.292 

.000 

11.938 

14.681 

.063     1 

620 

.371 

17 

3.000 

4  .000 

87.459 

-1.000 

11  .491 

14.681 

.062     1 

618 

.  352 

18 

3.000 

4  .  250 

88.465 

-2 .000 

11.467 

14.681 

.060     1 

617 

.333 

19 

3.000 

4.500 

88.831 

-3.000 

11.636 

14.681 

.061     1 

616 

.328 

20 

3.000 

4.750 

89.728 

-5.000 

11.852 

14.601 

.062     1 

615 

.315 

21 

3.000 

5.000 

89.913 

-6.000 

12.  185 

14.681 

.059     1 

614 

.309 

22 

3.000 

5.  250 

90.264 

-6.000 

12.537 

14.680 

.055     1 

614 

.299 

23 

3.000 

5.  500 

89.872 

-6.000 

12.973 

.14.680 

.054     1 

614 

.  304 

129 


LIST  OF  REFERENCES 


1.  Roane,  D.  P.,  The  Effect  of  a  Turbulent  Airstream  on  a  Vertically- 
Launched  Missile  at  High  Angles  of  Attack,  M.S.  Thesis,  Naval 
Postgraduate  School,  Monterey,  CA,  December  1987. 

2.  Rabang,  M.  P.,  Turbulence  Effects  on  the  High  Angle  of  Attack 
Aerodynamics  of  a  Vertically-Launched  Missile,  M.S.  Thesis,  Naval 
Postgraduate  School,  Monterey,  CA,  June  1988. 

3.  Lung,  Ming-Hung,  Flowfield  Measurements  in  the  Vortex  Wake  of  a 
Missile  at  High  Angle  of  Attack  in  Turbulence,  M.S.  Thesis,  Naval 
Postgraduate  School,  Monterey,  CA,  December  1988. 

4.  Dahlem,  V.,  Flaherty,  J.  I.,  Shereda,  D.  E.,  High  Angle  of  Attack  Missile 
Aerodynamics  at  Mach  Numbers  0.3  to  1.5,  Technical  Report  AFWAL- 
TR-80-3070,  Rutgers,  NJ,  November  1980. 

5.  Ericsson,  L.  E.  and  Reding,  J.  P.,  Asymmetric  Vortex  Shedding  From 
Bodies  of  Revolution,  Tactical  Missile  Aerodynamics,  American  Institute 
of  Aeronautics  and  Astronautics,  Inc.,  1986. 

6.  Gregoriou,  G.,  Modern  Missile  Design  for  High  Angle  of  Attack, 
AGARD/VKI  lecture  series  no.  121,  High  Angle  of  Attack  Aerodynamics, 
March  1982. 

7.  Keener,  E.  R.  and  Chapman,  G.  T.,  Similarity  in  Vortex  Asymmetries 
Over  Slender  Bodies  and  Wings,  AIAA  Journal,  v.  15,  no.  9.,  pp.  1370- 
1372,  September  1977. 

8.  Deffenbaugh,  F.  D.,  and  Koerner,  Asymmetric  Vortex  Wake  Development 
on  Missiles  at  High  Angles  of  Attack,  Journal  of  Spacecraft,  jv.  14,  no.  3, 
pp.  155-162,  March  1977. 

9.  Clark,  W.  C,  and  Nelson,  R.  C,  Body  Vortex  Formation  on  Missiles  at 
High  Angles  of  Attack,  AIAA  Paper  76-65,  January  1976. 

10.  Achenbach,  Elmar,  Influence  of  Surface  Roughness  on  the  Crossflow 
Around  a  Circular  Cylinder,  Journal  of  Fluid  Mechanics,  pp.  321-335, 
v.  46,  part  2,  1971. 

11.  Ericsson,  L.  E.  and  Reding,  J.  P.,  Dynamics  ofForebody  Flow  Separation 
and  Associated  Vortices,  AIAA  Paper  83-2118,  August  1983. 


130 


12.  Reding,  J.  P.  and  Ericsson,  L.  E.,  Re-examination  of  the  Maximum 
Normalized  Vortex-Induced  Side  Force,  Journal  of  Spacecraft,  v.  21, 
no.  5,  September-October  1984. 

13.  Reding,  J.  P.  and  Ericsson,  L.  E.,  Maximum  Vortex-Induced  Side  Forces 
on  Slender  Bodies,  AIAA  Paper  77-1155,  1977. 

14.  Lamont,  P.  J.,  Pressure  Measurements  on  an  Ogive-Cylinder  at  High 
Angles  of  Attack  With  Laminar,  Transitional,  or  Turbulent  Separation, 
AIAA  Paper  80-1556,  1980. 

15.  Jorgensen,  L.  H.  and  Nelson,  E.  R.,  Experimental  Aerodynamic 
Characteristic  for  a  Cylindrical  Body  Revolution  with  Various  Noses 

atAOA  from  0°to  58°  and  Mach  Numbers  from  0.6  to  2.0,  NASA  TM  X- 
3130,  March  1975. 

16.  Wardlaw,  A.  B.,  Jr.,  and  Morrison,  A.  M.,  Induced  Side  Forces  on  Bodies 
of  Revolution  at  High  Angles  of  Attack,  NSWC/WOL/TR  75-176, 
November  1975. 

17.  Pick,  George  S.,  Investigation  of  Side  Forces  on  Ogive-Cylinder  Bodies  at 
High  Angles  of  Attack  in  the  M=0.5  to  1.1  Range,  AIAA  Paper  71-570, 
1971. 

18.  Keener,  E.  R.  and  others,  Side  Force  on  Forebodies  at  High  AOA  and 
Mach  Number  from  0.1  to  0.7;  Two  Tangent  Ogives,  Paraboloid  and 
Cone,  NASA-TM-X-3438,  Moffet  Field,  CA,  February  1977. 

19.  Kruse,  Robert  L,  Keener,  Earl  R.,  and  Chapman,  Gary  T.,  Investigation 
of  the  Asymmetric  Aerodynamic  Characteristics  of  Cylindrical  Bodies  of 
Revolution  With  Various  Variations  in  Nose  Geometry  and  Rotational 
Orientation,  NASA-TM  78533,  September  1979. 

20.  Yongnian,  Y.,  Xinzhi,  Y.,  and  Jianying,  L.,  Active  Control  of  Asymmetric 
Forces  at  High  Incidence,  Journal  of  Aircraft,  v.  25,  no.  2,  pp.  190-192, 
February  1988. 

21.  Keener,  E.  R.  and  others,  Side  Forces  on  a  Tangent  Ogive  Forebody  with 
a  Fineness  Ratio  of  3.5  at  High  AOA  and  Mach  From  0.1  to  0.7,  NASA- 
TM-X-3437,  Moffet  Field,  CA,  February  1977. 

22.  Bradshaw,  P.,  An  Introduction  to  Turbulence  and  Its  Measurement, 
Pergamon  Press,  1971. 


131 


23.  Tieleman,  H.  W.,  A  Survey  of  the  Turbulence  in  the  Marine  Surface  Layer 
for  the  Operation  of  Low -Reynolds  Number  Aircraft,  Virginia  Polytechnic 
Institute  Report  VPI-E-85- 10,  Blacksburg,  VA,  March  1985. 

24.  Castro,  I.  P.,  Effects  of  Free  Stream  Turbulence  on  Low  Reynolds  Number 
Boundary  Layers,  Journal  of  Fluids  Engineering,  Volume  106,  pp.  298- 
306,  September  1984. 

25.  Meier,  H.  U.  and  Kreplin,  H.  P.,  Influence  of  Freest  re  am  Turbulence  on 
Boundary-Layer  Development,  AIAA  Journal,  Volume  18,  Number  1, 
pp.  11-15,  January  1980. 

26.  Deane,  J.  R.,  Missile  Body  Vortices  and  Their  Interaction  With  Lifting 
Surfaces,  AGARD/VKI,  lecture  series,  no.  121,  High  Angle  of  Attack 
Aerodynamics,  March  1982. 

27.  Arya,  S.  P.,  Atmospheric  Boundary  Layers  Over  Homogeneous  Terrain, 
Engineering  Meteorology,  Elsevier  Scientific  Publishing  Company,  1982. 

28.  Healey,  J.  V.,  Simulating  the  Helicopter-Ship  Interface  as  an  Alternative 
to  Current  Methods  of  Determining  the  Safe  Operating  Envelopes,  Naval 
Postgraduate  School  Report,  NPS  67-86-003,  Monterey,  CA,  September 
1986. 

29.  Hancock,  P.  E.  and  Bradshaw,  P.,  The  Effect  of  Free -Stream  Turbulence 
on  Turbulent  Boundary  Layers,  Journals  of  Fluids  Engineering,  Vol.  105, 
p.  285,  September  1983. 

30.  Gregoriou,  G.  and  Knoche,  H.  G.,  High  Incidence  Aerodynamics  of 
Missiles  During  Launch  Phase,  MBB  GMBH  Report  UA-523-80,  January 
1980. 

31.  Department  of  Aeronautics,  Laboratory  Manual  for  Low-Speed  Wind 
Tunnel  Testing,  Naval  Postgraduate  School,  Monterey,  CA,  1983. 

32.  Velmex,  Inc.,  User's  Guide  to  8300  Series  Stepping  Motor 
Controller  I  Drivers,  East  Bloomfield,  NY,  January  1988. 

33.  United  Sensors,  Inc.,  5-Hole  Probes  Calibration  Manual,  Watertown, 
MA,  June  1988. 

34.  Hewlett-Packard,  Inc.,  PC  Instruments  System  Owner's  Guide  Using  HP 
610618  System  Interface,  February  1986. 

35.  Kindelspire,  D.  W.,  The  Effects  of  Freestream  Turbulence  on  Airfoil 
Boundary  Layer  Behavior  at  Low  Reynolds  Number,  Master's  Thesis, 
Naval  Postgraduate  School,  Monterey,  CA,  September  1988. 


132 


INITIAL  DISTRIBUTION  LIST 

No.  Copies 

1.  Defense  Technical  Information  Center  2 
Cameron  Station 

Alexanderia,  VA  22304-6145 

2.  Library,  Code  0142  2 
Naval  Postgraduate  School 

Monterey,  CA  93943-5002 

3.  Chairman  1 
Department  of  Aeronautics,  Code  67 

Naval  Postgraduate  School 
Monterey,  CA  93943-5000 

4.  Commander  1 
Naval  Surface  Warfare  Center,  Code  G205 

Dahlgren,VA  22448-5000 

5.  Standard  Missile  Program  Office  2 
PMS  422,  G205 

ATTN:      Thomas  McCants 

Dr.  Jesse  East 
Naval  Surface  Warfare  Center 
Dahlgren,VA  22448-5000 

6.  Commander  1 
Naval  Weapons  Center,  Code  406 

Chi   a  Lake,  CA  93555 

7.  Commander  1 
Naval  Surface  Weapons  Center 

Silver  Spring,  MD  20903-5000 

8.  Commander  1 
Pacific  Missile  Test  Center 

Point  Mugu,  CA  93041 

9.  NASA  Ames  Research  Center  1 
Technical  Library 

Moffett  Field,  CA  94035 


133 


10.  Prof.  R.  M.  Howard 

Department  of  Aeronautics,  Code  67Ho 
Naval  Postgraduate  School 
Monterey,  CA  93943-5000 

11.  Prof.  J.  V.  Healey 

Department  of  Aeronautics,  Code  67Ho 
Naval  Postgraduate  School 
Monterey,  CA  93943-5000 

12.  LT.  Lung,  Ming-Hung 

No.  114-5  Chung-Ching  Rd. 
Taichung  Taiwan  400 
Republic  of  China 

13.  LT.  J.  J.  Viniotis,  USN 
379  Ocean  Avenue 
Massapequa  Pk,  NY  11762 


134 


The 


^6399 


Viniotis 

pIowfieid 
launch  on 


V 


effects  of 


laui*hed  -HaJaef0111^- 


Thesis 

V6&99 

c.l 


Viniotis 

flowfield  effects  of 
launch  on  a  vertically- 
launched  missile. 


&7r»*5. 


